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Foreword 


Tis book is based on a series of lectures delivered to 
experimental physicists by one of the authors (L. Landau) in Moscow 
in 1954. 

In maintaining the lecture form in the printed edition we are empha- 
sizing the fact that the presentation makes no pretense at completeness 
and that the choice of subject matter is purely arbitrary. 

Since there is, at the present time, no rational theory of nuclear forces, 
we have limited our conclusions concerning nuclear structure to those 
which can be reached from an analysis of the available experimental 
data, using only general quantum-mechanical relations. 

No attempt has been made to give a bibliography of the literature; 
rather we have indicated only new experimental results. (A rather 
complete list of references is given in Blatt and Weisskopf, Theoretical 
Nuclear Physics, Foreign Literature Press, 1954).* 


L. LANDAU AND YA. SMORODINSKY 


* (Blatt and Weisskopf, Theoretical Nuclear Physics, John Wiley and Sons, Inc., 
New York, N.Y., 1952.] 
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LECTURE ONE 


Nuclear Forces 


Win the discovery of the neutron it became clear that the 
structural units of which the atomic nucleus is composed are two kinds of parti- 
cles: protons and neutrons, usually called nucleons. 

Of these particles only the proton is stable. The second particle—the neu- 
tron—is unstable. In the free state the neutron is radioactive; it is transformed 
into a proton, emitting an electron and a neutrino 


n — p+ eT + . 


The half-life of the neutron is approximately 10 min. 

It can be shown that the neutron is basically not an “elementary” particle and 
that it is more appropriately described, in contrast with the proton, as a com- 
pound particle, since the neutron can decay into “simpler” particles. On the 
other hand, although in the free state the proton is a stable particle, in the bound 
state (inside the nucleus) it can decay into a neutron, positron and neutrino: 


p—>n +et +r. 


On this basis we may consider the proton as a complex particle which is trans- 
formed into the “simpler” neutron. 

Essentially the foregoing means that the question of which of the two particles 
is the more elementary is not physically meaningful. Both particles are of 
equal rank; the question of which particle is more susceptible to decay depends 
only on energy considerations. The free neutron is heavier than the free proton, 
hence only one of the reactions—neutron decay—can be realized for free nucleons. 
Within the nucleus both reactions are possible and the nature of the decay is 
determined by the mass of the decaying nucleus and the masses of the nuclei 
which are possible decay products. 

These properties of nucleons mean that it is valid to consider both as ele- 
mentary particles which can be transformed into each other. In those phe- 
nomena, however, in which 8-decay of the nuclei does not take place (because of 
its small probability), we can ignore the neutron-proton transformation. 

We shall start with the characteristics of the proton and neutron. The basic 
difference between these particles lies in the fact that the neutron, as indicated 
by its name, is a neutral particle whereas the proton has a positive electric 
charge, equal in magnitude to that of the electron (4.80-10—Y cgs units). Both 
of these particles have the same spin, 1/2 (in units of %, in which this quantity is 
usually measured). The masses of both particles are almost the same: the 
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mass of the proton is 1836 electron masses (1.6:10~?4 g), and the neutron is 2.5 
electron masses (about 1.3 Mev) heavier. In the decay of the neutron 0.51 Mev 
goes into the formation of the electron and 0.78 Mev into the kinetic energy of 
the particles which are formed. 

This small difference in the masses (less than 0.2%) and the fact that the spins 
are the same are more than coincidental. If one disregards the difference in 
electric properties the proton and neutron are very much the same, and this 
similarity is of fundamental significance in the physics of the nucleus. 

This similarity is the basis for many nuclear phenomena and will not be 
labored. The effect of the similarity is most vividly seen in so-called “mirror 
nuclei.” 

If we consider two nuclei, the second of which is the same as the first with all 
protons replaced by neutrons and all neutrons replaced by protons, we have a 
pair of nuclei which are known as mirror nuclei. Since the number of particles 
is the same in both nuclei, the atomic weights of both nuclei are also almost the 
same. If the first nucleus has an atomic number Z, obviously the second has an 
atomic number A — Z (A is the mass number). 

Mirror nuclei of this type are well known in the light mass region (it is apparent 
that the heavier nucleus must be unstable). The first such pair is comprised of 
the neutron and proton themselves. Examples of subsequent pairs are He’, H?; 
Be’, Li’; B®, Be®; C14, O14; and so forth. It has been found experimentally 
that both mirror nuclei have very similar properties—almost identical binding 
energies, similar excitation spectra for the excited levels, identical spins. This 
symmetry in the properties of mirror nuclei is undoubtedly due to the properties 
of the interaction between protons and neutrons, i. e., the symmetry properties 
of nuclear forces. 

From the foregoing, it is seen that, aside from the relatively weak electric 
forces which act between two protons, all the available experimental facts 
indicate that the forces between two protons are very similar to the forces 
between two neutrons, a feature which is known as the charge symmetry of nuclear 
forces. As we shall see later, the charge symmetry of nuclear forces is actually 
the manifestation of a still more fundamental characteristic, the so-called 
tsotopic invariance. 

Contemporary nuclear physics still has no theory of nuclear forces. Hence, 
the only presently available sources of knowledge on nuclear forces are the 
experimental results. Thus, in formulating the properties of nuclear forces it is 
natural to rely on the experiments from which this knowledge is obtained. 

The simplest nuclear system is a system composed of a proton and a neutron 
and the most elementary experiments are those involving scattering of neutrons 
by protons, and deuteron experiments. The very existence of the deuteron—a 
compound particle consisting of a proton bound to a neutron (with a binding 
energy of 2.23 Mev)—indicates that the nuclear interaction is rather strong and 
that it is, to some degree at Icast, an essential factor in the attraction between 
these particles. The latter statement follows because the particles could not 
come together to form a stable compound particle if only repulsive forces were 
acting. 
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It does not follow, however, that this interaction is responsible for the poten- 
tial energy, 2.23 Mev, as would be the case for two particles at rest in classical 
mechanics. The classical picture of two bound particles at rest cannot be 
carried over into quantum mechanics. Actually, according to the very signifi- 
cance of the words “‘bound system,” in such a system the particles are close to each 
other, and consequently the uncertainty in the mutual position Ar is small; 
whence, from the laws of quantum mechanics it follows that the momentum 
cannot be small, and must be at least of order#/Ar. In particular, this situation 
means that the binding energy is not determined exclusively by the interaction 
energy but also by the effective distances over which the nuclear forces must act. 


Vr) 


= V, 


Fig. 1. The potential as a function of the distance between particles in the 
“square-well’’ model. 


In order to investigate the quantum-mechanical properties of the interaction, 
we consider a simple model for the proton-neutron interaction—the so-called 
“square-well potential.” In particular, we assume that the potential due to the 
nuclear forces is zero at distances r > Ry, where Ry is some constant; at smaller 
distances it is assumed that the potential is constant and equal in absolute value 
to Vo (since the nuclear forces are attractive forces this potential is negative). 
A curve showing the potential V(r) is given in Fig. 1, from which the origin of the 
designation “square-well”” becomes clear. 

Recalling that the force is the derivative of the potential, we see that the 
square-well potential gives rise to a configuration of two particles which interact 
only when the distance between them is exactly Ro, since the potential V(r) is 
constant forr < Ro; at the point r = Ro, however, the interaction force becomes 
infinite. It is obvious that such a potential is not satisfactory from a physical 
point of view. However, this simplified model is useful for making calculations 
and gives a fairly good picture of the basic properties of the quantum-mechanical 
interaction; for these reasons it is frequently used for illustrative purposes. 

First of all we shall consider the effect of increasing the depth V, keeping the 
radius of the well equal to Ro. Atsmall values of V it turns out that the system 
can have no stable states. In order for such a state to exist, that is, in order for a 
system to be formed, the depth of the well must satisfy the inequality 
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where u is the reduced mass of the system: 1/4 = (1/mı) + (1/me); if m = ma, 
for example in a system of two nucleons, a = m/2. This inequality (aside from 
the numerical coefficient) follows immediately if one notes that the momentum 
is of order #/R, which, as has already been mentioned, is to be associated with a 
particle in a potential well of radius R, and leads to a kinetic energy of order 
h?/pR?. Roughly speaking, this inequality leads to the requirement that the 
depth of the well must be greater than the kinetic energy of the particle. 

When the depth of the well exceeds (#?/8) (A?/uR?), the system has a single 
level. With a further increase in the depth of the well the energy of the level 
increases in absolute value. However, the level energy is always considerably 
smaller than the depth of the well. 
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Fig. 2. Energy of the level for a square-well potential. The radius of the 
well is taken as 2.80 - 10713 cm. 


The level energy remains relatively small as the depth of the well increases to 
high values. This situation 1s illustrated in Fig. 2, which shows the energy of 
the level as a function of the well depth. It is obvious from this figure that care 
must be used in estimating the strength of the interaction on the basis of binding 
energy data because the interaction energy can exceed the binding energy by a 
large margin. This effect is well illustrated in the deuteron. The bound state 
of the deuteron lies close to the top of the potential well so that the binding 
energy of 2.23 Mev is considerably smaller than the interaction energy of the 
neutron and proton. This property of the deuteron is extremely important in 
analyses of scattering of neutrons by protons. We now consider this phenome- 
non, recalling certain features of the quantum theory of scattering which will be 
required in what follows. 

An analysis of the scattering of neutrons by protons is somewhat complicated 
by the fact that these particles have spin. However, an analysis of scattering of 
spinless particles reveals the main features of the more complete theory. Thus, 
we shall use the results of the more simple theory of scattering of spinless particles 
in a spherically symmetric force field (i. e., a field in which the potential is 
independent of angle). In the quantum-mechanical formalism scattering in 
such a force field is described by a superposition of plane waves, corresponding 
to the incident particle, and spherically scattered waves, representing particles 
scattered by the force center. 
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A plane wave, propagating along the positive z-axis, is represented by a wave 
function of the form e***? [with an amplitude of unity, this function represents a 
flux of v particles per second through 1 cm?; v = (1/m)hk is the particle veloc- 
ity]. At large distances the spherical wave due to the presence of the scat- 
tering center is given by a function of the form f(0)(e**" /r), where r is the distance 
from the scattering center and f(0) is the scattering amplitude; f(#) is a function 
of polar angle, and describes the angular distribution of scattered particles. 
The positive exponential means that the direction of scattering is from the 
scattering center toward infinity. 

The scattering amplitude introduced in this way is a quantum-mechanical 
feature of the process. It is easy to show that f(#) has the dimensions of length. 
This follows from the fact that f(6)/r must have the same dimensionality as e**?, 
that is to say, it must be dimensionless. Hence, f(0) is sometimes known as the 
scattering length (usually this designation is applied at small particle energies 
since f(9) becomes a constant which is independent of angle and energy). The 
experimentally measured scattering cross sections are related to f(@) by the simple 
relation 


do = | f(e) |*do 


(do is an element of solid angle). 

Quantum mechanics furnishes a method for relating the scattering amplitude 
to the properties of the system. In this analysis it is necessary, first of all, to 
resolve the incident wave into spherical waves: the incident flux with a given 
value of particle momentum, Ak, is given in the form of a sum of fluxes each term 
of which is characterized by a particular value of angular momentum. This 
procedure is valid because in a field of central symmetry the angular momentum 
is an integral of the motion, being described by the azimuthal quantum number L 
(the square of the angular momentum IL 2 = L(L + 1)43). Each spherical 
wave with a given L will have its own angular dependence, which is described by 
the factor P,(cos 9), the Legendre polynomial. The wave with L = 0 is 
spherically symmetric since Po(cos 9) = 1; the wave with L = 1 contains the 
factor Pı(cos 3) = cos Y, the wave with L = 2 contains the factor Pa(cos 3) = 
(3 cos? 4 — 1)/2 and so on. 

It is shown in quantum mechanics that in the resolution of a plane wave into 
spherical waves, at large distances from the origin, the term associated with 
angular momentum L+ is of the form 


L 
won | wR + | Ps (cos 0) (r > ©), 
r 


If there is a scattering center in the path of the incident beam the form of the 
function is changed. However, at large distances from the origin the depend- 
ence on distance remains virtually unchanged because at these distances the 
field has almost vanished and the wave function must satisfy the wave equation 
for a free particle. The scattering process manifests itself in the fact that the 
phase is changed by an amount 6,: 
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1 L 
Y’ Y -— sin [kr + = + | 
Y 


Thus the change of phase (or more simply, the phase shift), ö,, is a parameter 
which describes the scattering of the Lth spherical wave. 
The phase shift, ö,, is connected with the scattering amplitude by a simple 
relation. Specifically: 
218, _ 1 
(0) = iA (QL +1) — Picos ð), 
L 


where 


a | = 


or, in somewhat different form, 


f) =1 Y (QL + 1)ez sin 6,P1(cos 9). 
L 


This expression is complex; to obtain the scattering cross section it is necessary 
to compute the square of the modulus f(P). 
It is apparent that the total scattering cross section is given by the expression 


or = S | f(8) |.d0 = 4rx? Y, (2L + 1) sin? 87. 
L 


For example, for waves characterized by L = 0, we obtain: 


o(L = 0) = 4rX? sin? ôo. 


Scattering of particles with zero angular momentum is called S-scattering. 
It is obvious that the cross section for S-scattering cannot be arbitrarily large. 
Since sin 69 cannot be greater than unity, 


oL(L = 0)<4rA?. 


For waves characterized by L = 1 (P-scattering), L = 2 (D-scattering), etc. 
the limiting cross section is 4r(2L + 1)X?. 

We see that the scattering is completely characterized by the phase shifts. 
Hence, the first problem in an analysis of the results of an experiment is the 
determination of the phase shift; the phase value determined in this way yields 
information on the essential features of the interaction between the particles. 
Strictly speaking, a knowledge of the scattering phase allows one to calculate the 
dependence of the interaction force on distance; this calculation, however, 
requires extremely detailed information on the phase shifts, and has not been 
carried out at the present time. It should be emphasized that the determination 
of phase shifts from measured cross sections is in itself a complicated problem; 
this problem is particularly difficult if the particles have spin. In this case it is 
necessary to carry out experiments with polarized particles to obtain complete 
information. 
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We now consider the possible values which may be assumed by the phase. 

It is obvious that the phases 6, and 6, + ~r yield the same wave function. 
Thus, the phase is determined only to within an arbitrary factor, a multiple of r. 
To eliminate the ambiguity, it is convenient to limit the range of variation of the 
phase. Thus, it is usually assumed that either 


O<6< 7, 


or 


u S0 A 


NIA 


z 
2 
Usually the range of variation is chosen in accordance with the second relation- 
ship. 

Although the dependence of phase on energy is determined by the properties 
of the interaction in an extremely complicated way, it is possible to obtain 
certain general results at low energies. In particular, in systems in which there 
are no Coulomb forces (for example, a system composed of a neutron and a 
proton) and in which the forces act only over small distances, the phase shift 
depends on the energy of the system E (or on the wave vector k) in accordance 
with the formula 


õp ~ kH Ww ~ EL+UVr, 





A2L +1 


According to this formula the phase shifts fall off with energy and the rate of 
reduction is greater at higher values of Z. The S-phase falls off slowest of all 
(i. e., the phase 59) and at very low energies, is the only one which remains; 
thus this quantity characterizes the low-energy scattering completely. Physi- 
cally, this result means that at low energies the scattering becomes almost 
spherically symmetric and the cross section becomes independent of energy. 
The latter result follows from the fact that the cross section is proportional to 
the product of A? and sin? ôo, or A? and 4p”, since ôo is small. However, ôo itself is 
proportional to 1/4, so that the product 4769? is independent of energy. 

The general considerations which have been presented above give no in- 
formation as to the sign of the phase shift at low energies. This sign depends on 
whether the interaction force is one of attraction or repulsion. If the particles 
attract each other, ô> 0, and if they repel, ô< 0. However, these simple con- 
siderations relating the sign of the phase to the type of interaction are valid only 
when the system has no bound states. In a potential well in which there are no 
bound states, 1. e., a comparatively ““shallow” well, the phase shift is positive. 
If the depth of the well is increased the phase shift increases, tending toward 7/2 
as the depth of the well approaches the value at which a level appears. Since we 
are limiting the range of phase variation to — (1/2) < 6< +(r/2), the phase 2/2 
is identical with the phase — (7/2); with further increases in the depth of the 
well the phase shift increases from — (2/2) to 0, remaining negative. 

In determining phases from experimental data it should be emphasized that 
it is more difficult to determine the resultant sign of the phase than its magni- 
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tude. If the signs of all the phases in the scattering cross section are changed 
the expression itself is not changed; in this transformation the scattering ampli- 
tude becomes its complex conjugate and the square of the modulus remains the 
same. Thus, the sign of the phase shift cannot be found directly from scattering. 
To determine signs the scattering phase shifts must be related to a potential for 
which phase shifts can be found theoretically. Such a potential is the Coulomb 
potential. If we analyze a system in which there are Coulomb forces as well as 
the forces which we wish to investigate (nuclear forces in the present case), 
because of interference effects the scattering cross section can be related to the 
sign of the phase shifts of the nuclear forces since the phase shifts associated with 
the Coulomb scattering can be determined theoretically both as to magnitude 
and sign. In this way one determines the phase shifts in proton-proton scatter- 
ing. 

However, in the case of the scattering of neutrons by protons the sign of the 
phase shift cannot be determined in this way; the phase shift can be found only 
by using the fact that there is a bound state in a system comprising a proton and a 
neutron—the deuteron. We shall return to this problem later. 

We now consider the possible forms of the interaction between nucleons when 
spin is taken into account. At the outset we shall limit ourselves to the case of 
nucleons which move with velocities small compared with the velocity of light 
since it is these velocities with which we are concerned in the majority of phe- 
nomena in nuclear physics. For such velocities it is reasonable to assume that 
the interaction potential is independent of velocity. Under this assumption we 
can limit the possible forms of the potential. In the absence of spin one expects 
that the interaction will depend only on the distance between the particles. 
However, in nucleons the interaction can depend on and, as has been shown 
experimentally, does indeed depend on the relative orientation of the spins. 

We now determine the possible forms which the potential can assume. We 
know that the interaction must be determined by the distance between the 
particles r and the spin vectors sı and sz. It is apparent that the potential 
expression will contain one term which depends only on the coordinate Vy(r) and 
other terms which depend on the spins. Owing to the special quantum-mechan- 
ical properties of the spin vector s, for spin !/2 the components can enter only to 
the first power; any quadratic combination of the components of this vector can 


be reduced to one component or to a constant through the use of the following 
relations: 


a ee E : , . 
So? = Sy2 = sè = Vy, SeSy = 1/2 is, SySz = 1/2 152,552 = 1/2 1Sy. 


Hence, only one scalar quantity si-s2 can be formed from two spin vectors and 
the potential energy will contain a term Va(s: - s2), where V2 is some new function 
of the coordinates. The remaining terms can contain the projection of the spin 
along the direction of the radius vector s,-n or s:n where n is the unit vector 
(r/r). However, these products cannot enter individually because the spin 
vector sı is a pseudovector (or, as it is sometimes called, an axial vector). This 
means that in making a transformation from a right-hand system of coordinates 
to a left-hand system this vector reverses direction (as does any angular momen- 
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tum vector). Consequently, the product sı’ n is not a scalar; it is a pseudo- 
scalar and its sign depends on the choice of coordinate system. 

Thus the potential energy, which is a scalar by definition, can contain only the 
product (sı: n)(s:-n) since this is a real scalar quantity. It can be shown that 
it is impossible to form any other scalar expression which is independent of 
velocity; thus, we have the following expression for the interaction energy : 


V(r) = Vir) + VArds: 82) + Var Cs n(s: a). 


This expression includes all possible types of interactions for two protons. 
However, it does not include all possible types of interactions which may exist 
between protons and neutrons. 

It has been shown experimentally that one other process, which is unknown in 
classical mechanics, takes place in scattering. The particles can exchange 
charge during the interaction. The existence of exchange forces is related to the 
high degree of similarity between the proton and neutron—it is believed that 
when these particles are in close proximity a light charged particle is transferred 
from the proton to the neutron (or from the neutron to the proton), thereby 
changing the charge states. The existence of these charged particles is now 
unquestioned: these are the so-called r-mesons, about which more will be said 
later. 

In describing the exchange interaction we proceed by analogy with the pre- 
ceding potential, introducing three additional scalar functions: V4(r), Vs(r) and 
Ve(r) : 


Vexen(r) = [Va(r) + Vs(r)(Si*Se) + Vo(r)(s::n)(Se:n)}P. 


Here P is an operator which exchanges the positions of the ncutron and proton. 
Thus, a description of the interaction in a proton-neutron system requires a 
knowledge of six functions; this is the basic problem of nuclear physics and it is 
far from solved. If the particles involved in the scattering process are spinless 
and do not exchange charges, the problem reduces to the determination of a 
single function; spin and charge, however, make the solution much more com- 
plicated. 

It has been found experimentally that all six terms are of the same order of 
magnitude; thus we have no justification for neglecting any of these interactions. 

The expressions for the potential can be written in a somewhat different form 
in which the spin of the system appears explicitly. The two particles, the proton 
and neutron, may exist in a state in which the total spin is O or in a state in 
which the total spin is 1. The first state is called the singlet state (multiplicity 
equal to 1), the second state is the triplet state (multiplicity equal to 3). We 
may recall that in general the multiplicity of a system is given by the number 
25 + 1, where S is the spin (the multiplicity is equal to the number of possible 
orientations of the spin vector S). In states with S = 0, it is obvious that the 
spin cannot play a role, and in these states the potential for each of the interac- 
tions (simple and exchange) can depend only on the distance; we designate this 
case by UN + Urexcn(7)P. In the state with S = 1 the spins are parallel and 
the potential can be divided into two parts: one is independent of the direction 


10 Lectures on Nuclear Theory 


of the total spin in space and is written in the form U(r) + Ue exen(r)P, while the 
second depends on this direction. The dependence on direction in a system of 
two particles implies a preferred direction—that of the vector which connects the 
two particles; hence, the directional dependence of the potential is given by the 
pseudoscalar S-n. 

Because of the pseudoscalar nature of the product S-n this quantity can 
enter the potential only in quadratic form (S «n)?. Higher powers are excluded 
by the quantum-mechanical properties of the vector S for S = 1. 

It is customary to write the potential in a more convenient form. In par- 
ticular, the term which depends on the spin direction is written so that its 
average over direction becomes zero. Averaging the expression (S:n)? = 
D,5,S,nın, over all directions of the vector n and using the relations 


nin = */305% 


(these follow directly if it is recalled that the n; are the direction cosines of the 
lines which connect the particles), we obtain: 


2 SiSinin = 1/3 2, SiSrdır = 1/3 2 S: = 1!/S(S + 1) =?/3. 
t, 1, ? 
In the last expression we have used the properties of the square of the quantum- 
mechanical angular momentum vector for S = 1. 

Using the expression which has been derived we write the factor which deter- 
mines the directional dependence of the potential in the form 


[Udr) + Usexcn(r)P] (Sn)? — ?/s). 
Thus, in the singlet state the potential 1s 
Using = Ur) + Ui exch(7)P, 
while in the triplet state 
Utrip = [U2(r) + Usexen(r)P] + [Us(r) + Usexcn(r)P][(S°n)? — 2/5]. 


It is clear that the number of functions required to describe the interaction 
(namely, six) remains the same, regardless of the representation. There are no 
general relations which can be used to reduce this number. It is impossible to 
establish any general relation between the interactions in the singlet state and 
the triplet state—the connection between the functions U,(r) and U2(r) arises as 
a result of the properties of the actual nuclear system, consisting of the proton 
and neutron in the present case. Experiments indicate that the last term—the 
potential due to tensor forces *—is as important as the other terms. 

One of the first experimental findings which indicated the role of the tensor 
forces was the existence of the quadrupole moment in the deuteron. The spin 
of the deuteron is unity—this is to be expected as the result of adding the two 
parallel spins of the proton and neutron, each equal to 1/2. If there were no 


* The designation ‘‘tensor’’ is due to the fact that the spin enters the potential in 
the form of a tensor with components S;S;. 
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tensor forces, the spin of the deuteron would be independent of its axis—the line 
which connects the proton aud neutron. This would imply that the charge 
distribution in the deuteron is independent of the spin and the spin would not 
define a preferred direction; it follows that the charge distribution would be, on 
the average, spherically symmetric. It is well known that a spherically sym- 
metric charge distribution gives rise to a field which is identical with the field of a 
point charge, 1. e., such a system has no electric moments; in particular, it has 
no quadrupole moment. For this reason the existence of a quadrupole moment 
in the deuteron is to be interpreted as an indication of an asymmetry in the 
charge distribution. 

In terms of its electrical properties the deuteron can be considered a charged 
ellipsoid of revolution with axis of rotation parallel to the spin; the ellipsoid is 
elongated rather than flattened since the quadrupole moment of the deuteron is 
positive. 

The deuteron quadrupole moment is actually quite small. However, this 
situation does not arise because the interaction forces are weak but rather be- 
cause the deuteron dimensions are large. Because of the small binding energy 
the proton and neutron are, on the average, relatively far apart (long range of 
nuclear forces). It can be shown by calculation that the tensor forces must be 
quite large to produce even a small quadrupole moment in a system with large 
distances between particles. 

Thus, it turns out to be difficult to describe even the simplest nuclear system— 
a system of two particles. For this reason the study of nuclear forces is being 
carried out at the present time with considerably less information than is avail- 
able in the case of atomic systems. In the latter the Coulomb field gives rise to 
comparatively simple features which are well known. 


LECTURE TWO 


Nuclear Forces 


(SCATTERING OF NUCLEONS BY NUCLEONS) 


A, we have already indicated, experimental data on the proton- 
neutron system are acquired by studying the scattering of neutrons by protons 
and by investigations of the properties of the deuteron. An analysis of these 
data, the goal of which is the determination of phase shifts, requires an expansion 
of the scattered waves in terms of spherical functions. Because the proton and 
neutron have spins this procedure is more complicated than in the case of spinless 
particles, which have been considered in the last lecture. First we must introduce 
a system for classifying the possible states of the system. In the case of spinless 
particles the state is determined by the azimuthal quantum number £. If the 
particles have spin the state of the system is also a function of the spin quantum 
number S. If the spins and coordinates are independent, the orbital angular 
momentum and the spin angular momentum are conserved individually. Any 
coupling between the spin and orbital motion means that the vectors L and S 
are not conserved individually. Instead the total angular momentum J is 
conserved;* this quantity is the sum of the orbital and spin angular momenta. 

In nuclear systems, as in all other quantum-mechanical systems, in addition 
to the momentum J there is one other quantity which is conserved: this is the 
parity, a quantity which has no classical analog. The parity determines the 
properties of the wave function of the system with respect to reflection through 
the origin of coordinates. The parity is positive if the wave function of the 
system does not change sign under such a reflection; in this case the system is 
said to be even. The parity is negative if the wave function changes sign under 
such a reflection (an odd system). The division of systems in terms of parity is 
similar to the division of vector quantities into real and pseudo quantities which 
was discussed earlier. 

The states of a system composed of two particles are designated as in atomic 
spectroscopy. Capital letters are used to designate the orbital moment: the 
azimuthal quantum numbers L = 0, 1, 2, 3, 4, 5, etc., are denoted by the letters 
S,P,D,F,G,.... The multiplicity index 2S + 1 is given above and to the left 
while the value of J is given below and to the right. Thus, ¿Dz designates a 
triplet state with Z = 2 and a total moment of 3. 


* In what follows the quantum number designating the total angular momentum 


will be called J. 
13 
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First we shall write out the singlet states. Since the spin in these states is 
zero, J coincides with L and the classification is the same as for particles without 
spin: 150, 1P,, De, IPs, 1G4, etc. 

In the triplet states the classification scheme is somewhat more complicated. 
When L = 0, we have only one state since the spin, equal to 1, added to L = 0 
yields unity. This is the êS; state. In the state with L = 1, we have two vec- 
tors—the orbital moment and the spin, each of which is unity. Adding these, 
we obtain three possible values for the total moment, 0, 1, and 2. Thus, there 
are three states (corresponding to the multiplicity of three): *Po, ?Pı, and ?P;. 
Likewise, when L = 2 we obtain the triplet 3D,, ?D,, and ?D3. Continuing this 
process we arrive at the triplet system: 

¿Po `D; 3 Fo 3G3 
351, {3Pır, {3Da}, (3F3}, 3*G4) etc. 
¿P, 3D; 3 F; 3Gs 


The 3S, state is called a triplet state to maintain the convention although there 
is only one state when L = 0. 

Having established the possible values for the angular momentum (the quan- 
tum number L) and the spin of the neutron-proton system, we can now carry 
out the classification of the states of this system. A good quantum number 
for the system is J; this quantity characterizes the total momentum and the 
parity of the system. In a system of two particles the parity is determined 
uniquely by the azimuthal quantum number L. Since the wave function of the 
system with azimuthal quantum number L contains the Lth Legendre poly- 
nomial P; (cos 9) and since the Legendre polynomial is multiplied by the factor 
(—1)” under transformation from a right-handed coordinate system to a left- 
handed system,* the parity of the system is determined by L and is (—1)?. 
Thus we see the origin of the name for this quantum characteristic. 

Thus, the S, D, G states are even states while the P, F states are odd states. 
If we now consider the possible values of J and L we find that in the triplet 
states the value of J and the parity, in the general case, do not uniquely deter- 
mine the value of L. Lis determined uniquely only in triplet states such as the 
¿D, state, in which J coincides with L, and in all singlet states. In triplet 
states in which J = L + 1 there are always two states with the same value of J 
and the same parity but with different L (aside from the value J = 0 (odd) 
which gives rise to only one state, the ?P, state). These two states “combine” 
with each other. This means that the actual state of the system is a super- 
position of two states with different values of L (obviously these must differ by 
2 units). We may say that the system spends only part of the time in each 
state with different L. Thus, for example, the deuteron is a system with J = 1 
(even) and spends approximately 96% of the time in the 3S, state and 4% of 
the time in the êD; state. 


* In this transformation Y goes over to r — % and cos Y to — cos 8. In other 
words, a Legendre polynomial with even L is an even function of its argument while a 
polynomial with odd L is an odd function. This is the basis for the statement which 
follows. 
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It is easy to enumerate the possible states for an n- system. These states are 
designated by a number (the value of J) and by a superscript + or —, denoting 
the parity of the system. Thus, in this classification system the deuteron state 
is designated by (1+). If the state is actually a superposition of two states with 
different L we designate it symbolically by the sum of these states: for example, 
35, + 3D; for the deuteron. 

The possible states of the n-p system are shown in Table 1. 


TABLE 1 
Singlet Triplet 
(O*)1S9 (07)Po 
(1-)1Pı (1*)3S1 + °Dı 
(2*)1D, (17)? P: 
(3—)!F; (2+)39D, 
etc. | (2)%P, + 3 Fy 
| (34D, + 3G; 
| etc. 


A division has been made between singlet and triplet states. However, it 
does not necessarily follow from the above that singlet and triplet states cannot 
combine (for example, 1D. and 3D: states). However, it can be shown that if the 
energy of the interaction depends only on the coordinates and spin (i. e., if the 
interaction is velocity independent) these states cannot in fact combine. 

One of the important features of the n-p system is the existence of a bound 
state, the deuteron. The deuteron has a binding energy of 2.23 Mev. As has 
already been mentioned, this quantity does not reflect the strength of the inter- 
action between the neutron and the proton. The actual “potential well’’ is 
considerably deeper than 2.23 Mev. If we describe the interaction in terms of 
the potential-well model the depth of the well must be greater than 30 Mev. 
The relatively small binding energy indicates merely that the depth of the “well” 
is not much greater than that required for the existence of the first level. If 
the binding energy were strictly zero, the following relation between the depth 
of the “well” and the radius would be satisfied : 


2 À 
T 
Vro? ~ z — = 10-24 Mev:cm?, 
p 


For the radius ro = 2 107! cm (which must be taken to agree with experiment), 
V = 25 Mev; actually, V is one and a half times larger. These numbers, 
among other things, illustrate the extent to which the well must be “deepened” 
in order to increase the binding energy from 0 to 2.23 Mev. 

How should the deuteron state be classified? Experiment indicates that the 
spin of the deuteron, i. e., the quantum number J, is 1. The table given above 
contains three states with this value of J. These are the two states (1—)!P; and 
(1-)?P, and the state (1*)3S; + *D,. How is the parity of the deuteron to be 
determined? The most reliable information is obtained from experiments on 
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the scattering of neutrons by protons. This effect has now been investigated 
very carefully over a wide energy region, ranging from hundredths of an electron 
volt to hundreds of millions of electron volts. 

It has been shown experimentally that this scattering is spherically symmetric 
(in the center-of-mass system) for neutron energies ranging from zero to approxi- 
mately 20 Mev. This region is undoubtedly to be associated with S-scattering. 
Using the magnitude of the scattering phase shift at low energies and an analysis 
of the experimental data we can determine whether or not the system has a 
shallow S-level and can calculate the energy of this level if it exists. Such 
analyses of the scattering of neutrons by protons have shown that the system 
has a level corresponding to the 3S, state and that this level has an energy equal 
to the binding energy of the deuteron. Actually, there is a small admixture of 
the 3D, state which does not affect scattering at low energies. 

Scattering cxperiments also determine the phase shift associated with the 'So 
state; this phase shift turns out to be large at very small energies because the 
scattering at these energies is of resonance nature. This means that the inter- 
action is very similar to that in the system in which there is a level at E = 0. 
The answer to the question of whether there actually is a level in the system 
depends on the sign of the phase shift; if a level does exist this phase shift should 
be negative. However, it is impossible to determine the sign of the phase shift 
from the scattering of neutrons on free protons—the cross section does not 
change if the sign of the phase shift is reversed. To determine the signs it is 
necessary to investigate effects in which interference takes place between scat- 
tering in the 1S) and 4S; states. Such an effect is the scattering of slow neutrons 
on hydrogen molecules in which the spins of both nuclei are parallel (ortho- 
hydrogen) or anti-parallel(para-hydrogen). It has been established from these 
experiments that the sign of the phase shift is different in the 35, and 15, states. 
Since there is a level corresponding to the 3.5, state (deuteron), it follows that the 
phase shift is positive in the 15, state and there is no level. 

Although there is no level associated with the tSo state, in analyzing phase 
shifts it is convenient to introduce the idea of a “virtual” level. 

Lct Vo be the depth of the well for which a level exists (i. e., the depth for 
which the level energy is zero). Then, a well of depth Vo + ôV, where ôV is 
small, will have a level, while a well of depth Vo — ôV will not. The scattering 
phase shift at zero level-energy for a potential Və is by definition 7/2. If the 
well is now made somewhat decper (potential Vo + ôV) the phase becomes some- 
what larger. According to the convention which has been adopted the phase 
shift must lie within the interval [-(#/2), 7/2]. Hence, an increase in phase 
means that it becomes negative, assuming a value close to —(7/2). On the 
other hand, the potential Vo — ôV obviously leads to a positive phase with a 
value close to 7/2. 

A more exact analysis shows that in both of these wells the scattering at low 
energies is of the same magnitude, but of opposite phase. We can imagine that 
the real level in the Vo + ôV well corresponds to a “virtual” level of the same 
magnitude in the V — ôV well. This terminology is especially convenient 
since the scattering cross section at low energies is independent of the sign of the 
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phase shift; that is to say, it does not depend on whether the level is real or 
virtual (scattering from a Vo + ôV well and a Vo — ôV well is the same). 

It turns out that the depth of the “virtual” level for the 15, state is approxi- 
mately 70 kev. 

Thus, at low energies a scattering analysis makes it possible to find a param- 
eter which characterizes the interaction from the measured cross section. It 
also turns out that an investigation of the energy dependence of the scattering 
allows us to find one more parameter. Hence, two parameters can be deter- 
mined from experiment. These parameters acquire more significance when used 
to describe the interaction in terms of some concrete model. For example, if 
one chooses a rectangular well, the depth and width can be found experimentally. 

Scattering theory indicates that it is impossible to determine any more than 
two constants from the scattering of low-energy neutrons. Hence, the true 
shape of the potential cannot be found in experiments with low-energy neutrons 
even in the simple case of the singlet state in which the interaction, as has already 
been pointed out, is described by a single function.* The experimental data 
can be fitted to almost any potential as long as the interaction vanishes at large 
distances. Different forms of the two-parameter potential are used in different 
work. Typical examples are a rectangular well (with parameters radius r and 
depth V, as described above) and a potential of the form 


ear 


Vo 





r 


with the parameters being Vo anda. The actual shape of the potential has still 
not been determined because even experiments with high-energy neutrons pro- 
vide no data for solving this problem. Consideration of the concrete form of the 
potential can, at present, serve only to give a picture of the interaction. For 
instance, it can be shown that scattering in the 15, state can be described by a 
rectangular well with a radius of 2.8 : 1071% cm and a depth somewhat greater 
than 10 Mev. 

The situation is much more complicated in the triplet state 35,. In the triplet 
states, as has been pointed out above, the tensor forces play an important role 
and the 3S, state is mixed with the ?D, state. If, for the sake of argument, we 
neglect the tensor forces, this state can be described by scattering in a rectangular 
well. In this case the triplet-state ‘“‘well’’ is much narrower than the singlet well 
and the depth is approximately three times greater. These numbers give some 
idea of the effect due to the relative orientation of the spins of the colliding 
particles. 

Slow-neutron experiments give no information which can be used to estimate 
the magnitude of the tensor forces, Only one quantity is available for making 
such an estimate—the quadrupole moment of the deuteron, and this does not 
provide sufficient information for a quantitative description. It can only be 
stated that the tensor forces are of approximately the same order of magnitude 


as the ordinary forces. 


* In the 1S, state this function is the sum of the potentials V; + V1 exch. 
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We now turn to a second system: the p-p system. It differs from the n-p 
system primarily in that there is also an electric interaction between protons. 
The presence of the nuclear interaction changes the distribution for proton- 
proton scattering from that which is calculated on the basis of a Coulomb law 
(Rutherford scattering). 

Identity considerations must be taken into account when two protons interact. 
This is in accordance with the requirement imposed on such a system by the Pauli 
principle. According to the Pauli principle the wave function for a system of 
two protons must be anti-symmetric with respect to interchange of coordinates 
and spins for both particles (it must change sign in such an interchange). In 
triplet states of the system the wave function does not change sign upon inter- 
change of particle spins; in the singlet states the wave function does change 
sign in such an interchange. To make the total wave function anti-symmetric 
it must change sign upon interchange of particle coordinates in triplet states (the 
states must be odd states) and must not change sign for interchange of coordi- 
nates in singlet states (the singlet states must be even states). It follows from 
these considerations that the p-p system has available only half the states which 
are allowed for the n-p system (in which there is no limitation on the parity of the 
state and the symmetry of the wave function). 

The allowed states for a p-p system are shown in Table 2. 


TABLE 2 
Singlet Triplet 
(O*)150 (0~)*Po 
(21D, (17)3Pı 
(4*)1G, (2-)8P, + 3F, 
etc. (37) Fs 


etc. 





Proton-scattering experiments at energies of several millions of electron volts 
yield important results concerning the interaction in the IS), state. A scattering 
analysis allows us to reduce the experimental data to the determination of two 
constants. An important experimental fact, discovered in proton-proton 
scattering experiments, is that the proton-proton and proton-neutron inter- 
actions are virtually identical. 

If we write the interaction potential for two protons in the 15, state in the form 


e2 
te MR), 
Y 


i. e., as a sum of the Coulomb potential and the nuclear potential, the parameters 
for the potential V(r) virtually coincide with the parameters for the potential 
Vi + Uy, exchP in the +5, state of the n-p system. Hence, it follows that there is 
no stable level in the p-p system (just as in the 1So state of the n-p system). In 
other words, the nucleus He? cannot cxist in nature. 


Nuclear Forces 19 


The similarity of the interaction in p-p and n-p systems is of fundamental 
significance. This similarity indicates the deep-lying symmetry between these 
particles. This symmetry goes far beyond mere charge symmetry; it appears 
in the identical spectra of mirror nuclei and gives rise to the identical forces be- 
tween protons and between neutrons. All the presently available data on 
nucleon behavior at all energies indicate that p-p, p-n, and n-n systems are very 
similar, not only in scattering in the 15, state but in all effects in which the 
Coulomb force is relatively unimportant. This property is called isotopic (or 
charge) invariance. 

It should be emphasized that isotopic invariance by no means implies com- 
plete equivalence in the interactions for the three different particle pairs; for one 
thing the n-p system has states which cannot exist in identical-particle systems. 
It is clear that the interactions in these states are not directly related to the 
properties of the p-p or n-n systems. In speaking of isotopic invariance we are 
discussing the interaction in identical states of different systems. 

In considering isotopic invariance it is convenient to introduce a formalism 
which will be found very useful in obtaining information as to the nature of 
different nuclear systems. We present the basic features of this formalism, 
which is purely an abstract description. 

The starting point of this scheme lies in considering the neutron and proton 
as two different “charge states” of the same particle. These two charge states 
of the nucleon are described by means of a new concept—+sotopic spin.* In this 
scheme the particles are characterized by a new vector—the isotopic spin vector 
7; its absolute value is 1/2. As for ordinary spin */, its projection along a 
preferred axis in isotopic spin space—the { axis—can have only two values (the 
other two axes are designated by ¢ and y). The value +1/, indicates that the 
particle is a proton, and the value —!/2 indicates that the particle is a neutron. 
The transition from neutron to proton is denoted by a change in the value of the 
projection from —1*/, to +*/2. By analogy with ordinary space this operation 
can be described as a rotation through 180° about the & axis in isotopic spin 
space. (In ordinary space this would be a rotation about the x or y axes.) If 
the system consists of several nucleons the isotopic spins are added in accordance 
with the rules for quantum-mechanical addition of vectors. In such a system 
the projection of the total isotopic spin on the ¢ axis is equal to the neutron 
excess! N — Z multiplied by —!/,, since each proton contributes +!/2 to the 
projected value and each neutron contributes — 1/2. 

The concept of an isotopic spin space is purely formalistic and abstract. 
Hence, it is meaningless to discuss the relation between this space and ordinary 
space; similarly, it is quite pointless to try to extend the interpretation given 
above. 

The convenience of isotopic spin space lies in the fact that the isotopic invari- 
ance features can be formulated compactly. If we neglect Coulomb forces, a 
system consisting of two protons differs from a system consisting of one proton 


* Isotopic spin = isobaric spin in more recent American usage.—Translator. 
+ Z is the atomic number and N is the number of neutrons. 
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and one neutron only by a rotation in isotopic space. Analogous states of these 
two systems (that is, states with the same quantum numbers) can be considered 
as states corresponding to different projections of the isotopic spin of the two- 
nucleon system along the ¢ axis in isotopic space. In these terms isotopic in- 
variance is described in terms of the invariance of a system with respect to 
rotation in isotopic space. We may note that charge symmetry—invariance 
of the properties of the system under replacement of all neutrons in the system 
by protons, and vice versa—appears as a particular case of isotopic invariance; 
in this formalism we say the system is invariant under a rotation of 180° about 
an axis in the En plane in isotopic space. 

Obviously, isotopic invariance does not extend to the Coulomb interaction. 
This can also be shown in a formal way: the charge is given by the projection 
of the isotopic spin; consequently the Coulomb interaction depends on the 
projection of the isotopic spin, a quantity which is obviously not invariant 
against rotation of the system of coordinates. 

The introduction of isotopic spin makes it possible to extend the classification 
system for the states of two-nucleon systems through the introduction of a new 
quantum number, the isotopic spin. Since each of the nucleons has an isotopic 
spin equal to one-half, according to the rules for adding quantum-mechanical 
vectors, the vector sum of these must be either unity or zero. Hence, all states 
of the system can be divided into two groups: one is associated with isotopic 
spin T = 1, the other with T = 0. 

An isotopic spin T = 1 can have three values for its projection. In accordance 
with rules given above, the projection Ty = —1 is to be associated with a 
charge O for the system, the projection T; = O with a charge +1 (in units of e, 
the charge of the electron) and the projection 7; = 1 with charge +2. It is 
apparent from these considerations that a system with T = 1 can be realized by 
two neutrons, two protons, or a neutron and a proton; in the latter case we can 
associate with states T = 1 only those states of the p-n system which can also 
exist in a system of two identical particles, 1. e., states which satisfy the Pauli 
principle. 

The state with T = O can have only one projection, T} = 0; this describes all 
states of the n-p system which do not satisfy the Pauli principle. It is apparent 
that these states are described by a symmetric wave function (these states are 
even triplets and odd singlets). 


Finally, we arrive at the following classification of states: 


for T = 1 150, Po, 1, 2, etc.; 
for T = 0 3So, Pa, etc. 


It is interesting to note that in this scheme the triplets and singlets (for a given L) 
always refer to different values of T, and hence cannot combine with each other. 

Using the idea of isotopic spin, it 1s now possible to formulate exactly the 
properties of isotopic invariance of nuclear forces as applied to a system of two 
nucleons: the isotopic invariance now expresses the fact that the interaction 
is the same for any two nucleons in states with T = 1. 
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Before turning to a description of the experimental results with high-energy 
nucleons, something should be said about nucleon interactions at low energies in 
states with L > 0. 

We have already indicated that it is only at low energies that the interaction 
is basically described in terms of the S-state. Theoretically, it is impossible to 
predict at what energy other states start to appear. It has been shown experi- 
mentally that even at 20 Mev the role of these other states is still small and a 
determination of the residual phase shift remains difficult. At 40 Mev other 
states play an important role in the p-n system. It is unfortunate that it is 
precisely this energy region (20 to 40 Mev) which is almost completely unin- 
vestigated. At an energy close to 20 Mev, the phase shift in the P-state is 
estimated as approximately 1° in the p-p system; because of the small value of 
this phase shift it is impossible to determine which of the three P-states is re- 
sponsible. 

Nothing more can be said with regard to this state. Experiments, which we 
will discuss later on, indicate that the interactions in all three P-states (Po, ?P,, 
3P.) are considerably different, and hence any statement as to the average inter- 
action in the three P-states would have very little validity. Still less can be said 
concerning the parameters for this interaction. The interaction is characterized 
here by the product of the well depth and the fifth power of the radius (in the S- 
state it would be the square of the radius); hence, the well depth is very sensitive 
to the choice of radius. 


LECTURE THREE 
Nuclear Forces 


(SCATTERING OF NUCLEONS AT HIGH ENERGIES) 


W. have seen that experiments on nucleon-nucleon scattering 
at low energies yield information only on the interaction in ?S, and 15, states 
and are not very useful in determining the interaction potentials between nucle- 
ons. The reason lies in the fact that at energies of the order of 10 Mev the 
nucleon wavelength is of the order of the interaction radius and the scattering is 
insensitive to the detailed dependence of the interaction on distance. In order 
to “probe” the potential well the wavelength must be considerably smaller than 
the dimensions of the well. Hence, more information is obtained from the 
experiments with high-energy nucleons (several hundreds of millions of electron 
volts) which have been carried out in the last several years in various laboratories. 

These experiments are all the more interesting because quantum mechanics 
shows that it is possible to obtain direct information on the potential which acts 
between the particles from the scattering data at these energies without recourse 
to a phase-shift analysis. This possibility is based on the assumption that the 
interaction energy is small compared with the particle energy; this assumption 
is certainly valid at high energies. However, we may not be justified in assum- 
ing that the interaction is “weak”; the picture at high energies may not be the 
same as that which would be expected on the basis of ordinary perturbation 
theory. 

Nevertheless, at the outset we base our analysis of scattering at high energies 
on perturbation theory since this approach turns out to be useful in the later 
analysis. 

When the kinetic energy of the colliding particles is much greater than the 
interaction energy, the scattering amplitude in the center-of-mass system (c.m.) 
is given by the following expression, which is derived in quantum mechanics 
(the Born approximation): 





H / 
wo) = st $ reta 
W) = zin J vo 
In this expression u is the effective mass, V(r) is the interaction potential, k and 
k’ are the wave vectors before and after scattering, Y is the angle between K and 
k’ (scattering angle) and dr is a volume element. 

23 
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We may recall that the scattering cross section is given in terms of the scatter- 
ing amplitude by the expression 


do(9) = | f(a) |?do 


and that 9 and q = k — k’ are related by the expression (for elastic scattering) 
0 
q? = (k — k’)? = 42? sin? a 


It is easy to show that the scattering amplitude (to within a multiplicative 
constant) is simply the Fourier amplitude of the interaction potential. Thus, 
if the scattering theory formulas apply, having determined the scattering ampli- 
tude, we can find the potential through the properties of the Fourier integral. 

These formulas describe the angular distribution of the scattered particles. 
For purposes of illustration we choose a potential which differs from zero only 
in a region of dimensions a. 

In the scattering amplitude formula the factor e**”* under the integral sign is a 
periodic function with period (in the direction of the vector q) 27/q. If this 
period is small compared with the dimensions of the well, this oscillating factor 
changes sign many times within the boundaries of the well. Thus the integral 
in the expression for f(#) will be small since V(r) varies at a much slower rate 
than e% and the different parts of the integral will approximately cancel. The 
larger the value of g, the more effective is this cancellation effect; hence when 
q > 2r/a the value of the integral may be taken as zero. This integral will be 
different from zero only for 


iq-r 


~~ 2r 
AA 
a 
On the other hand, at small 4 the magnitude of the vector q is related to the 


scattering angle by the expression 
q = kd. 


Whence it is obvious that the integral will be nonvanishing only for 9 < 27/ka. 

Thus, the angle 3, within which the particles are scattered, is inversely propor- 
tional to the wave vector, i. e., the square root of the particle energy. In this 
case, the amplitude corresponding to scattering in the direction Y = ( is inde- 
pendent of particle energy and is given by the expression 


2u 
10) = — = f V(r)redr. 


It can be shown from perturbation theory that at high energies the scattering 
cross section is constant only within a small solid angle, 32, which (in the c.m. 
system) is inversely proportional to the energy of the particle. Under these 
conditions, the total scattering cross section is obviously inversely proportional 
to E. 


The scattering picture described above is changed if the forces which act 
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between the particles have exchange character. The nature of this change is 
easily understood if one considers, in the problem which has just been discussed, 
the angular distribution of the recoil particles rather than the scattered particles 
(the recoil particles are those which are at rest in the laboratory system before 
the collision). In the center-of-mass system the recoil particles move in a 
direction opposite to that of the scattered particles and thus the angular distribu- 


tion of these particles is described by the same scattering cone but in the back- 
ward direction (the region Y ~ r). 





_ mbarns/sterad 
Y 
sa i 





ge 


0° 10° 0° 50° 60° 70° 80° 90° 
g 
Fig. 3. Elastic proton-proton scattering (from data reported by N. P. Boga- 
chev, I. K. Vzorov, M. G. Meshcheryakov, B. S. Neganov and E. V. 
Piskarev). Proton energy 460 Mev (almost isotropic scattering) and 657 
Mev (highly anisotropic scattering); Y is the proton scattering angle in the 
center-of-mass system. 


Exchange scattering means that charge is exchanged between particles during 
the scattering process; that is, an incoming neutron is scattered as a proton, 
while the recoil particle proves to be a neutron. 

It is clear that the neutron distribution in exchange scattering is the same as 
the distribution of recoil particles in ordinary scattering. In the general case the 
scattering diagram consists of two directional cones: these are in the forward 
and backward directions (in the center-of-mass system), and the relative intensi- 
ties are determined by the fraction of exchange force in the interaction between 
the particles. 

The experimental results published in the last few years are found to be com- 
pletely unlike the picture which has been given. A clearcut contradiction is 


26 Lectures on Nuclear Theory 


found in experiments on the scattering of protons by protons. These experi- 
ments have shown that over an extremely wide energy range (150 to 350 Mev) 
the cross section is independent, within the experimental errors (approximately 
10%), of both angle and energy. It is only in the region of extremely small 
angles that the cross section increases and this effect is due to ordinary Coulomb 
scattering. The cross section per unit solid angle is approximately 


o(9) = 4 mbarns/sterad. * 


It is only at higher energies (430 Mev) that a noticeable angular dependence is 
found. In recent experiments carried out with the Soviet synchrocyclotron the 
scattering cross section has been investigated at 657 Mev (Fig. 3). At these 
energies the cross section shows a pronounced angular dependence. 

These experimental results indicate that the basic assumption of perturbation 
theory—that the interaction is weak at high energies of the incoming particles— 
is not justified. The experimental results indicate that extremely powerful 
forces come into play when the distance between particles is small. To carry 
out a more detailed analysis of these forces, we must make a phase-shift analysis 
of the angular distribution. At the present time, however, the data which have 
been accumulated are still inadequate for a complete analysis of this type. 
Nevertheless, 1t is possible to show why the cross section is independent of scat- 
tering angle over a wide angular region. 

It has been shown above that there is a limitation on the magnitude of the 
scattering cross section. In particular, the cross section for S-scattering (that is, 
for L = 0) per unit solid angle cannot exceed 


domax _ 


do 


2 


(this formula applies both for neutron scattering and proton scattering in S- 
states). 

At 380 Mev, X? is approximately 2.5 mbarns; this value is smaller than the 
observed cross section. It might be said (and this statement is encountered in 
the literature) that there is a clearcut inconsistency. However, such is not the 
case. There is one other state (in addition to the S-state) which is characterized 
by spherical symmetry in scattering—the *Po state. The maximum cross section 
for this state is also 1?; thus, the states So and *Po together can give rise to an 
isotropic cross section of 242, which is consistent with experiment. It is obvious 
that the phase shifts of the other states with L = 1 (èP, and *P» states) must 
either be small or have highly anisotropic cross sections. From these consider- 
ations we conclude that the interaction is essentially different in the *Po, 3P, and 
®P, states; that is, the spin-orbit interaction is large in a system of two protons. 
This behavior in a system of two particles has already been encountered in the 
deuteron. In that case it is responsible for the small quadrupole moment 
associated with mixing of Sy and 3D, states. While the effect in question is 
small in the deuteron, at higher energies it becomes predominant. 


* 1 barn = 10724 cm’. 
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Further information on the proton-proton system can be obtained from 
investigations of scattering of polarized particles on nonpolarized particles and 
scattering from polarized targets. The problem of making a target of polarized 
protons is an extremely difficult one and has still not been solved. It is only in 
the last few years that the first experiments on the scattering of polarized 
particles on unpolarized targets have been reported. 

When a polarized beam is scattered on a target (analyzer) there is a difference 
in the scattering intensity in the forward and backward directions. The differ- 
ence in the forward and backward cross sections also yields information as to the 
polarization of the beam. At the present time the results of a number of experi- 
ments on the scattering of polarized protons by protons have been reported. 
Although these data are as yet not very accurate, the importance at these 
energies of states other than 15, and *P, states is apparent. At about 200 Mev 
the phase shift of the êP» state reaches 7-10°, corresponding to a polarization of 
approximately 20%. At higher energies the polarization increases and the 
interpretation of the data at 400 Mev requires that we consider states with L > 1 
in addition to states with L = 1. These results are still far from reliable; 
further refinement will require further experimental developments. Additional 
information on proton-proton interactions should be given by experiments in 
which polarized targets are used in conjunction with polarized proton beams. 
However, work with unpolarized targets can be useful if the experiments are 
improved so that observations on the polarization of the recoil proton are 
carried out together with observations of the polarization of the scattered pro- 
tons. The results of such experiments should make it possible, for example, to 
distinguish between scattering in 15, and *Po states; such a distinction is not 
possible in ordinary scattering experiments (as we have seen, both of these 
states are characterized by isotropic scattering). 

An important feature of high energy scattering (higher than 400 Mev) is 
inelastic scattering, i. e., Scattering accompanied by the production of r-mesons. 
Because of the production of r-mesons the proton-proton interaction changes 
markedly at energies in the vicinity of 1000 Mev. The total cross section for the 
interaction (elastic plus the production of r-mesons) increases, reaching a value 
of approximately 50 mbarns at 800 Mev and remaining constant [within the 
limits of the errors (10%)] up to 1000 Mev, the highest energy at which these 
measurements have been carried out. The elastic scattering, however, remains 
virtually unchanged at approximately 25 mbarns (the value at the lower ener- 
gies). Thus the elastic and inelastic cross sections are essentially equal and 
constant in this region. It can be shown theoretically (we present these calcu- 
lations in the Eighth Lecture) that this behavior is characteristic of scattering of 
particles on an absolutely black sphere, that is to say, a sphere which absorbs all 
particles incident on its surface. Because of the wave nature of the proton flux 
this absorption is accompanied by diffraction at the edges of the sphere (analo- 
gous to the diffraction of light at the edge of a black obstacle); this effect is also 
observed in elastic scattering. 

The absorption cross section (production of r-mesons) in this model is +R; 
this is also the value of the diffraction cross section for elastic scattering. If this 
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Fig. 4. Elastic neutron-proton scattering (from data reported by V. P. 
Dzhelepov and Yu. M. Kazarinov). Neutron energy (in the laboratory 
system) 380 and 590 Mev; # is the neutron scattering angle in the c.m. 
system. 


model is used to describe the scattering of protons by protons, taking the sum 
of the “radii” of the two protons as the collision parameter, the proton radius is 
found to be 


fp = 4.5: 1071 cm. 


It is interesting to note that these results are in agreement with the results of 
experiments on the scattering of electrons by protons. Such experiments 
indicate that at distances smaller than approximately 7 - 1014 cm the electric 
field of the proton does not correspond to that of a point charge. 

Elastic scattering of particles on a black sphere should also be characterized 
by an angular distribution typical of diffraction. Such a distribution is observed 
experimentally. It will be extremely interesting to examine the behavior of the 
cross section at higher energies. 
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We now turn to the scattering of neutrons by protons. This type of scattering 
has been studied over a wide energy range (up to 590 Mev) and yields a picture 
which is considerably different from that described above. The general energy 
dependence of the cross section may be described as follows: the cross section 
is a strong function of angle, increasing by several times at 180° as compared 
with the cross section at 90° (Fig. 4). Thus, for example, at 380 Mev the 90° 
cross section is 2 mbarns while the 180° value is 12 mbarns. At higher energies 
the 90° cross section falls off, becoming 1 mbarn at 590 Mev; however, the 
cross section at 180° remains almost constant (measurements reported by V. P. 
Dzhelepov and Yu. M. Kazarinov). 

Below 90° the cross section does not increase so rapidly and the curve is asym- 
metric (till quite recently it had been assumed, without justification, that the 
scattering curve was symmetrical about 90°). 


6/¥)} mbarns/sterad 
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Fig. 5. Approximate dependence of the nucleon scattering cross section 
in states with isotopic spin T = 0 and T = 1. The nucleon energy is 
approximately 400 Mev; ð is the scattering angle in the c.m. system. 


There is good reason to believe that high-energy experiments are completely 
in agreement with the hypothesis of isotopic invariance of nuclear forces. The 
isotopic spin formalism allows us to make some progress in an analysis of the 
results of these experiments. 

As indicated above, we can assign an isotopic spin of !/2 to a nucleon. The 
proton corresponds to an isotopic spin projection of +1/, while the neutron 
corresponds to a projection —'/,. A system of two nucleons can have an iso- 
topic spin of either unity or zero. In a system consisting of two protons the 
isotopic spin is unity (the projection is +1). An n-p system may be found ina 
state with isotopic spin of either 1 or 0. In states with T = 1 the interaction 
should be the same as that for the interaction of two protons (neglecting the 
Coulomb interaction). States with T = 0 do not have an analog in a system of 
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two protons and the interaction in these states applies only to two different 
particles. Hence, even within the framework of isotopic invariance one would 
not expect that the scattering of neutrons by protons would be the same as the 
scattering of protons by protons. Moreover the neutron cross section at a given 
angle should not be equal to the sum of the scattering cross sections for the states 
with T = O and T = 1 since these states interfere with each other. It might be 
possible, however, to find a method for identifying the interference terms, so 
that it would be possible to distinguish between these states. For this purpose, 
instead of using the neutron cross section at Y in the c.m. system we take the 
sum of the cross sections at 4 and r — Y. The physical significance of this pro- 
cedure is that the neutron cross section is replaced by the particle cross section, 
that is to say, We consider the total number of neutrons and protons scattered 
at a given angle. This procedure is a natural one in the framework of isotopic 
invariance, and the interference terms vanish. 
The quantity 


(9) = ony F) + Ong — 9) 


then represents the sum of the scattering cross sections for the states with T = 0 
and T = 1, 


(9) = o8) + olè). 


Since 0,(9) is nothing more than the proton-proton scattering cross section 
(neglecting the Coulomb scattering), knowing opp(®), from the experimental data 
we are now able to calculate ou(d). 

First of all we must note that the cross section vo(Y) calculated in this way is 
always positive. 

In Fig. 5 we show the behavior of o0(9) and oı(d) at energies near 400 Mev. 
The fact which is most striking is that these curves are completely different. 
The oo(#) curve is very much like the curve given by perturbation theory. 
Although the accuracy for this curve is small (since it represents the differences 
between two curves, each of which is characterized by an error of the order of 
10%), it is reasonable to assume that the interaction in states with T = 0 is 
weak and can be described, to some degree at least, by perturbation theory. 

To close this brief survey of nuclear forces we would like to call attention to 
the fact that the interactions between protons and neutrons in a nucleus may 
be different from thc interactions of free particles which are studied in scattering 
cxperiments. In ordinary atomic systems the electronic interaction is made up 
of the mutual interactions of all particle pairs. For example, the potential 
energy associated with three electrons is given by the function 


ee (7 Bee de aon): 
nsr) [re — rol lrs — ri] 


i.e., the interaction of the threc particles is equal to the sum of the mutual inter- 
actions of the thrce particle pairs. This “additivity” of the Coulomb forces is a 
consequence of the linear equations of electrodynamics. We have no basis for 


Nuclear Forces 31 


believing that nuclear forces are also additive; the interaction of two nucleons 
may not be independent of the presence of a third nucleon in the immediate 
vicinity. It is possible that nuclear forces are not additive and that the inter- 
action of particles in the nucleus is not simply the sum of the pair interactions. 
If such is the case the study of complex nuclei will lead to information on nuclear 
forces which cannot be obtained from investigations of two-particle systems. 


LECTURE FOUR 


Nuclear Structure 


(INDEPENDENT PARTICLE MODEL) 


I t is clear that no simple model can represent all the proper- 
ties of an extremely complicated quantum system such as the nucleus. Any 
model, of necessity, must have limited application. We should not be surprised if 
different effects require different models for their description; sometimes these 
models may even have mutually exclusive properties (for example, the inde- 
pendent particle model and the optical model). 

Many properties of the nucleus can be explained by the independent particle 
model. This model is based on an analogy with the electronic shells of atoms. 
The analysis of electrons in an atom is based on the assumption that each 
electron moves in some average field produced by all the other electrons. 

In a similar way we may use the notion that each nuclear particle moves in 
the average field produced by all the other nucleons in the nucleus without 
interacting with any of the other particles individually. 

We have already directed attention to the fact that the nucleon dimensions 
are of the order of 4.5°10714cm. The distance between nucleons in a nucleus is 
approximately 1.8-1071% cm. Thus, roughly speaking, the nucleons occupy 
1/59 of the volume of the nucleus. It is not surprising to find that nucleon 
properties are maintained inside the nucleus. In particular, this situation is 
responsible for the fact that the magnetic moments of nucleons inside nuclei are 
the same as for free nucleons. 

If, in accordance with the above, we describe the motion of each nucleon 
individually in a common field, it is most natural to take this field as spherically 
symmetric and to describe the states of the nuclear particles in terms of the well- 
known quantum-mechanical relations for the motion of particles in such a field. 
We will discuss the validity of this assumption later; in any case, this assumption 
certainly holds in light nuclei. 

In a field of spherical symmetry a particle has a definite angular momentum. 
The orbital momentum of the particle (in units of %) is denoted by the symbol 1 
(in contrast to the symbol L, which is used to denote the orbital momentum of 
the system). As is well known, the corresponding quantum number / can 
assume integral values 0, 1, 2,..., and these states are denoted, respectively, by 
the symbols s, p, d, f, and so on. 

33 
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In addition to its orbital momentum, a nucleon has spin. Thus, the total 
angular momentum of the nucleon j is composed of two vectors 


j=13+s. 


In the analysis of the interaction of two particles, we have seen that the spins 
of the nucleons interact with the orbital moment. This interaction can be 
described in terms of tensor forces, 1. e., forces whose potentials are proportional 
to the expression 


(s, n)(s: n) 
where n is a unit vector along the line connecting the particles. 


It is impossible to form an expression of this type for a nucleon which moves 
in a central field. The field itself has no spin and thus we have at our disposal 
only the nucleon spin s and its radius vector r or the corresponding unit vector 
n=r/r. Itisimpossible to form a scalar from these two vectors because, as we 
have seen, $n is a pseudoscalar while (s:n)? is equal to 1/4, since the projection 
of s in any direction (in particular in the direction of n) is +*/z. But this does 
not mean that the energy of the nucleon is completely independent of the spin 
direction. 

The spin of a nucleon in a nucleus is related to its orbital motion; the relation, 
however, depends on the velocity of the nucleon v. Hence, we can introduce 
one more vector—the nucleon velocity—in the energy expression. It is now 
possible to form a scalar from the three vectors s, n and v. This scalar will be 
of the form 


Jar X v)s 


where f(r) is a scalar function of the coordinates. Since 
1 
rxv=- _-1l, 
m 
the energy expression for the spin-orbit coupling can be written in the form 
1 
— f(r)l-s. 
m 


This interaction is relativistic (it vanishes when v — 0) and is weak compared 
to the interaction of the nucleon with the avcrage field of the nucleus. 

The experimental data indicate that the spin-orbit interaction is characterized 
by an energy of the order of 2 Mev whereas the interaction with the nuclear field 
is of the order of 20-30 Mev. 

In a system of several nucleons the spins of the different nucleons interact 
with each other as well as with the orbital moments. 

It is well known that in light atoms the spin-orbit interaction for a given 
electron 1s small compared with the interaction of electrons between themselves. 
In a system of this type the orbital moments of the particles are added, forming 
the total orbital momentum of the atom 
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The particle spins are added the same way since the coupling due to the identity 
of the particles is also stronger than the spin-orbit intcraction 


S= Š sy 
4 


If there is no spin-orbit interaction, in general the state of the system does not 
depend on the mutual orientation of the vectors L and S, that is, the vector sum 


J=L-+S. 


The weak interaction between the spins and the orbital motion means that 
different values of the quantum number J correspond to different states of the 
system. 

This coupling scheme is called the Russell-Saunders scheme. 

There is another type of coupling scheme; this scheme corresponds to the case 
in which the spin-orbit interaction of each individual particle is large compared 
with the interaction of the spins and orbits with those of other particles. In this 
scheme the total angular momentum of a particle is given by 


ji = l4; +s: 


and is unique since the coupling between the vectors 1 and s is not disturbed by 
the weak interaction with other particles. 

The jz vectors of the individual particles are added to form the total angular 
momentum of the system 


J=} j 


This scheme is called the 7-7 coupling scheme. It has been shown experi- 
mentally that this type of coupling obtains in certain heavy atoms. In nuclei 
the 7-7 coupling scheme is much closer to the actual picture than the Russell- 
Saunders scheme. 

The notation of the j-7 coupling scheme is now described briefly. 

The state of each particle is given by its total momentum and parity. In 
general, the orbital moment is not known in the case of arbitrary spin. How- 
ever, in the casc of a particle with spin !/, the orbital momentum is uniquely 
determined. It is easy to show that the vector j can be formed either from a 
state with / = 7 — 1/2 or from a state with l = ¡+ 1/2. Fora given value (half 
integral!) each state has different parity; hence the parity and 7 determine the 
quantum number /. 

It is obvious that the following states can exist (with different ¡ and J): si, 
Pr Paja, dija, ds/, andsoon. Thes,d,... states are even, while the p, f,. . .states 
are odd. Several different states of the same type, that is, several states with 
the same l and j, can exist in the nucleus. They differ in terms of the numbering 
sequence. There are two notation systems. One of these follows the example 
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of the principal quantum number in the hydrogen atom* and starts the enumer- 
ation of states with n = ¿+ 1. In this scheme the first s-level is denoted by Is, 
the first p-level by 22, the first d-level by 3d, and so on. 

In the other scheme (which we will use) the numbering of the levels starts 
with unity. 

Because we have so little information about the average nuclear field in which 
the nucleons move we cannot determine theoretically the order and disposition 
of the levels. This information can be obtained only by experiment. Hence 
an analysis of the nuclear properties can only establish a number of gross features 
of the arrangement of the levels. 

It turns out that the level energy increases with increasing orbital quantum 
number. This rule is a result of the fact that the centrifugal potential of the 
particle increases with increasing 7, and this leads to a reduction in binding 
energy. 

Further, the spin-orbit coupling in the nucleus turns out to be such that the 
level with 7 = / + 1/2 (that is, 1 parallel to s) is lower (has greater binding 
energy) than the level with j = / — 1/2. There does not seem to be any excep- 
tion to this rule. 

The following empirical rule pertains to the isotopic spin of nuclei. 

We recall that the projection of the isotopic spin of the nucleus Tr, defined as 
the sum of the projections of the isotopic spins of the neutrons and protons 


Lp 2 Tit = 1/(Z — N), 


is equal to half the neutron excess of the nucleus with sign reversed. 

It is obvious that the projection Ty can be associated with any isotopic spin 
vector greater than IT |; that is, a nucleus with an isotopic spin projection Ty 
can be characterized by states with isotopic spins which satisfy the inequality 


T > |Tel. 


Data on light nuclei (in which the isotopic spin is a good quantum number be- 
cause of the small electric interaction) lead to the conclusion that the ground 
state of the nucleus is characterized by the smallest isotopic spin compatible 
with this inequality, that is to say, the isotopic spin of light nuclei in the ground 
state is given by 


T ground = Tel = 1/(N — 2). 


This rule is clearly connected with the nature of the interaction between the 
neutron and proton. We have seen that in an n-p system, of the two possible 
values of the isotopic spin, T = 1 (the state characteristic of two protons) and 
T = 0, the T = 0 state (deuteron) has a larger binding energy. (In the latter 
state the attraction forces do not bring about a real level.) This rule also 
applies to complex nuclei; for example, the Be” nucleus and the CY nucleus have 


* The principal quantum number in the hydrogen atom is related to the energy of 
the level by En = —(R,/2n?). 
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ground states with isotopic spin T = 1. The BY nucleus has the same number 
of neutrons and protons and its isotopic spin in the ground state turns out to be 
T = 0. However, the level T = 1 in this nucleus, which is analogous to the 
level in Be and C*, lies 1.74 Mev above the ground state. 

It is possible to formulate a number of rules pertaining to ground state spins 
of nuclei. These rules determine the manner in which the angular momenta 
of the individual nucleons j are added to form the total moment of the nucleus. * 

If the nucleus has an even number of protons and neutrons, that is, if the mass 
number A and the atomic number Z are both even, the spins add in such a way 
as to make the total momentum of the system zero 


> jı = 0 (for nuclei with even N and Z). 
$ 


Furthermore, if all the nucleons (beyond the closed shells) are in the same state, 
in nuclei in which A is odd, that is, in which either the number of protons or the 
number of neutrons is odd, the spins add in such a way as to make the total 
momentum of the system equal to the momentum of one particle. 

For a long time this feature of nuclei has been used as a basis for the incorrect 
model of odd nuclei in which one (odd) particle moves in the field of all the other 
particles, which form a closed core. 

Finally, in nuclei in which the number of protons and the number of neutrons 
are both odd, that is, in which A is even while Z is odd, and in which the neutrons 
and protons are in the same states (the same j and parity) the total momentum 
is equal to twice the momentum of one nucleon. As is well known there are only 
four such stable nuclei. These are H?, Lif, Be! and N. Other nuclei of this 
type are radioactive. 

The spin is connected with the magnetic moment of the nucleus. From 
analyses of experiments on atomic spectra it is well known that atomic magnetic 
moments (for example, those of the rare earths) are fairly well explained by the 
Russell-Saunders vector model. It is of interest to determine the degree to 
which the independent-particle vector model is able to explain the magnetic 
moments of light nuclei. 

In determining the magnetic moment we may recall that the magnetic moment 
associated with a particle in the nucleus is made up of two parts: the inherent 
moment of the nucleon, parallel or anti-parallel to the spin of the nucleon, and 
the orbital moment, produced by the motion of the charged particle in its orbit— 
obviously, this part of the magnetic moment can arise only from protons. 

In accordance with well-known formulas of electrodynamics the orbital mo- 


ment of the proton 1s: 





e e 
=f v= r Xp, 
re E 


* At all times we shall be considering nucleons outside of closed shells, because the 
nucleons which form these shells make no contribution to the spin and magnetic 
moment of the nucleus. 
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where v is the velocity of the proton, mp is its mass and pits momentum. Since 
r Xp = 1 is the angular momentum of the proton, 








_ en i 
Ag 2mp 
The factor 
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ENT 2M pC 


is called the nuclear magneton; it is smaller than the Bohr magneton 


= 


Lo = 
2mc 


(m is the mass of the electron) by a factor m,/m and is equal to 5.05-107*4 
erg: gauss !, 

Nuclear magnetic moments are usually measured in nuclear magnetons so 
that in the following we will drop the factor uy from the formulas. 

The spin magnetic moments of the proton and neutron have been measured 
experimentally. The magnetic moment of the free neutron is —1.91, while the 
magnetic moment of the free proton is 2.79. 

In speaking of measurements of the magnetic moment of a nucleus one 
always understands the value of the magnetic moment averaged over the motion 
of particles in the nucleus. This mean magnetic moment of a nucleus is directed 
along the spin direction since this is the only preferred direction in the nucleus. 
Thus, we may write 


(u) = gij; 
the factor g; is called the gyromagnetic ratio. We may note that the magnetic 


moment usually given in tables is the maximum value of the projection of (#) 
in the direction of the magnetic field, 1. e., the quantity 


u = gij. 
The angular momentum of a nucleon j is made up of two parts 
jolts, 
and the magnetic moment of a nucleon can also be written as the sum 
u = gil + 88. 


In this expression gı and gs are called, respectively, the orbital gyromagnetic 
ratio and the spin gyromagnetic ratio. It is apparent that these factors are 
equal for the proton and neutron. Specifically, assuming that the inherent 
value of s is +1/,, we have 


forthe proton g, = land 
g. = 5.58, 
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for the neutron g; = Oand 
gs = —3.82. 


Our problem is to determine gy for the proton and neutron, expressing this 
quantity in terms of gz, gs and 7 (or J). 
We rewrite the expression for ¥ in the form 


u = '/(g + gl +s) + 1/2: — gl — 8). 


The average value of u must be along j. In order to determine the gyromagnetic 
ratio we multiply this expression by j = 1 + s. To compute the scalar product 
u : j we take the dot-product of the vector Y and the vector j. Thus 


Uj = 21? = 2379 +1) = 2G + 1); 


here we have used the definition 4 = g;j and the identity j? = j(j7 + 1). 
Using these expressions we find 


1 1 I—sXl+s+1 
g= zlate) +5 (er ae TY, 


HIT) 
Substituting s = 1/2 andj = / + !/,, we have 
2s — gı 
. == + |, 
gi £1 ol + 1 


For the proton this expression yields 


288 (l 
2.29 ; 1 
n= (1 +); = 7 + 2.29 (i -1 +3) 


while for the neutron we have 


191 e :) 
ee, J 9 , 
1.91 1 
„= -j = -191 G -1+2). 
J 2 


In Table 3 we present the values of g; and u = gy] for the proton and neutron 
in various states. 


TABLE 3 
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In computing the magnetic moment of a nucleus we must take the vector sum 
of the magnetic moments of all the nucleons and project it in the direction of the 
nuclear spin because, as in the case of a single nucleon, experimentally we can 
determine only the mean value of the nuclear magnetic moment. In this case 
nucleons in closed shells contribute nothing to the total magnetic moment and 
in the calculation we take account only of nucleons in unfilled levels of the 
nucleus. 

If there is only an odd number of neutrons with the same 7 beyond the closed 
shells, 


(W )aucı. = > Zij = gi > j = gjj 


by virtue of the addition rule for nucleon spins given above. The same relation 
applies for nuclei in which there are only protons beyond closed shells. 

In these cases, both the spin and the magnetic moment can be described in 
terms of a model consisting of a single particle which moves in the field of the 
remaining nucleons. However, such a model is no longer valid in nuclei in 
which there are both neutrons and protons beyond the closed shells. In this 
case it is impossible to express the magnetic moment of a nucleus in terms of the 
spin since the gyromagnetic ratios of the neutron and proton are not the same. 
We shall carry out a calculation of the moments of a number of light nuclei in a 
later lecture. 

Another characteristic electric property of the nucleus is the quadrupole 
moment. The quadrupole moment of a system of charges is given by the tensor 


Qik = a Cal SX igX key — di? a?), 


A 


where the summation extends over all particles (denoted by the index a). If the 
charged particles have the same charge e, this quantity can be taken out from 
the summation sign. 

Usually the quadrupole moment is measured in barns (1072 cm?). In de- 
fining the quantity in this way we neglect the factor e and measure distance in 
units of 1071? cm. 

From the definition of the quadrupole moment it is apparent that the sign of 
the component Qzz for a nucleus 


Qa = > (324° = Ya”) 


(summed over all protons) is determined by the shape of the charged body. 
If the body is elongated along the z-axis, Ozz > O (since the mean value of the z- 
coordinate is larger than the mean value of the square of the radius); if the body 
is compressed along this axis, Oz < 0. 

In a quantized system, the quadrupole moment, like the magnetic moment, 
is related to the spin vector. To find the relation we make use of the following 
properties of the tensor Q:x a) it is symmetrical Qix = Qx, b) the sum of the 
diagonal elements is zero Qrz + Qu + Oz. = 0. These properties allow us to 
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express Qix uniquely in terms of the spin vector (quantized) of the nucleus, I. 
Specifically, 


Qir = A(IiỌlr + Lil; — 2/P8:5x), 


where A is aconstant. (We may recall that the components of the spin vector 
do not commute.) 

The quadrupole moment Q given in tables refers to the component Oz, of this 
tensor in the state in which the projection of the spin on the z-axis is equal to 
the spin itself. Setting J, = Z and I? = /(I + 1), we have 


O = 2%AIQI — 1). 


It is apparent from this formula that when J = 0 or when J = ?!/, the nucleus 
cannot have a quadrupole moment, since Q = 0 identically. 

Using this same formula we can express the constant A in terms of Q and, 
making use of the value which is obtained, write the expression for the quad- 
rupole moment tensor in the form 
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Qu = sar — 1) 


E + Luli — = II + Dea |: 


LECTURE FIVE 


Structure of the Nucleus 


(LIGHT NUCLEI) 


W. turn now to an analysis of the magnetic moments of nuclei, 
starting with light nuclei. We shall see that the measured values of the mag- 
netic Moments are in agreement with shell theory up to nuclei consisting of 
more than 20 particles. 

n, p. We start with individual nucleons. The magnetic moments, measured 
in nuclear magnetons, are 2.79 for the proton and —1.91 for the neutron. It is 
obvious that speaking of states in one nucleon is meaningless; however, for 
consistency of notation we may speak of a single particle as being in a 1s state. 

d. The deuteron is primarily in a state with moment L = 0. We will speak 
of the two particles as being in the 1s? state, using the 1s? designation for this 
state by analogy with atomic nomenclature for particles in a spherically sym- 
metric field. 

The magnetic moment for the 1s? configuration should be equal to the sum of 
the moments of the two nucleons, i. e., 2.79 — 1.91 = 0.88. Actually the mag- 
netic moment of the deuteron is 0.86. This may be taken as good agreement. 

H*. The three nucleons which make up tritium can all be in 1s levels without 
contradicting the Pauli principle, since only two of them are identical. This 
structure results in a spin of !/2 for tritium. Thus, the H? configuration may be 
designated by 1s*. The magnetic moment of this configuration is easily deter- 
mined if one notes that there is one unfilled proton 1s state. It is clear that a 
system consisting of tritium plus a proton, in which this state is filled, will not 
have a magnetic moment. Hence, a system in a 1s* configuration (one “hole” 
in the 1s shell) should have the same moment as the proton.* Again referring 
to this property the “hole” state is also sometimes designated by 1s7?*. Thus, 
the theoretical value of the tritium magnetic moment is 2.79. The measured 
magnetic moment is 2.98. The discrepancy between theory and experiment is 
0.2 magneton. Below we shall see that the discrepancy of 0.2 magneton 1s 
characteristic of the accuracy of the shell model with 77-coupling in magnetic 
moment problems. 

He’. The next nucleus is the mirror image of tritium. Its spin is also 1/2. 
Here we are dealing with the 1s? configuration with the neutron site unfilled. 


* The sign is the same since the sign of the magnetic moment is defined with 
respect to the spin direction. 
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Hence, the magnetic moment should be equal to the magnetic moment of the 
neutron, —1.91. Experimentally the value is found to be —2.13. As in trı- 
tium, the discrepancy is approximately 0.2 magneton, except that in this case the 
theoretical value is larger than the experimental value. 

He‘. The spin and magnetic moment are both zero. All four 1s states are 
filled; thus we arrive at the first closed shell (15%). 


He? and Li®, These nuclei do not exist in nature. They are unstable against 
separation into He‘ and a nucleon. This is evidence of the fact that the 1s* 
shell, which becomes filled at Het, is a closed system and that an additional 
nucleon cannot be added. 

Data for proton and neutron scattering on Het reveal a resonance associated 
with the intermediate nucleus He? (or Li), and indicate that the nucleon which 
is not part of the He? occupies a 1P»,, state. 


Li‘. The spin of this nucleus is 1. If, as might be expected from the simple 
scheme for filling shells, Lif had the configuration He*1p3,,?, according to the 
empirical rule for forming moments in odd-odd nuclei its spin should be 3. This 
spin is actually observed in B*, which, as we shall see later, has the configuration 
1ps/,7?. It may be assumed that this empirical rule is not obeyed; in the Li? 
nucleus the spins of the two nucleons add up in such a way as to form a total 
spin of 1. 

This system could have a magnetic moment of 0.63, whereas the experimental 
value is (Lif) = 0.82. This discrepancy is not very great; it is 0.2 magneton, 
the discrepancy already noted for H? and He?. However, a marked discrepancy 
arises when the quadrupole moment of the system is calculated. Itis well known 
that the quadrupole moment of Li® is no larger than 2% of the quadrupole 
moment of Li’ (—4.6-1074 barns). This means that the quadrupole moment of 
Li® is essentially zero. On the other hand, two nucleons in a 19s,,? state should 
have a large quadrupole moment. It is apparent that in the Li” nucleus we are 
dealing with an irregularity in the system for filling levels. This effect, which is 
frequently considered anomalous, is actually not strange. The same effect is 
encountered in electron shells in atoms. Here, the shells also filled in an irregular 
manner. It is well known, for example, that in the rare-earth region there is 
“competition” between the 4f, 6s and 5d shells, and that in the uranium region 
the 5f, 7s and 6d shells are filled together. Hence, one should not expect absolute 
regularity in the order of filling nuclear levels. The simplest hypothesis as to 
the structure of the Li? nucleus is based on the assumption that both nucleons 
(for brevity we will frequently not bother to write the nucleons in the closed 
shells) occupy the new s-state, 251/ The configuration 2sı,,? for Li® is in agree- 
ment with the empirical rule of spin 1. The quadrupole moment of this state is 
obviously zero (spherical symmetry of the s-state!). The magnetic moment, as 
in the deuteron, equals the sum of the magnetic moments of the proton and 
neutron (the value is 0.88), and is also in good agreement with the experimental 
value. 

If both nucleons were in the 11, state rather than 2sı,,, this would also give 
no quadrupole moment (the nucleon moment is 1/2). This Lif structure, how- 
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ever, would disturb the spin-orbit level scheme, and gives an incorrect value for 
the magnetic moment, namely 0.40. 

Li’. The spin of this nucleus is 3/2. The magnetic moment is 3.26. The 
nucleus consists of 3 nucleons—1 proton and 2 neutrons (beyond He‘). These 
three particles should be in a 1p:⁄} configuration, as is indicated by the spin 
value. We cannot use the simple considerations which have been used up to this 
point to compute the magnetic moment of this system. The wave function of a 
system of 3 particles is not unique even when the spin of the system is assigned. 
To determine the state uniquely it is necessary to determine the isotopic spin of 
the system. A general empirical rule indicatcs that the lowest state has the 
minimum isotopic spin. Hence, we assign to the ground state of Li’ an isotopic 
spin of 1/2. The magnetic moment computed for this state (19,3, T = Ya. 
I = 4/,) is found to be 3.07, which is 0.2 smaller than the experimental value. 
These results seem to verify the validity of the above scheme. 

Thus, in the Li’ nucleus we start filling the p3,, shell; this shell continues to be 
filled in subsequent nuclei. 

Be’. The spin of this nucleus is */2 and we shall assume that all 5 particles 
are in 19:,, states and that the spins add up in accordance with the rule which 
states that the total spin is equal to the moment of one particle. This con- 
figuration can also be written as a configuration in which three particles (2 
protons and 1 neutron) are lacking for filling the C*? shell. Thus, the Be? 
configuration can be designated either by He*lp3,, or by C!21p:,”?. The 
magnetic moment of this system also depends on its isotopic spin. Assuming 
that the isotopic spin of Be? is 1/2, we compute the magnetic moment. It turns 
out to be —1.14. The experimental value is —1.18. The agreement is rather 
good. 

The Be? nucleus is interesting in another respect. Using this nucleus, we 
can show the shortcomings of another model which is frequently discussed in 
the literature—namely, the a-particle model. From the point of view of the a- 
particle model the Be? nucleus consists of 2 a-particles and 1 neutron. The 
magnetic moment due to 1 neutron in a state with j = °/2 is either 1.15 (if this is 
the ds, state) or — 1.91 (if this is a psy, state). Both of these values are in dis- 
agreement with experiment. 

BY The spin of B is 3; this value is in agreement with the rules and 
corresponds to 2 nucleons of different types in the 1p:/, state. In the present 
case we are discussing unfilled configurations, that is, C*?1p3,,7?. The calculated 
moment of this state is 1.88 and the experimental value is 1.80. 

Bu The spin of B" is 3/2. We would expect a C!?1P:,,”! configuration for 
this nucleus and a corresponding magnetic moment, equal to the magnetic 
moment of the proton in the p3,, state. According to Table 3 of the previous 
lecture, this value is 3.79 and is in sharp disagreement with the experimental 
value 2.69. Thus, we again encounter an irregularity in filling the shells. It is 
reasonable to assume that the nucleons can avail themselves of the 2s state. 
Agreement with experiment can be reached only by assuming that 2 nucleons go 
into the new 2sı,, state while the remaining 5 remain in the 1Ps,, state. Un- 
fortunately, we cannot compare theory and experiment since the level is not as- 
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signed uniquely. The ambiguity is due to the fact that, although the total 
isotopic spin of the system is given, it is impossible to ascertain the way in which 
it is distributed between the two groups of particles. We can say only that the 
proposed scheme does not violate the magnitude of the quadrupole moment of 
this nucleus. Thus the following configuration can be assigned for B*!: 


He‘1pz/,° 251/72. 


C12, This nucleus has neither spin nor moment, and consists of two closed 
shells, 1sı,,* and 1p3,2. The filling of the 1:,, shell is not characterized by the 
same strong reduction in the binding energy as is the case for the Isı,, shell in 
He‘. 

An idea of the degree to which a shell is closed can be obtained from the 
nucleon affinity in a nucleus, i. e., the energy evolved in adding a proton or 
neutron to this nucleus. It is found that the nucleon affinity in a nucleus with a 
closed shell is considerably smaller than in neighboring nuclei. With an increase 
in atomic weight, however, this effect is reduced because of irregularities in 
filling of the shells and the fact that all shell schemes become less reliable as the 
number of nucleons in the nucleus is increased. In He‘ the affinity is negative, 
since the nuclei He? and Li? do not exist. In C!? the neutron binding energy 
is 4.9 Mev, while that of the protonis 1.9 Mev. Inthe B" nucleus the neutron 
and proton binding energies are, respectively, 11 and 9 Mev. 

C13, The spin of the CY nucleus is 1/2. This spin can be obtained if we add a 
nucleon in the 1ġ:ı/, state or in the 2s1,, state. In the first case the magnetic 
moment is 0.64 (cf. Table 3); in the second case it is —1.91. The experimental 
value is 0.70. The result determines the state uniquely as C!?1pı,,. 

N! The N! nucleus corroborates the finding that the 1p:,, state begins to 
be filled after C*?. The spin of N!‘ is 1, corresponding to a neutron and proton 
in the state j = !/}. The magnetic moment of this state is 0.40, a value which is 
in agreement with the experimental value. Thus, there is no doubt that the 
configuration for this nucleus is C!71p1/,?. 

At this point we should note the extremely surprising long lifetime for the f- 
active nucleus CM, Reasoning by analogy with N, this nucleus should have 
spin 0, corresponding to a configuration C!*1p:/,? and the decay reaction C! > 
N1 should be an allowed one. Actually, this reaction is strongly forbidden (the 
half-life 1s 5600 years, instead of several hours as is typical of an allowed decay). 
This anomaly indicates that there is a difference in the shell structure in C' and 
N, 

N'. The N! nucleus furnishes new evidence that the 1Pı,, level becomes 
filled. The spin is !/2 and the magnetic moment is —0.28. If this configuration 
is interpreted as one “hole” in a 1f1,, configuration (lacking one proton for O18): 
O161p1/,7*, the moment should be — 0.24. 

O*. In this nucleus the 1p»,, shell is filled. The structure is 


1s1/,1p3/21p1/,. 


Like the C*, it is closed. The binding energy for a neutron in O! is 4.1 Mev; 
the binding energy for a proton is 0.6 Mev. This quantity is considerably 
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smaller than, say, in the O” nucleus, in which the binding energy of a neutron is 
7.9 Mev and that of a proton, 5.7 Mev. 

O", The spin of the O" nucleus is 5/2}. This means that a shell with / = 2 
starts to be filled and we are encountering the 1ds/, state for the first time. The 
neutron in the ds/, state has a magnetic moment of —1.91, which is in excellent 
agreement with the value — 1.89 measured experimentally. 

F*% Contrary to our simple assumption, the spin of F!? is not 5/2, as would 
be the case if the 1ds,, state were filled; rather, it is *'/z. This value apparently 
means that 2sı,, states enter, displacing the 1ds,, state. 

The magnetic moment of the system Ne?(2s./,)~! should be equal to the 
moment of the proton, 2.79. Experimentally the value is found to be 2.63. 

Ne”. In accordance with what has been said above this nucleus consists of 
four closed configurations: 

1si/,4 1ps/2 1p1/,* 281/4. 

Ne?!. The spin of this nucleus is 3/s, indicating a 2p3,, configuration. In this 
case (Ne? plus a p3,, neutron) the magnetic moment should be —1.91. There 
are indications that the magnetic moment of Ne?! is negative. Further study 
of this nucleus will be of great interest. 

Na??. This is a radioactive nucleus in which both the spin (3) and the mag- 
netic moment (1.75) are known. It is most reasonable to assume that the two 
nucleons (beyond Ne”) are in 23/, states. Thus, the Na?? configuration is very 
similar to a B! configuration, in which there are two “holes” in the 1ps,, level. 
The spin of this system should be 3, and the magnetic moment, as we have seen 
in the B” case, 1.88. We may recall that the magnetic moment of Na? (experi- 
mental) is 1.79, that is, almost equal to the B® moment; these results tend to 
support our picture of the structure of these two nuclei. 

Mg”. This is the last nucleus whose structure and configuration can be 
analyzed without introducing a large number of assumptions. The spin is 5/2. 
The magnetic moment is —0.85. These quantities are reasonably well explained 
if we assume a configuration Ods}? = Si*%ds,, ~? in Mg®. With T = 1%, this 
configuration has a magnetic moment which agrees with experiment. It is 
curious that in this case even those 251/,4 levels which are being filled seem to 
vanish. 


Further comparison of the theoretical values of the magnetic moments with 
the experimental values is not fruitful because the large numbers of particles in 
the shells mean that the computations must be ambiguous. However, the 
success of the independent particle model in describing light nuclei is apparent 
even from the cases which have been considered. “To make these considerations 
clearer we have collected all the data in Table 4, in which the values of the spins 
and magnetic moments of the light nuclei are tabulated. 

The independent particle scheme makes it possible to understand the magnetic 
moment and spins of the light nuclei and to develop a number of characteristics 
of those nuclear reactions in which excited states and ground states of the 
nuclei are both of importance. A direct measurement of the spin and magnetic 
moment of excited states of nuclei is beyond the capabilities of present-day 
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TABLE 4 





Magnetic 


Spin Magnetic 
Nucleus moment Configuration at 


Measured values (theoretical) 


n 1/2 — 1.91 1517, 

H! 1/2 2.79 lsı), 

H? 1 0.86 Isiy,? 0.88 
H3 Ua 2.98 lsı/,? 2.79 
H? a —2.13 ls1/,? —1.91 
He Closed shell 1s1/,' 

Li® | 1 0.82 Het 2s1/,? 0.88 
Li? | 3), 3.26 Het 1pz/,3 3.07 
Be? | 3/5 —1.17 He! 1p:/,® —1.14 
B10 3 1.80 Het 1p3/,5 1.88 
BH 3/2 2.69 Het 1p3,,% 251/,? Not unique 
E | Closed shell lsı/,* 1 pay, 

ES 1/2 0.70 C2 Ipi 0.64 
N14 1 0.40 Cl2p17,? 0.40 
N16 ‘My —0.28 Cl pi," —().24 
O's Closed shell 1s1/,* Lpay,? 1p1/,* 

Or 5), —1.89 Old; /, —1.91 
p19 1/2 2.63 O12517/,° 2.79 
Ne” Closed shell si 193,8 TP 2817," 

Na2 | 3/3 <0 l Ne 237, —1.91 
Na?? 3 1.75 | Ne? 2p3/,? 1.88 
Mg? 5/a —0.85 | Ot 1ds/,° —1.06 
Siz | Closed shell lsi Ipar 1p1/,* 1d,/,*? 


experimental techniques; hence, the identification of these states must be carried 
out from studies of the selection rules and the angular distributions of the 
various particles. Analyses of this type have already yielded a large amount of 
information as to the levels in light nuclei, and this information is in agreement 
with the configuration scheme considered above. We present several examples 
below. 

The light nuclei H?, H? and He? have no low-lying excited states and the 
particle interaction in these nuclei is so small that the “potential wells” have 
only one level, corresponding to the ground state. 

Similarly, there are no levels in a-particles. However, there is evidence that 
an unstable level does exist in a-particles; this level can decay into tritium plus 
a neutron. It is believed that this level appears in the reaction 


H? + p— Hei + n. 


A very wide resonance is observed in this reaction at an energy of 2.6 Mev (in 
the c.m. system). This resonance can be interpreted as the appearance of a 
level with an energy of approximately 22 Mev and a width of approximately 
1 Mev in the a-particle. Such a level arises if one of the nucleons in the a- 
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particle makes a transition from the Isı,, state to the 12;,, state. In this case 
the total spin of the system is determined by the vector sum 1/2 + 3/2, and is 
equal to either 2 or 1 (states which are obviously odd). It is possible that an 
experiment can reveal the two states 2 and 1 separated by approximately 1 Mev. 
We may note that the absence of a stable level in the a-particle makes it possible 
to draw certain conclusions as to the stability of H*. One proton and three 
neutrons cannot all be in the 1sı,, state because of the Pauli principle. Hence, 
11 H‘ exists, one of the neutrons must be in the 1:/, state, and the H* nucleus 
must be similar to an excited state of Het However, there are no excited states 
in He? which are stable against emission of a nucleon; consequently, there are 
none in the H‘system. It may also be noted that there is no guarantee that the 
excited state of He? has an analog in the H? system. The H* system has an 
isotopic spin of at least 1 (projection 7x = 1). The excited state of He? can 
also have an isotopic spin of zero. 

An interesting effect is observed in the He? system. It has already been noted 
that in scattering of neutrons by nuclei there is a resonance corresponding to the 
formation of He! in the intermediate state. This level is unstable against decay 
to Het + p (it has an energy of 1 Mev), and corresponds to a 1pz/, neutron con- 
figuration. 

The same pattern is observed in the scattering of protons on helium, which 
leads to the production of the Li? nucleus in the intermediate state. In this case, 
because of Coulomb forces, the level lies somewhat higher than in He? (1.8 Mev). 

There is one other level in He? which is of interest. It is well known that the 
reaction 


H? + H? > Het + n 


has a very sharp resonance at 15 kev (in the c.m. system). The value of this 
maximum corresponds to the maximum possible cross section in a system of spin 
3/2, The surprising thing about this reaction is that the resonance is very 
narrow (width of approximately 40 kev) in spite of the fact that the neutrons in 
this reaction have high energy. The narrowness of this resonance indicates a 
comparatively high stability for the nucleus at this high level (16.6 Mev above 
the ground state of He®), which is strange. Actually, it is easy to explain this 
reaction if it is assumed that two particles in the He? system are excited to the 
1p3,, level. Since the excitation of He‘ requires an energy of approximately 20 
Mev, the appearance of a level at this height is completely reasonable. On the 
other hand, the decay of this excited He? nucleus into an a-particle and neutron 
should be characterized by transitions of two particles, one neutron—the 
emitted free neutron (continuous spectrum)—and another nucleon, to the 151, 
ground state. It is reasonable that a transition of this type will have a consider- 
ably smaller probability than the single-particle transitions which are encoun- 
tered in the majority of reactions. This situation explains the relative stability 
of this state. 

A similar effect is observed in the Li® system which is produced in the reaction 


He? + He? — He? + H!. 


50 Lectures on Nuclear Theory 


Here the resonance lies at a somewhat higher energy (260 kev) and, correspond- 
ingly, has a smaller magnitude. [As is well known, the maximum cross section 
is proportional to X? ~ 1/E.] 

In interpreting nuclear reactions one must also take account of the selection 
rules. The best known example of this kind appears in the decay of the Be? 
nucleus into two e-particles. As is well known, the two a-particles can be only 
in symmetric states (as for any two identical particles with no spin). This 
requirement means that they can separate only in s, d, or f states (1 even). 
Hence, the Be® nucleus can decay into two a-particles only from those states 
which have even spin and positive parity (isotopic spin zero). This selection 
rule is apparent in the Li” + p reaction. This reaction can occur with the 
emission of two a-particles, 


Li? + p > 2a, 


or, by the above reasoning, with the emission of y-quanta. In this case the 
reaction occurs in two stages: 


Li? + p > Bel + y. 


—2a. 


An interesting feature of this reaction is the fact that the y-ray yield curve has a 
maximum at a proton energy of 441 kev (which can be associated with the Be? 
level), whereas the «-particle yield curve has no maximum at this energy. This 
means that the Be? nucleus, in the state which is formed at this proton energy, 
cannot divide into two a-particles. This forbiddenness is explained by the fact 
that the Be? level in this reaction (excitation energy of 17.63 Mev) has spin 1 and 
is odd. This level cannot decay into two a-particles, and the excitation can be 
dissipated only by radiation. 

The isotopic spin selection rule appears in an interesting way in reactions. 
If, for example, we consider the inelastic scattering of deuterons or a-particles 
on nuclei, in which the nucleus is left in an excited state, the isotopic spin is zero 
for both the deuteron and a-particle, and we can excite only those levels which 
have the same isotopic spin as the ground state. For example, in the elastic 
scattering of deuterons on B" we excite levels in B* with energies of 0.72, 2.15 
and 3.58 Mev. In inelastic proton scattering, however, in addition to these 
levels we excite the level with an energy 1.74 Mev. Since the isotopic spin of 
BY is zero, it follows unambiguously that the isotopic spin of the level at 1.74 
Mev must be unity. This is the level which is similar to the ground state of 
B” and CY, 

The simplest manifestation of isotopic invariance occurs in the mirror nuclei 
Be’-Li’, CN, OY%.-FY, and so on. In these nuclei each level of one nucleus 
corresponds to a level of another and thus their spectra differ only by the 
Coulomb field (and mass differences due to the difference in the mass of the 
proton and neutron). 

A comparison of isobaric nuclei is sometimes desirable for determining the 
isotopic spin of a level. For example, consider the nuclei N! and O'. The N$ 
nucleus has a projected isotopic spin Ty = —1, while the O% nucleus has a 
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projection Tè = 0. This is apparent from the fact that the neutron excess 
N — Zin N" is 2, while this same quantity is zero in O*®. Thus, the levels in N' 
can have isotopic spin T = 1, 2,..., while the O! nucleus can have only levels 
with T = 0. Alllevels of the N" nucleus should be encountered in the spectrum 
of the O nucleus, but not all the O* levels are included in the N' spectrum. 
Those levels with T = 0 will not appear in the N! spectrum. The transition 
from a level in the O! nucleus to a corresponding level in the N** nucleus is 
accomplished by replacing one proton by a neutron. Hence, the Nt? nucleus in 
this state will have a mass differing from the mass of the O}* nucleus in the 
corresponding state by the Coulomb energy of the missing proton minus the 
difference of the mass of the neutron and proton. In light nuclei this quantity is 
small.* It is equal to zero (more precisely, 15 kev) in the pair H3, He? and in- 
creases with increasing nuclear charge. A comparison of the spectra of neighbor- 
ing nuclei is also desirable for identification of the levels. It is possible, in 
particular, to state that O' in the ground state has an isotopic spin of zero (01 
in the ground state has a mass which is smaller than the N!$ mass). 

This type of analysis can also be carried out for other isobars. Supplementing 
this information with other data on nuclear reactions we are now in a position to 
describe a large number of the lowest excited states of light nuclei. 


* As is customary in nuclear physics, we speak of the mass of the neutral atom 
rather than the mass of the nuclei. 


LECTURE SIX 


Structure of the Nucleus 


(HEAVY NUCLEI) 


I n the light-nucleus region we have assumed that the neutron 
and proton are “isotopically” similar; in the hcavy-nucleus region, in which an 
important role is played by Coulomb forces, this assumption is no longer valid. 
Whereas the properties of light nuclei depend on the total number of nucleons in 
a shell, the properties of heavy nuclei depend on the number of neutrons, N, and 
the number of protons, Z, individually. 

The explicit dependence of nuclear properties on N and Z is manifest in the 
existence of so-called “magic numbers” and corresponding “magic nuclei.” 

Magic nuclei are nuclei in which either the number of protons or the number of 
neutrons corresponds to a filled shell. The properties of magic nuclei are rather 
different from those of other nuclei. 

It is now generally believed that the numbers 50, 82, 126 and 152 are magic 
numbers. Other magic numbers which are encountered in the literature refer 
either to extremely light nuclei, in which these numbers have little significance, 
or to numbers whose existence is open to doubt. 

The best example of a magic nucleus is the lead isotope of atomic weight 208. 
This nucleus is, so to speak, a doubly magic nucleus, since it contains a magic 
number of protons (82) and a magic number of neutrons (126). Pb? exhibits 
properties which are characteristic of a closed system. Whereas the first excited 
levels in ordinary heavy nuclei are never higher than 200 kev, the excitation of 
Pb? requires an energy of 2.6 Mev. Furthermore, the binding energy for an 
additional nucleon is anomalously small. Neutron capture is usually character- 
ized by the excitation of the resulting nucleus to about 7 Mev; neutron capture 
in Pb2% leads to an excitation of only 4 Mev. 

Probably the best way of determining the magic nuclei is an investigation of 
the affinity of nuclei for a-particles (i. e., the energy released in capture of an a- 
particle by a nucleus). The dependence of the affinity for a-particles on the 
number of nucleons in the nucleus is a more reliable characteristic than the 
affinity for neutrons or protons, which depends on the number of particles in an 
irregular way, changing discontinuously in going from even to odd elements. 
Studies of the systematics of a-active elements and reactions involving the emis- 
sion of a-particles are extremely useful in this connection. 
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Fig. 6. Spectrum of the nucleus »Hf!® with an even number of protons 
and neutrons (from A. Bohr, Rotational States of Atomic Nuclei, Copenhagen, 
1955). The energies corresponding to the formula for rotational levels are 
given in parentheses. The level 11+ which appears at the highest energy is 
not rotational in character. The arrows indicate observed electromagnetic 
transitions (E2 electric quadrupole, M3 magnetic octupole). It is apparent 
that the quadrupole transitions occur between rotational levels (radiation 
of a charged drop). It is also clear that the last level is not a rotational level. 


It is unfortunate that, having established the magic numbers, there is still very 
little we can say regarding the nucleon states. Many more experiments will be 
required before the pattern of states for heavy nuclei becomes in any way com- 
plete. The level scheme and the order of filling which frequently appear in the 
literature are based on a simplified model of the potential well andan idealized fill- 
ing system which, as is apparent from the light-nucleus examples, does not 
correspond to reality. 

Moreover, experiment indicates that in heavy nuclei the model of a particle 
moving in a central field requires essential modification because the original 
assumption as to spherical symmetry is no longer valid. 

It is undoubtedly true that magic nuclei are sphericalinshape. However, this 
shape is extremely unstable, especially for heavy nuclei; the Coulomb forces, 
which tend to elongate the nucleus, and the surface tension forces, which tend to 
reduce the surface in heavy nuclei, almost completely cancel each other. Hence, 
even under the effect of relatively weak forces the nucleus becomes a sphere. 

The instability of heavy nuclei is very well illustrated by nuclear fission. 
With comparatively small excitation (approximately 6 Mev) the U? nucleus, 
which is formed as a result of thermal neutron capture, becomes unstable and 
somewhat elongated; eventually this situation leads to its division into two 
fragments. In the nuclei U*%8, Pu?®, ete., the instability is even larger—these 
nuclei fission spontaneously, being found in the ground state. Because of this 
instability many heavy nuclei in the ground state are elongated rather than 
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Fig. 7. Spectrum of the nucleus 3Ta!%! with an odd number of protons and 
an even number of neutrons (same reference as Fig. 6). The spin of the 
ground state corresponds to a nucleon spin projection Q = %/z. The first 
rotational level has an energy which almost corresponds to the same moment 
of inertia as the neighboring nucleus „Hf!80. The second level is also 
described by the rotational level formula. The energy ratio Eu,,:Es,, = 
2.21 + 0.02 is in agreement with the theoretical value 2.22. The rotational 
nature of the levels is corroborated by the observed quadrupole transitions, 
which, however, compete here with the magnetic transitions—dipole 
transitions associated with a change of Q. 


spherical. The elongation of the nucleus is related to the interaction between 
the surface and the nucleons outside closed shells. Without an analysis of the 
experimental data it is impossible to say in which nuclei the elongation becomes 
sufficiently large to make the nucleus nonspherical. The light nuclei are un- 
doubtedly spherical—this observation is based on the success of shell theory, 
which is based on states in a field of spherical symmetry. From what has been 
indicated above it is clear that a particularly marked nonsphericity is to be 
expected in nuclei at the end of the periodic table. 

The basic ideas concerning nonspherical nuclei have been most completely 
described by A. Bohr.* We now examine the properties of nonspherical nuclei. 

The most important distinction between nonspherical and spherical nuclei is 
that the former can have rotational levels. According to the basic ideas of 
quantum mechanics, it is meaningless to discuss rotation of a spherical nucleus. 
A spherically symmetric system, 1. e., a system whose properties are independent 
of angle, cannot have a rotational energy spectrum. The concept of rotation in 
such a system is meaningless; similarly, it is meaningless to speak of rotation in 
a spherically symmetric nucleus. However, in an elongated nucleus the concept 
of rotation is meaningful. If, for simplicity, it is assumed that the nucleus is an 
elongated ellipsoid of rotation, the nucleus can rotate about an axis perpendicular 
to the axis of symmetry. Rotation about the symmetry axis of the nucleus, 
however, is meaningless for the reasons given above. 

A nonspherical nucleus is characterized by the state of the “quiescent” 
nucleus and by the moment of inertia about the axis perpendicular to the 
symmetry axis of the nucleus. 


* We refer the reader to the collection ‘‘Problems of Contemporary Physics’’ No. 9, 
Foreign Lit. Press, 1955, in which there are translations of the two basic papers on 
nonspherical nuclei. (Cf. also A. Bohr, Rotational States of Atomic Nuclei, Copen- 
hagen, 1955.) 
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The concept of angular momentum for the individual nucleons loses signifi- 
cance in the field of a nonspherical nucleus because this quantity is a good 
quantum number only in a spherically symmetric field. The nucleon state in 
this field must be classified by analogy with a diatomic moleculc. 

In a field of axial symmetry the quantity which is conserved is not the angular 
momentum, but rather its projection on the symmetry axis. The projection of 
the angular momentum of the nucleon is made up of the projection of its orbital 
moment and the projection of its spin; it is not meaningful to speak of the total 
moment of a nucleon in a field of axial symmetry. 

We shall designate the sum of the projections of the nucleon angular momenta 
on the nuclear axis by the symbol Q. 

Let us consider a rotating nucleus. Rotation of a nucleus is characterized by 
the angular momentum K which, as we have seen, has no component along the 
symmetry axis of the nucleus. If we use the symbol n to denote a unit vector 
along the symmetry axis of the nucleus, the total spin of the nucleus can be 
written in the form: 


I= Q2 +K. 


The vector K is perpendicular to the symmetry axis, and thus 
Ken = 0. 


From this it follows that I:n = Q, 1. e., the projection of the nuclear spin on the 
symmetry axis is always equal to 2. 

From the obvious consideration that the projection of a vector cannot exceed 
its length, it follows that for a given Q the possible values of a nuclear spin are 


TI = 2,2+ 1,02 + 2, ..., etc. 


Special attention is required for a nucleus in which Q = 0, that is, a nucleus 
with zero spin. As the theory indicates, such nuclei are similar to diatomic 
molecules with identical spinless nuclei; hence the angular momentum can 
only be an even number—0, 2, 4,.... 

A study of nuclear spectra with the purpose of identifying rotational levels 
leads to interesting results. The most striking rotational spectra appear in the 
heavy a-active nuclei. The rotational levels are best seen in even-even nuclei, 
in which the excitation of individual nucleons requires a comparatively large 
energy and the nucleon excitation spectrum is not superposed on the rotational 
spectrum. An example of such a spectrum is given in Fig. 6. 

It is also possible to separate out rotational levels in odd nuclei, using the 
features described below. In Fig. 7 we show the rotational spectrum in the odd 
nucleus 73 Ta*!, 

It should be kept in mind that experimentally it is difficult to separate rota- 
tional levels with high moments. However, at the present time, rotational 
levels with moments of 8 have been observed. 

The energy associated with the rotational levels of a nucleus is given by the 
well-known formula 
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where M is the angular momentum and J is the moment of inertia of the nucleus. 
In a quantized system M? = #77([ + 1), and the energy of a level with moment J 
is given by 


N? 
Er == ; 
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J, the moment of inertia of a nucleus, has little in common with the moment 
of inertia of a solid ellipsoid. In the case of spherical symmetry the latter is a 
finite quantity—the moment of inertia of a sphere. However, rotation of a 
spherical nucleus, as we have seen, is not physically meaningful, and the most 
important property of the moment of inertia J is that J — 0 as the nuclear shape 
becomes spherical. 

Thus, nuclear rotation cannot be described in terms of the rotation of a solid 
body. A theoretical analysis of a nonspherical nucleus indicates that rotation 
can be described approximately by equations which are similar to the hydro- 
dynamic equations of the potential motion of an ideal liquid in a rotating asym- 
metrical shell. However, the suitability of this approximation is still not com- 
pletely established. 
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Fig. 8. Probable scheme for the first levels of 9¿Am?1? identified by B-decay 
of curium (from the data of S. A. Baranov and K. N. Shlyagin). The 
nucleus 9,Am?*? has an odd number of protons as well as an odd number of 
neutrons. The spins of the levels arise as a result of addition and subtrac- 
tion of the projections of the moments of both nucleons, each of which is 5/2. 
It is apparent from the scheme that the addition of the projections is to be 
associated with the lower level. 


The solution of this problem 1s characterized by precisely the same features 
as those indicated above. The usual conditions of hydrodynamics require that 
the normal component of the velocity of a liquid be zero at a boundary defined 
by a fixed wall. If the wall is not fixed, but moves with a given velocity, the 
boundary conditions require that the normal (to the wall) component of the 
velocity be the same as the normal component of the wall velocity. 

We consider a liquid enclosed in a spherical container rotating about an axis 
through its center. It is apparent that in such motion the velocity of each 
point of the surface is tangential and that the normal component of the wall 
velocity is zero everywhere. 

Only a liquid at rest can satisfy the hydrodynamic equations for potential 
motion in such a container, since v = 0 satisfies both the equations and the 
boundary conditions. In a nonspherical container, however, the walls have a 
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normal vclocity component and the liquid is entrained by the walls; moreover, 
the kinetic energy for a given angular velocity of the shell becomes greater as the 
container becomes less spherical. 

We shall use as a model a nucleus in the shape of an ellipsoid. The semi-axes 
of the ellipsoid are denoted by c (symmetry axis) anda. It can be shown that 
the kinetic energy associated with potential rotation of a liquid ellipsoid is 


M (02 — a?) 
10 c*2 + a? 


9 
Wo", 


where wo is the angular velocity of rotation of the ellipsoid. Comparing this 
expression with the general expression for the kinetic energy of a rotating body, 


1/2 Jo, 
we have 
Ja AE, 
5 c?+ a? 


We note that the moment of inertia of a solid ellipsoid is 
M 
Jo = 5 (c? + a?), 


and that the ratio of these quantities is 


J ($ A 
Jo u c? + a? l 
This ratio approaches zero as (c/a) — 1 (sphere), and unity when (c/a) > œ. 


The volume of the ellipsoid is 


V = rar. 


The mean radius of the nucleus, which is the quantity generally considered in 
nuclear physics, is obviously 
R = (a%)'/s, 


We now consider the rotational energy levels. In a nucleus without spin 
(evcn-even) the energy levels are given by the expressions: 


I(T + 1) 312 10%? 
E eg, Ey = 0; Ea = i ze 
21%? 

Es = ——, etc. 
6 I etc 


We compare these expressions with the experimental energy levels for heavy 
nuclei. We must first verify the interval rule—the ratio of the energy differences 
of the levels—which should be 
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(E: — Ey): (E, — E): (Es — Es): ... = 3:7:11:..., 
or, computing the energy from the ground state of the nucleus, 


Es: Ei: Es: Es: ... = a A 2: eit aus 


This is in excellent agreement with the levels of two groups of even-even 
nuclei—in the rare earth region and in the region of heavy a-active nuclei. It is 
extremely interesting that there are no rotational levels in Pb?®, which is a 
doubly-magic spherical nucleus; it is only as one goes away from lead toward the 
heavier nuclei that rotational levels arise and the ratio Es: Ez approaches the 
theoretical value 3.3. The distance to the first excited level is approximately 
80 kev in the rare earth region and approximately 40 kev in a-active nuclei. 
It is thus possible to determine the moments of inertia. In the rare earths 
J = 1.9:107* g-cm?, while in the heavy a-active nuclei J = 3.8:107% g - cm?. 

The interval rule for the levels is best illustrated by examples. One of these is 
the U?*4 spectrum, which is measured from the energy of the a-particles of Pu? 
decay. Here there are three rotational levels with energies of 43, 143 and 297 
kev with spins and parity, respectively, 2*, 4* and 6*. The ratio of the energy 
levels is 1:3.3:6.9, in strict agreement with the theoretical values. 

A second example—the spectrum of the rare earth element 7,Hf!™ (Fig. 6)— 
appears in the decay of the isomeric state (5.5 hours). Here, four levels are 
apparent—2*, 4+, 6+ and 8+; these have energies of 93, 309, 642 and 1085 kev, 
which are in the ratio 1:3.3:6.9:11.7. The spins and parities of the levels are 
such that all transitions between these levels can be identified as electric quad- 
rupole transitions. 

The expression for distances between levels is different in odd nuclei. Since 
the spin in these nuclei is not zero in the ground state, the levels are determined 
by the expression 


22 
= — IE 1), 
E op ii +1) 


where I = lo, lo + 1, Io + 2, and so on (Jo is the spin of the ground state of the 
nucleus); an example of such a spectrum is given in Fig. 7. 

It is necessary to have one other experimental quantity to determine the 
elongation of the ellipsoid, since the shape of the ellipsoid 1s specified by two 
quantities, a and c. This can be done by introducing other data on nonspherical 
nuclei. A comparison of this type is carried out below. 

We now consider the individual nucleon states in a nonspherical nucleus. 

First of all, a good quantum number for a nucleon in an axially symmetric field 
is A, the projection of the nucleon orbital angular momentum on the ellipsoid 
axis. Since the spin of the nucleon can have a projection equal to ="/, along 
the nuclear axis, the total projection of the nucleon angular momentum w can be 
either A + %/2 or A — 1/2. We may note a single obvious exception. IfA = 0, 
w can assume only one value: w = !/2. 

In the ellipsoidal field the nucleon levels are designated in spectroscopy by the 
value of the projection A. Just as the values l = 0, 1, 2, 3,. . . are denoted by the 
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symbols s, p, d,..., the levels with A = 0, 1, 2, 3,... are denoted by the symbols 
TT, RER 

Inasmuch as a nucleon in a nonspherical nucleus does not have a definite 
angular momentum vector, we cannot use vector addition for the momenta. 
In a system of several nucleons we can only add the projections of the momenta 
on the nuclear axis. In view of this situation, the scheme for filling levels in a 
nonspherical nucleus proves to be completely different from that which applies 
for the central field in light nuclei. 

The chief difference is in the following. Particle levels with momentum j 
in a central field are highly degenerate. Each value of j corresponds to 2(27 + 
1) states with the same energy, differing either in charge or projection on an 
arbitrary axis. Thus, ina state with given j there can be 27 + 1 neutrons and as 
many protons. 

When a particle is in an axially symmetric field the degeneracy is considerably 
smaller. Ina level characterized by a given projection we can have two nucleons 
with +w, respectively. In this case, because protons and neutrons occupy 
different levels in heavy nuclei, the charge degeneracy (isotopic spin) does not 
exist. Hence, each pair of identical nucleons forms a closed shell which has 
neither spin nor magnetic moment. On the other hand, in light nuclei, we need 
2(27 + 1) nucleons to form a closed shell. Hence, to a considerable degree the 
properties of the heavy nuclei are determined by the properties of the last odd 
particle, a fact which leads to a somewhat paradoxical situation: heavy nuclei, 
with a large number of particles, can be described in terms of the single-particle 
model, whereas, as we have seen, this model is unsuitable for light nuclei. 

If there are two different particles beyond the closed shells—one proton and 
one neutron—these nucleons occupy different levels, and the behavior of the 
nucleus becomes more complicated. If the interaction between the nuclecns is 
small the energy of the level depends to some extent on the sign of the projection 
ot the angular momentum of the nucleus. Thus, one would expect the appear- 
ance of two levels, lying close together, with spins equal, respectively, to the 
sum of the projections of the moments of the two nucleons and the difference. 
Such a spectrum has been observed, for example, in the nucleus Am? * (Fig. 8); 
the lower level apparently has spin 5, while the spin of the upper level is zero. 
This picture corresponds to two nucleons, each with a spin projection of 5/2. 
The lower level corresponds to the sum of the projections, while the upper level 
corresponds to the difference. The spacing between the levels is approximately 
40 kev; this value gives an idea of the magnitude of the interaction between the 
nucleons. 

We start our analysis of the properties of heavy nuclei by considering spin. 
Disregarding odd-odd nuclei for the moment, we may say that the spin of the 
nucleus is determined by the state of the odd nucleon. The nuclear spin coin- 
cides with w, the projection of the total angular momentum of the nucleon on the 
nuclear axis. The projection of the orbital moment X is either w — 1/2 or w + 1/2. 


* Baranov and Shlyagin, Report to the Conference on the Peaceful Uses of Atomic 
Energy, Academy of Sciences USSR, July 1955. 
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Since each pair of identical nucleons comprises a closed system, in the region 
of heavy nuclei the spins of neighboring nuclei follow each other in an irregular 
way. 

The magnetic moment of a heavy nucleus is made up of two parts. In addi- 
tion to the magnetic moment due to the individual nucleons in the “quiescent” 
nucleus, there is an additional magnetic moment due to the rotation of the 
nucleus asa whole. The ratio of the magnetic moment to the mechanical mo- 
ment due to the rotation of the nucleus as a whole is given by Ze/2A mc, where Ze 
is the charge and Am the mass of the nucleus. In nuclear magnetons the gyro- 
magnetic ratio is given by 


a 
SE A 
Thus, the total moment u of the nucleus is the sum: 
Z 
= — K, 
u = um + F 


where n, as before, is a unit vector along the nuclear axis, K is the angular mo- 
mentum of the nucleus and po is the projection of the magnetic moment of the 
odd nucleon on the symmetry axis. 

We now determine the projection of the magnetic moment in the spin direction 
u. As we know, it is this quantity which is called the magnetic moment of the 
nucleus. By definition: 


whence 


MO Tad 
To compute the scalar product w-I, we multiply the expression for u by the spin 
of the nucleus I: 


Z 
yl = pon T+ FKL 


The nuclear spin I consists of two terms: I = Qn + K. Multiplying the last 
expression by I and recalling that Q = J in the ground state, we have: 
II+V)=R+Kl 


Here we have made use of the fact thatn-I = J. Whence K:l = /. 
Thus, 


7 
wel = pol + gE 


Finally, the expression for the magnetic moment of the nucleus is 
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The projection uo has an especially simple form in the case in which the odd 
particle is a neutron. Since the spin of the nucleon has a projection of either 
+1/, or —!/, along the nuclear axis, yo = +1.91. 

Assuming that Z/A = 0.45 for all nuclei (2/4 changes from 0.5 in light 
nuclei to 0.4 in heavy nuclei), for nuclei with even Z and odd A we have the 
formulas 


I I 
or u = 2.4——_, 
I+1 | ae | 





p= —15 


depending on the direction of the neutron spin with respect to its orbital moment. 

The formulas become more complicated in nuclei which contain an even 
number of neutrons and an odd number of protons because of the contribution 
of the orbital moment of the proton. The magnetic moment of the odd nucleon 
is equal to the sum of the projections of the intrinsic magnetic moment of the 
proton up = 2.8 and of its orbital magnetic moment, which is numerically equal 
to the orbital momentum in nuclear magnetons. The sum of the projection of 
the orbital moment of the proton and the projection of the spin is equal to the 
spin of the nucleus, J. Hence, the projection of the orbital moment is either 
I+ 1/, or I — 1/9, depending on the sign of the spin projection along the nuclear 
axis. Whence we have: 


or 


mw = 2.3 + Í and pw = —2.3 + I. 


Finally, the total magnetic moment of the nucleus is 


a det 
AAA — 
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To make a comparison of the calculated values of the magnetic moment with 
the experimental value we must assume that nearly all nuclei are nonspherical. 
We must presume at the outset that the level spectrum of a nucleus is rotational, 
i. e., that the nucleus is nonspherical. As has already been indicated, rotational 
level structures have been established only in a-active heavy nuclei and in nuclei 
in the rare-earth region. Hence, strictly speaking, a comparison of the formulas 
which have been derived can be made only for these nuclei. In the region of 
a-active nuclei, the magnetic moments have been measured, although very 
roughly, for several nuclei, and the known values of the magnetic moments of 
these nuclei are given in Table 5. 


Structure of the Nucleus 63 


The nuclei U? and Pu**! have one missing neutron. For spin 7/2 the negative 
magnetic moment is —1.1, which is close to the experimental value. With spin 
5/2, the magnetic moment can be either — 1.1 or 1.7, also in fairly good agreement 
with experiment. With spin !/2, the Pu?*% nucleus should have a magnetic mo- 
ment —0.5. Finally, in the Np?” case, which has one missing proton, with spin 
5/2 the magnetic moment should be 3.1, a value which lies within the limits of 
experimental error. 


TABLE 5 
Magnetic moment 
Nucleus Spin 

Experimental Theoretical 
920235 1/2(5/2) —0.8 0.2 —1.1(1.7) 
g3 Np 237 5/2 6 +42.5 3.1 
94P 11239 1/9 i (—)0.4 +2 —0.5 
„Pu? 5/2 | =(1.4 6) —1.1 (1.7) 


The agreement is somewhat poorer in the rare-earth region. However, in this 
region very little is known about the nuclear spectra or rotational levels. In 
hafnium, the rotational nature of the spectrum is well established, and the 
magnetic moment (only the absolute value is known) is 0.6 for spin Ya or 3/2, 
and is close to the theoretical value 0.5 or 0.9. 

It is apparent that further studies of the spectra and magnetic moments of 
the heavy nuclei will be required before more detailed conclusions can be reached 
as to the applicability of a model based on a nonspherical nucleus. 

An instructive example is furnished by the nucleus B1%%, which is comprised 
of the doubly-magic nucleus Pb% plus one proton. If we assume that this 
nucleus is spherical, its magnetic moment should be determined by the last 
proton. With spin 9/2 the orbital moment of the proton should be either 4 or 5 
(the go, state or the As,, state). In this case the value of the magnetic moment 
would be 6.8 or 2.6, respectively; both values are in disagreement with the 
experimental value, 4.08. This situation indicates that the addition of even one 
nucleon to a magic nucleus is sufficient to disturb the stability of the surface. 

Interesting results as to the shapes of heavy nuclei can also be obtained from a 
consideration of quadrupole moments. In the region with A > 50, the majority 
of nuclei have positive quadrupole moments which are several times greater 
than the quadrupole moment for a single particle. 

We now compute the quadrupole moment of a nonspherical nucleus. The 
quadrupole moment of a uniformly charged system with charge Z, which has an 
axis of symmetry, is given by the expression 


Q = 2 (322 — x? — y? — z?)dy, 
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where is the volume of the ellipsoid. Because of the axial symmetry 


Q = 2z: f (22 — x2)dv. 


In the case of a uniformly charged ellipsoid the following formula obtains: 


Q = (ct — a), 


However, the quadrupole moment computed from this expression is not 
equal to the quadrupole moment of the nucleus as measured experimentally. 
The quadrupole moment Qe is the O3; component of the quadrupole tensor in a 
coordinate system which is rigidly connected to a nucleus (in which the “3” axis 
is parallel to the symmetry axis of the nucleus). 

Experimentally, however, one measures the average value of the component 
of the quadrupole tensor (Oz. = O) of a rotating nucleus in a coordinate system 
which is fixed in space and usually defined by the direction of the external field. 
The relation between O and Qoe is given by the expression 


I 2I-—1 
= Q — 
I+12/-+3 
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This expression vanishes when J = 0 and J = 1/2. 

The factor I/(I + 1):(22 — 1)/(2I + 3) strongly reduces the quadrupole 
moment and approaches unity very slowly. When J = 1 this factor is 1/1) and 
when J = 3/, this factor is 1/5; even when J = 9/» this factor is still far from 
unity, being $/ı. 

Recalling the expression for the moment of inertia 


ee ee 
O 5 +g?’ 


we can find the semi-axis of the ellipse for every nucleus in which the quadrupole 
moment and rotational levels are known. 

The quadrupole moments Q of the odd heavy nuclei are known from spectro- 
scopic observations. However, these measurements determine only the inter- 
action energy of the nuclear quadrupole moment with the electric field due to the 
electron shell, and this calculation is not very reliable. There is another more 
direct way of determining the magnitudes of quadrupole moments; this method 
is independent of the atomic properties. If it is established that the levels 
exhibit rotational structure (by determining the intervals between levels), it 
can be shown that the lifetime for each level* is proportional to the square of the 
true quadrupole moment of the nucleus Qo (not the spectroscopic Q). This 
method also makes it possible to determine the quadrupole moments in even- 


* Also the probability for excitation by an electric field (in collisions with charged 
particles). 
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even nuclei, in which there is no spin and in which O = 0.* If now, from the 
values of Qo and the moment of inertia, J, we determine the semi-axes of the 
ellipsoid (the ratio Q?/J yields the sum of the squares and Q the difference) it is 
found that the ellipsoid is very highly elongated. However, in this analysis it is 
found that the volume of the nucleus is extremely small; the volume which is 
obtained corresponds to a value of 0.7-10713 cm for the constant ro which 
appears in the expression for nuclear radius R = nA”. On the other hand, 
data from electron scattering experiments give a considerably higher value for ro: 
approximately 1.1-10713 cm. The data obtained from the spectra of mesic 
atoms are in agreement with this latter value. Moreover, if one computes the 
semi-axes of the ellipse from the radius of the nucleus and Qo, the moment of 
inertia calculated from the potential flow model, turns out to be considerably 
smaller (three times) than the value determined from the observed spacings 
between rotational levels. 

To some extent these discrepancies may be a result of the fact that the radius 
of the charge distribution in the nucleus (which determines Q) is somewhat 
smaller than the radius of the mass distribution (which is associated with the 
moment of inertia). However, the basic cause would seem to lie in the approxi- 
mate character of the hydrodynamic model. This problem will require a great 
deal of further investigation. 

Although there is no quantitative agreement, the experimental facts indicate 
that the heavy nuclei far from the magic nuclei are elongated ellipsoids of rota- 
tion, and that the elongation (the ratio c/a) is at least 1.5. Hence, it cannot be 
assumed that the field of a heavy nucleus is spherically symmetric. 


APPENDIX 


In this appendix we derive several formulas which have been used in the 
present lecture. 

A. The moment of inertia for potential motion of a liquid ın an axtally-sym- 
metric ellipsoidal container. We choose the z-axis as the symmetry axis of the 
ellipsoid and the y-axis as the axis of rotation. 

Potential motion of a liquid is described by a potential y which satisfies 
Laplace's equation 

Ap = 0. 


The velocity of the liquid v is given by 
y = grad p. 


The boundary conditions are that the normal component of the surface velocity 
of the liquid must be equal to the velocity of the container. 


* The quadrupole moment has been measured by both methods in the odd nucleus 
Tal8l, In this case the spectroscopic value O = 7 (or Qo = 15, since the spin of Ta!!! 
is 7/2) has been found to be approximately twice as great as Qo measured from the 
lifetime. Comparing these values with Qo for neighboring nuclei, it may be assumed 
that the lower value of Qo is the correct one. 
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We write the equation for the ellipse in normal form. 


x+y gh 
a? ch 
The velocity of the surface is © X r. 


The direction cosines of the normals to the surface of the ellipse are propor- 
tional to x/a?, y/a?, and 2/c?, respectively. 
Multiplying the components of the difference 


Veur — QU x r= (grad g )sur — Q) x r 


by the direction cosines of the normals and equating the expression which 1s 
obtained to zero (this is obviously the condition that the normal components 
vanish), we have: 


x [de y 0% z fee ) 
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This relation must be satisfied over the entire surface. We can satisfy this rela- 
tion by choosing as a solution of the Laplace equation the function 


g = Áxz, 


where A is a constant. Substituting this function in the boundary conditions 
we have: 


1 1 
(Au) +Z (4 +e) = 0. 
a C 


From this, having determined A, we find the potential 


c? — a? 


w x 
c+ a? 





g = 2. 


We now compute the kinetic energy of the liquid. Setting the density equal 
to unity, we write: 


r= f viv == f dime? ——*) f e 2d 
=> pan grad p)dv = 7 a. Z x?)dv. 


Making the usual substitution of variables x = af, y = an and z = ct, after some 
simple transformations, we have 


wt [c2 — a?N? 
Ts ə (@ 7 =) a’c(c? + a?) f E2dtdndt, 








where the integral is taken over a sphere of radius unity. Finally: 


a” M (c? = a2)? 


T = : 
25 +a 
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where M is the mass of the ellipsoid (equal numerically to its volume). It is 
apparent that the coefficient of w?/2 is the moment of inertia in question: 

7 M (2 — a2)? 

5.0 + a? 


J 


B. Quadrupole moment of a uniformly charged rotating ellipsoid. To deter- 
mine the factor by which we must multiply Qo, the quadrupole moment of the 
fixed nucleus, to find the quadrupole moment of the rotating nucleus, Q, we form 
the quadrupole tensor from quantities which characterize the fixed nucleus and 
then from quantities which characterize the rotating nucleus. 

We use the symbol n; to denote components of the unit vector directed along 
the ‘‘3’’ axis of the fixed nucleus. The quadrupole tensor can be written 


Qin = 3/2Qo(nine — */30:1). 


The coefficient 3/2 is chosen to make Qo = Os. 
For the rotating nucleus Qi, is expressed in terms of the spin vector I. As in 
Lecture Four, we can write: 


Qo 


al m1) {IIs + Il: — 2/841(1 + 1)). 


On = 


Setting this expression equal to the mean value of the quadrupole moment com- 
puted by the preceding formula, we have 


Q 


Tar way Wile + Ile — Yu + 1)}. 


Ol ninn — */38:x) = 


We multiply both parts twice by the spin vector and make use of the relation 
Iımı = Q, where Q is the projection of the spin on the symmetry axis of the 
nucleus (equal to the spin J), and the relations 

Y Lidl = 14 DP, doubt, = UI + 1D? U + 1); 

i, k i, 
the latter relation follows from the commutation rules for the components of the 
spin vector 

Laly — Iyle = 212, etc. 
As a result we find: 
Q 


1 4 A pa _ 
ode — 2 10+ vt gruen II + D? 


1 
3 QU + 1X21 + 8), 


whence 
302 — (1 + 1) 


T (I+ DOI + 3) 
With 2 = I we find the formula in question: 


_ I(2I — 1) 
(14 DQOI+3) 


Q Qo. 


Q Qo. 
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Nuclear Reactions 


(STATISTICAL THEORY) 


T.. concept of individual nucleon levels is meaningful only as 
long as the interaction between nucleons is smaller than the distance between 
levels (with the same moment and parity) corresponding to different nucleon 
configurations. When the excitation energy increases, the density of levels 
grows sharply and the division of nucleons into shells loses its significance. 
This situation arises when we deal with heavy nuclei and relatively high excita- 
tions. 

In these cases we must depart from our description based on individual- 
particle mechanics and go over to the other limiting case, that is to say, we must 
use an analysis based on statistical mechanics and thermodynamics. In doing 
so, however, we must keep in mind the fact that since statistical laws are valid 
only in systems characterized by a large number of degrees of freedom, the 
application of these laws to nuclei in which the number of particles is not very 
large obviously cannot yield accurate results. 

From the point of view of statistical mechanics an excited nucleus is character- 
ized, first of all, by an energy E, which may conveniently be measured from the 
ground state of the nucleus. On this scale E is simply the excitation energy. 
Since there are many levels in a system with a large number of degrees of freedom, 
the distribution of levels may be described conveniently by the mean spacing 
between levels at a given excitation energy D(£); similarly, this distribution can 
be given in terms of the density of levels at a given point in the spectrum w(E) = 
1/D(E). The level density w(E) changes with the excitation energy. As is 

rell known, the level density increases with excitation energy. 

The problem of a statistical analysis, then, is to establish a relation between 
the level density and the excitation energy. The relation between these two 
quantities is established by introducing the idea of entropy. 

In statistical mechanics one establishes a relation between the entropy S 
and N, the number of states having approximately the same energy in a given 
system. This relation is the well-known Boltzmann formula 


S(E) = In N 


(in our analysis the entropy S is measured in dimensionless units). 
The inverse relation then follows: the number of levels lying in a certain 
energy interval AZ is given by the expression 
69 
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N = exp [S(E)]. 


The energy interval AE is determined by the magnitude of the energy fluctu- 
ations in the thermodynamic system at constant temperature. This quantity 
determines the inherent uncertainty in the energy of the system and can be 
computed for any thermodynamic system. 

From the last expression we can also determine the mean spacing between 
levels D(E). 

Assuming AE = ND, we have: 


D(E) = AE exp [—S(E)] = exp [-S(E) + In AE]. 


Having determined the entropy of the system, without any further assump- 
tions we can introduce the temperature of the nucleus, using the well-known 


thermodynamic expression 
-1 
sa (28)7. 
DE 


To apply statistical considerations to the nucleus and to compute the entropy 
of a system of nucleons we must, first of all, know the way in which the excitation 
energy of the nucleus depends on temperature or, what is the same thing, the 
heat capacity of the nucleus. 

The system of nucleons which comprise the nucleus may be considered in 
terms of an ideal gas of particles which obey Fermi statistics. 

At first glance this assumption seems strange since there is a strong interaction 
between nucleons. However, it is well known that a system of interacting Fermi 
particles has an energy spectrum very similar to the spectrum of an ideal de- 
generate Fermi gas. A convincing example is the spectrum of electrons in a 
metal. It is well known that the heat capacity of metals at low temperatures 
(when the contribution of the lattice to the heat capacity is negligibly small) is 
proportional to the temperature, just as in an ideal gas. At the same time, it is 
obvious that the interaction of electrons in a metal is not small. Hence, there is 
a basis for assuming that the heat capacity obeys a similar relation in the nuclear 
case. 

We shall assume that the heat capacity of the nucleus depends on temperature 
in the same way as does the heat capacity of a degenerate Fermi gas* 


a 
6S E 
2 


where a is a constant. Since 


dE 
dT’ 

* The condition for “strong degeneracy” of a Fermi gas is given by the inequality 
T< (A2/M)p”°, where p is the density of particles. Ina nucleus p = (%/3rr03)71 
(ro = 1.2-10713 cm), whence the degeneracy condition assumes the form T & 10 
Mev; it may be assumed that the temperature of the nucleus is low enough since 


T ~ 1 Mev corresponds to an excitation energy of approximately 10 Mev in a heavy 
nucleus. 
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we can determine the energy: 


„DS _ 
ar ” 
which yields 
S=5T = (aE), 


The constant a should be proportional to the number of particles in the nucleus 
because at a given temperature the energy of the nucleus must obviously be 
proportional to the amount of matter. Hence, 


a KA and E = Z AT, 


where the coefficient K is relatively independent of the properties of the nucleus. 

Knowing the entropy of the nucleus, we can also compute the mean spacing 
between levels. However, for this purpose we must calculate the entropy in a 
more exact fashion than is possible using the formulas given above. In this 
calculation logarithmic terms appear in the entropy; in the e* expression these 
terms appear in the factor which multiplies the exponential. Thus, if we use the 


expression S = VaE for the entropy it is impossible to calculate the factor in 
front of the exponential by ordinary fluctuation theory. 

Since a rigorous calculation of the entropy of the nucleus has not been carried 
out up to this time, the quantity A£ in the formula 


D = AE exp [(aE)-"?] 


cannot be computed from theory. On the other hand, the present experimental 
data on excited levels in heavy nuclei are inadequate for carrying out a detailed 
comparison of these formulas with experiment. 

In discussing the experimental data on nuclear levels, we must not forget the 
manner in which these data are obtained. A significant number of levels are 
determined from experiments on slow-neutron capture. In reactions of this 
kind we cannot excite arbitrary levels since a slow neutron can penetrate the 
nucleus only when its orbital moment is zero. For example, if the initial 
nucleus has spin J = 0, the intermediate nucleus can have only one spin value, 
namely 1/2. Hence, out of the large number of levels which are available in a 
heavy nucleus, the slow-neutron capture reaction selects only levels with lowest 
spin. The mean level density determined by experiments of this kind is smaller 
than the true density. In other experiments, for example, in experiments on 
inelastic scattering of fast neutrons, there are virtually no limitations due to spin; 
the density found in these cases is closer to the correct density. 
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In this connection it is interesting to examine the level density associated with 
a given spin value. 

In accordance with the laws of statistical mechanics the probability that a 
system (in the present case, a nucleus) will be found in a state with a moment 
lying in the interval dM,dM,dM, about the value M is proportional to 


dW( Mr, My, Ma) e VU TIM? + My + Me) Y, dM, dM.. 


The proportionality factor is determined by the normalization. Thus, for the 
density of levels we have 


doM) = w(2QrIpT)~/* e7 “PIDA MAM, 


where w is the total density of levels. In quantum mechanics the moment is 
quantized and this formula becomes 


(21 + 1)2 (=)", 1 
vr \Qh)/ T’ 


where w,(E) is the density of levels with a given value of moment J. 

In discussing nuclear levels we have assumed above that these levels are inde- 
pendent. The formulas given determine only the mean distance between 
levels; the distribution function—that is, the probability of a given spacing 
between levels for a given mean spacing—has not as yet been investigated. 

There is no basis for assuming that nuclear levels with different spins are in 
any way related; the spacing between levels with different spins should be 
distributed according to a random law. 

However, this conclusion does not apply to levels with the same spin. The 
problem of the distribution for these levels has not been solved. At this point 
it is reasonable to assume that, in any case, levels with the same spin are not 
distributed according to a random law. From very general considerations one 
expects that levels with the same spin are distributed in such a way that the 
probability of a very small spacing will be very small—the levels “repel” each 
other. This problem requires further investigation. 

Statistical considerations prove to be extremely useful for studying both 
nuclear reactions and nuclear decay. We start with the simplest question—the 
probability of a-decay. Until recently contradictory theoretical expressions 
have appeared in the literature. The probability for alpha decay, as is known, 
is determined basically by the probability that an a-particle will penetrate the 
potential barrier. In the theory of this effect the penetrability of the barrier 
(the ratio of the square of the modulus of the wave function on both sides of the 
barrier) is given by the expression 


e- (R/U HN, 





w(E) = w(E) 


Poet, 


where T is a quantity which depends on the energy of the a-particle and the 
charge of the nucleus. The decay constant A (the reciprocal of the probability) 
should be 
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where v and A have the dimensions of frequency. It is the quantity v which is 
the source of misunderstanding. In the first papers by Gamow the nucleus was 
described as a potential box and the quantity v was interpreted as the frequency 
with which the «-particle kept “banging” against the edge of the box. Whence 
the following value was obtained: 


v 
‚Rp 


where R is a typical dimension of a nucleus and v is the “velocity” of the a- 
particle in the nucleus. This frequency corresponds to an energy hv ~ hv/R, 
which is of the order of tens of millions of electron volts. It is clear that this 
description is not an accurate one. 

Later, another point of view was suggested; this model led to a factor which 
corresponded to an energy which was too small—of the order of several electron 
volts. 

A statistical analysis of the nucleus makes it possible to obtain a rigorous 
formula for the factor in question. This formula is of the form: 


D -r 

dd 2rh ° 
The coefficient D/2rA has a simple physical meaning: it is the frequency of a 
linear oscillator which has the same spacing between levels (equal to hv) as the 
nucleus in the excitation interval being considered. This quantity is of the 
order of several tens of kiloelectron volts. 

In this form the formula relates the probability of a-decay to the total level 
density for the nucleus and not to the behavior of the «-particle inside the nu- 
cleus, as is the case in the Gamow formula. Hence, in deriving the formula there 
is no need for discussing the existence of the «-particle inside the nucleus. 

The thermodynamic approach proves to be especially useful in analyzing a 
second type of nuclear reaction—the emission of particles from an intermediate 
nucleus. An excited intermediate nucleus can be considered a heated body 
from which particles are evaporated. The evaporation probability can be 
calculated in exactly the same way as the evaporation rate for ordinary liquids. 

In this picture we consider the nucleus to be in thermodynamic equilibrium 
with the gas of evaporated particles. In a system of this kind the number of 
evaporated neutrons is equal to the number of absorbed neutrons; moreover, the 
number of absorbed neutrons can be calculated easily. Since wı are not inter- 
ested in the absolute probability, but only in the relative probability of emission 
of neutrons with different energies (neutron evaporation spectrum), we neglect 
all factors which do not depend on energy. The energy distribution of the 
neutrons in the gas which surrounds the nucleus is given by 


fis de ~ e E DEA. 


The number of neutrons with energy & which collide with the surface of the nu- 
cleus in unit time is proportional to the numbcr of neutrons f(8) and the neutron 
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velocity (1. e., 8. Thus, the number of neutrons which strike the nucleus 
and have energy 8 is given by 


v(s)de ~ e EDede. 


However, not every collision of a neutron with a nucleus results in the formation 
of an intermediate nucleus. We shall denote the probability for an event of this 
kind by &(&). The coefficient E is called the sticking coeficient (in molecular 
physics the analogous quantity is called the accommodation coefficient). Col- 
lecting all factors and using the fact that the number of absorbed neutrons is 
equal to the number of emitted neutrons, we can write the distribution relation 
for evaporated neutrons in the form 


n(g)d& ~ (e)se 7 &/ Dde. 


At energies of several million electron volts the sticking coefficient approaches 
unity. If we assume ¿ = 1, the mean energy of the evaporated neutrons be- 
comes: 


Smean = 2T. 


At low energies the coefficient & is considerably different from unity; at 
energies close to zero (thermal region), the sticking coefficient is proportional to 
velocity: 


¿~ elt (8>0) 


whence 
Emean TE; T (8 0). 
Thus, the mean energy of the evaporated neutrons satisfies the relation 
2T < mean < 2.91. 


It is important that the energy of the evaporated particles be considerably 
smaller than the excitation energy. Under these conditions the excitation 
energy is dissipated in the emission of several nucleons with comparatively smali 
energy rather than in the emission of a single fast nucleon. 

The evaporation of neutrons has one distinguishing feature by which it differs 
from evaporation in ordinary liquids. The energy of a molecule which is 
evaporated from the surface of a liquid is small compared with the energy of the 
total liquid, so that the temperature of the liquid remains constant (isothermal 
evaporation). In the case of nuclear evaporation, however, the neutron carries 
off (considering its binding energy) a considerable fraction of the available 
energy. Hence, to obtain really accurate formulas we must indicate the proper 
nuclear temperature. It is easy to show (inasmuch as the derivation of the 
formula was based on an analysis of the inverse problem—the equilibrium of 
the final nucleus with the neutron gas) that the quantity T which appears in 
the formulas is the temperature of the final nucleus. 

The evaporation of charged particles—protons, deuterons and heavier nu- 
clides—can be considered by an analysis similar to that used in the evaporation 
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of neutrons. Proton evaporation differs from neutron evaporation because the 
Coulomb repulsion of the protons means that the absorption probability ¿(8) is 
smaller at lower energies and that the limiting value ¿ = 1 is reached at higher 
energies. 

From this point of view the case of the deuteron is of special interest. It is 
known that there is a comparatively large number of deuterons in cosmic rays. 
These are undoubtedly the result of nuclear reactions. It seems strange, at 
first glance, that a nucleus can emit a comparatively unstable particle such as 
the deuteron in one piece. 

A consideration of deuteron evaporation from the nucleus shows that the 
probability for evaporation is determined only by the excitation energy (tem- 
perature) of the final nucleus and is, consequently, virtually independent of the 
binding energy of the deuteron. Even if the binding energy of the deuteron 
were zero its evaporation probability would not vanish. These same consider- 
ations would seem to explain the fact that the Be® nucleus is emitted from nuclei 
in nuclear reactions; this nuclide is extremely unstable and decays into two a- 
particles in a time of approximately 10714 sec (decay energy approximately 
50 kev). 


LECTURE EIGHT 


Nuclear Reactions 


(OPTICAL MODEL. DEUTERON REACTIONS) 


I n this lecture we shall consider nuclear reactions from an 
entirely different point of view. To eliminate the complications which arise 
because of the Coulomb field of the nucleus we shall consider neutron reactions 
only. 

Consider a neutron which collides with a heavy nucleus. There is consider- 
able probability that the neutron will pass through the nucleus without losing 
energy, suffering only a change in the direction of its motion. This process 
corresponds to elastic scattering. However, a neutron can also be absorbed, 
thereby producing an intermediate nucleus which can then decay by various 
modes. In particular, the intermediate nucleus can emit a neutron with smaller 
energy (inelastic scattering) or with energy equal to that of the initial neutron 
(re-emission of a neutron). 

If we are not concerned with the intermediate nucleus, all events which occur 
in the collision of a neutron with a nucleus can be divided into two types— 
scattering (without the formation of an intermediate nucleus) and absorption 
(formation of an intermediate nucleus). 

Using this picture we can consider the collision of a neutron with a nucleus in 
terms of the passage of the neutron through a certain medium—nuclear matter— 
and consider this entire problem by analogy with the passage of light through a 
medium which exhibits refraction and absorption. 

Just as in the case of an optical medium, we can introduce an index of refrac- 
tion, which is a complex quantity if there is absorption; thus, we require two 
quantities to characterize the “optical” properties of nuclear matter (in the 
following we will speak of the optical properties of the nucleus). This approach 
to nuclear reactions is called the optical model of the nucleus. 

In order to determine these parameters we turn to the Schrödinger equation 
for the neutron. According to quantum mechanics the motion of a nucleon is 
described by a wave function which is associated with a certain potential. If 
the potential is real we are dealing with neutron scattering; on the other hand, 
as we shall now see, neutron absorption is to be associated with a complex 
potential. 

To verify this statement we introduce the expression for a flux of particles, 
using a wave equation with complex potential. We will proceed exactly as in 
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deriving the expression for a flux of particles using an ordinary potential. 
We write the wave equation for a particle in a field described by a complex 
potential U — iV (U and V are real functions): 


hz 
Ey = — — Ay + Uy — iVy, 
2m 
where E is the energy of the particle. We now write the conjugate equation: 
h2 
Ey* = — — Ay* + Uy* + 1Vy*. 
2m 


Multiplying the first equation by y* and the second by y, we take the difference 
of these two expressions. After some elementary manipulation we have 


22 
— om Lv (y* grad y — y grad y*) — 2Vyy* = 0. 
m 
Since yy* is the density of particles p, and #/2mi(y* grad y — y grad y*) is the 
current density j, the last relation can be written: 


divj = —2pV. 


The left-hand side of this expression is the divergence of a current; from the 
physical meaning of this quantity, we see that it vanishes if particles are neither 
produced nor absorbed. Hence, the right-hand part of the equation is nothing 
else than the source intensity. Thus, the absorption of particles per unit time 
per unit volume is determined by the quantity 2Vp = 2Vy*y. 

The imaginary part of the complex potential is then the absorption coefficient. 

In order to formulate a theory for the passage of a neutron through the nucleus 
starting from these ideas, we must have some information for determining both 
parts of the potential (real and imaginary). There is no reason to assume that 
these quantities are independent of energy; on the contrary, they are un- 
doubtedly functions of the neutron energy, and must be determined from 
experiment. 

To make a clear distinction between processes connected with the formation 
of an intermediate nucleus and those connected with true elastic scattering we 
must first make sure that the probability for the re-emission of a neutron with 
the initial energy is small. In turn, this requirement implies that the inter- 
mediate nucleus must be capable of decay by many modes. This condition is 
satisfied if the neutron has an energy of the order of several million electron volts; 
under these conditions there are a number of possible reaction channels and 
excited states available for the final nucleus. 

Transforming the Schrédinger equation as written above, we can show that in 
the optical model the neutron motion is given by the equation 


2m 
Ay + yr (B -— U + ¿VW = 0. 
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This equation is analogous to the wave equation in vacuum 


Ay + ky = 0, 


\ 
where k is the propagation vector. 


Formally, we can say that the magnitude of the propagation vector, which is 
k = (2mE)"?/h 
in vacuum, becomes a complex quantity in the nucleus: 


2my/ 
K = Ki + ik, = © (E- U+ iv)”. 





By analogy with optics, we can speak of a complex index of refraction for nuclear 
matter: 


n = (Kı + Ka). 


=> | = 


This quantity, however, is usually not employed in nuclear physics. 

We now examine the physical meaning of the complex propagation vector. 
When the neutron is inside the nucleus the particular solution of the wave 
equation can be written in the same form as the vacuum solution, i. e., a plane 
wave (but with a complex propagation vector). 

Limiting ourselves to the one-dimensional problem, we can write 


Y = et Kk = et Kixy— Kox 


It is obvious from this expression that the particle absorption depends on Ka. 

In order to relate the complex propagation vector to the quantities which are 
determined from neutron scattering experiments, we must solve the problem of 
diffraction of a neutron wave on a nucleus with given optical properties and 
radius. 

We must determine three quantities from experiment at each value of the 
energy: the real and the imaginary parts of the potential energy, and the 
radius of the nucleus (independent of energy); thus, we require rather detailed 
experiments on both elastic and inelastic scattering to determine all quantities 1f 
additional assumptions are to be avoided. 

It should be noted that in general these constants can be different for protons 
and neutrons. This difference may be due, first, to Coulomb forces, and second, 
to the fact that the number of protons and neutrons is different in heavy nuclei. 
At the present time this very interesting question has not been investigated 
experimentally. 

In order to solve the problem of neutron diffraction on the nucleus we must 
find the solution of the wave equation inside and outside the nucleus. 

The inside and outside solutions must be matched at the boundary. The 
derivatives must also be matched at the boundary. Usually the solution is 
expanded in a series of spherical functions and the matching conditions yield 
algebraic equations for the coefficients of the series. The solution of this prob- 
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lem requires the calculation of a comparatively large number of terms and proves 
to be rather difficult. 

The problem can be simplified considerably if we consider high energy neu- 
trons; in this case the wavelength of the neutron is small compared with the 
dimensions of the nucleus. The neutron wavelength A, is given in centimeters 
by the formula 


where E is given in millions of electron volts. The radius of the nucleus R is of 
the order of 5-10-13 cm; hence, the condition An X R is satisfied at neutron 
energies of several hundreds of millions of electron volts. Using neutrons of this 
kind, we can solve the problem by a method which is very similar to that used 
in optics to solve Fraunhofer diffraction, for example, the diffraction of a 
parallel beam of light at a circular aperture. 

For simplicity we shall consider a spherical nucleus. If the nucleus is not 
spherical, the problem is considerably more difficult. Suppose that the neutron 
beam is incident along the x-axis. The incident neutron beam is described by 
the plane wave e***, where k is the magnitude of the neutron propagation vector 
in vacuum. The amplitude of the wave function is taken as unity. 

Consider a plane, perpendicular to the incident neutron beam and located 
directly behind the nucleus (i. e., on the side away from the incident beam). We 
will call this plane the “shadow” of the nucleus, i. e., that region which neutrons 
moving along rectilinear trajectories parallel to the propagation vector k can 
reach only by passing through the nucleus. Asan approximation we assume that 
at all points of the plane outside of the shadow region the wave function is a 
constant whose value is that which obtains in the absence of the nucleus. 

Since we can always multiply the wave function by an arbitrary phase factor 
we take the value of the phase to be zero at all points lying outside the shadow 
region. This means that we have chosen the z-coordinate in such a way that the 
plane corresponds toz = 0. Thus the wave function 


becomes unity outside the shadow region. 

To compute the wave function in the shadow region we note that inside the 
nucleus (over a path of length A, which is different for different parts of the 
shadow), the propagation vector is not k, but Kı + iK. If this segment of path 
were in vacuum the phase would be changed by an amount kA. Inside the 
nucleus, however, the actual phase change is (Kı + iK2)A. Thus, in passing 
through the nucleus the neutron wave function undergoes an additional phase 
shift, equal in magnitude to 


6 = (Ky + 1K? = k)A. 


We have assumed that the wave function is unity outside the shadow; in the 
shadow the wave function is then 
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Wshadow = er(Kı a k) Ap — KoA 

From a knowledge of the wave function we can also compute the cross sections. 
It is simplest to use the absorption cross section. The neutron density (equal to 
y*y) in various points of the plane is simply 


e—2K24 in the shadow region, 


1 at other points 


we =} 


If K: = 0, the density is unity at all points in the plane. 

It is apparent that the difference between unity and e~?*** corresponds to the 
absorption of a neutron by the nucleus. The attenuation of the beam which 
takes place as a result of absorption by the nucleus is thus given by 1 — e~?4*4. 
From this we easily obtain the cross section. Specifically, integrating the 
attenuation function over the cross section of the beam, i. e., over the plane, we 
have 


Öabs = f (1 = e™?K:A)dyxdy, 


where the integration is actually carried out over the shadow region; obviously 
A depends on the coordinates x and y. 

It is somewhat more difficult to obtain the formula for elastic scattering. In 
the absence of a nucleus the neutron wave function would be a plane wave and, 
in a coordinate system in which the z-axis is parallel to the wave vector, the wave 
function would be independent of the x and y coordinates (the wave function is 
simply e***). The presence of the nucleus changes the wave function in such a 
way that it becomes a function of the coordinates of the plane. Thus, the wave 
function is no longer a plane wave, 1. e., it no longer describes a parallel beam of 
particles. Inasmuch as the neutron energy is not affected (processes which lead 
to a change of energy are taken into account by the absorption), the new wave 
function describes particles which are scattered in different directions. In 
order to distinguish the different particles by direction, in accordance with the 
laws of quantum mechanics, we must expand the wave function in the xy plane 
in a Fourier integral. The terms in this expansion (more precisely, the squares 
of the moduli) then determine the particle flux in various directions. The 
Fourier expansion must be carried out in the two variables kz and ky, correspond- 
ing to the two coordinates of the xy plane. If there were no scattering this 
expansion would contain only one component, since kz = ky = 0. 

Thus we write the wave function in the form 


y = (27) f a(kz, k,)e Rex + hu dk,dky. 


Then, in accordance with the laws of quantum mechanics, la(ke, ky)|? is the 
quantity which determines the fraction of neutrons with propagation-vector 
projections on the x-axis and y-axis in the intervals kz, kz + dk, and ky, ky + 
dk,. From the Fourier theorem 


alka ky) = (2r) | yx, yet + dady, 
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The function y(x, y), which must be substituted in this formula as has been 
mentioned above, is exp [7(K, — k)A — K2A] in the shadow region, and unity 
at other points. The scattering cross section is given by the expression 





do(kz,ky) = (21)? f Ux,y) exp [—ilkex + kyy)ldxdy | ? dkzdky. 

To transform from the distribution over propagation-vector components 
(kz, Ry) to a distribution over scattering angle, we introduce the vector in the 
xy plane; the magnitude of this vector determines the scattering angle. Thus, 
kz = kôz and ky = ky. Since the scattering is independent of the azimuthal 
angle, the cross section depends only on the length of the vector 6. The integral 
which determines the scattering cross section can be computed, although the 
expression which is obtained is extremely complicated. 

To obtain the total scattering cross section we can integrate the differential 
cross section over angle or propagation-vector components kz and ky. However, 
there is a simpler method. 

We have seen that the wave function y differs from unity only over a small 
section of the plane. Thus we can write y in the form 


y =1+(Y — 1. 


Written this way, the unity term can be interpreted as the incident wave 
propagating along the axis (we may recall that we are speaking of the value of 
the wave function at the given plane). Then the difference y — 1 describes the 
change in the wave function due to diffraction at the nucleus. The above ex- 
pression can be interpreted as the sum of the incident and scattered waves. 
The total number of elastically scattered neutrons can be obtained if we integrate 
|y — 1|? over the plane. Thus we have 


Celastic = f ly = 112 dxdy. 


As in the case of absorption, the integration in this formula is actually carried 
out only in the shadow region. Since ly | < 1 for absorption, any absorption 
must be accompanied by elastic scattering. 

The total cross section for the interaction of neutrons with the nucleus is 
written as the sum 


Ototal = Celastic + Jabs. 


Combining the expressions for both cross sections, we have: 


Ototal = f (1 — |y]? T ly — 1|2)dxdy 
or, taking out the expressions for the moduli: 
1 — |yl?+ |y — 1|? = 2 — 2Rey, 


where Re denotes the real part. 
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Thus, 
Stotal = 2 f (1 — Rey)dıdy. 


We can also express ototaı in terms of the functions K, and Ke. Substituting the 
expression for y and integrating over the shadow region, we have: 


Ototal = 2 f (1 — e2K:A cos 2(Kı = K:)A)dxdy. 


Thus, we have obtained all the relations needed to describe the interaction of a 
neutron with a nucleus. At each neutron energy the values of the two cross 
sections, for example, ototaı and gaps, are related to the values of K, and Kat this 
energy by the above formulas. 

At the present time, the experimental data are not sufficient to determine all 
three optical characteristics of the nucleus: the complex index of refraction 
k-\(K, + iK) (which is a function of energy) and the radius of the nucleus 
(which is independent of energy). To determine these quantities, it is obviously 
not sufficient to know only the magnitudes of the total and elastic cross sections 
for any given nucleus; it is also necessary to know at least the elastic-scattering 
angular distribution. 

At this point we should add that the nuclear radii determined from neutron 
cross sections are considerably larger than those found in electron scattering and 
mesic atoms. However, one must remember that this difference does not 
necessarily imply a discrepancy, since we are actually speaking of effects in 
which the radius may not be defined in the same way. Moreover, the radius 
obtained from neutron scattering should be larger than the radius of the nucleus 
by an amount equal to the radius of the neutron itself (range of nuclear forces). 

Several conclusions can be reached from neutron scattering at high energies. 
If the nucleus were completely opaque to neutrons the absorption cross section 
would be rR?. Indeed, in this case y would become zero in the shadow region 
and the integral over 1 — |y|? would be equal to the area of the shadow. Ex- 
perimentally, the absorption cross section proves to be smaller than this value. 
This situation means that the nucleus 1s a semi-opaque medium for fast neutrons. 

We consider the case in which the nuclear refraction is large (this is the actual 
physical case), so that 


(Kı — k)A>1. 


In the shadow region the quantity Rey = e-R:A cos (Kı — k) A, which appears 
under the integral sign and determines the total cross section, changes sign 
many times. Since this function oscillates, in the integration over the shadow 
region the total integral is small. Hence, in cases of large refraction Rey be- 
comes effectively zero, and the formula for the total cross section assumes the 
form 


Ctotal = 27rR?. 


The ratio of cross sections for elastic scattering and absorption obviously de- 
pends on the magnitude of the absorption. 
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We have seen that at large absorption the absorption cross section approaches 
TR?; hence, in the limiting case 


Gabs = Gelastic = TR? for KeA >1. 


On the other hand, at small K: there is essentially no absorption, and the scatter- 
ing cross section becomes the same as the total cross section: 


Cabs = O, 
Gelastic = 2mR? (KA < 1). 


We have discussed only the simplest properties of the optical model. In an 
actual analysis of the experiments it must be considered in greater detail. First 
of all, the considerations given above must be extended to include the case in 
which the incoming particle is a charged particle. This situation leads to a 
certain complication 1n the formulas. 

The greatest difficulty appears if we give up the assumption that the nucleus 
ig a sphere with sharply defined boundaries. In an actual nucleus the density 
does not change sharply, but falls off more or less gradually at the boundary. 
Moreover, in heavy nuclei we must take account of the fact that the nucleus is 
not spherical. The introduction of these refinements makes the formulas 
complicated and makes an analysis of the experimental results very difficult. 

Interesting results are obtained if one takes account of the spin of the incoming 
particles. Because of spin-orbit coupling (which has already been discussed in 
our review of the light nuclei), particles with different spin orientations are scat- 
tered differently; experimentally this leads to a polarization of the scattered 
particles. An analysis of this effect can also be carried out within the framework 
of the optical model. 

There is one type of nuclear reaction which differs from the usual reaction 
since the nucleons involved exhibit a resonance effect. We are speaking of 
reactions induced by deuterons. 

The deuteron is distinguished from all other nuclei by its small binding energy. 
The binding energy of the deuteron is 2.23 Mev. Because of its small binding 
energy the deuteron is an “‘extended”’ system. The mean distance between the 
proton and neutron in the deuteron is approximately 4-10713 cm. This value 
exceeds the range of the interaction forces between nucleons; hence the deuteron 
(like a nucleon) cannot be considered a point particle. This situation leads to a 
number of characteristic features in deuteron reactions. 

Different processes can occur in the collision of a deuteron with a nucleus. 
In the first place a deuteron can be scattered elastically. The deuteron can also 
enter the nucleus as a whole, forming an intermediate nucleus. Both of these 
processes have no unusual features and will not be considered further. 

We consider reactions of the (d, n) or (d, p) type. In these reactions the out- 
going particle is part of the initial deuteron. These reactions can go via the 
usual channel—through an intermediate nucleus—but there is considerable 
probability (in many cases higher than the probability for formation of an inter- 
mediate nucleus) of another type of reaction. Because of its large size, the entire 
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deuteron may not enter the nucleus; it is sufficient if either the proton [(d, n) 
reaction | or the neutron [(d, p) reaction] does so. The second process has some 
unique features. 

In contrast to the deuteron as a whole, the neutron does not have to overcome 
Coulomb repulsion in order to penetrate the nucleus. Hence, at rather small 
deuteron energies (from several hundreds of kiloelectron volts to many millions 
of electron volts) the reaction being considered is more probable than the entry 
of the whole deuteron into the nucleus. A similar situation obtains in the (d, n) 
reaction. Although the Coulomb barrier for the proton is the same as for the 
deuteron, the proton mass is only 1/2 that of the deuteron and the probability for 
a proton to penetrate the barrier is greater than that of the deuteron. Deuteron 
reactions which take place without the formation of an intermediate nucleus are 
called Oppenheimer-Phillips reactions, from the names of the authors who first 
suggested (aithough not completely correctly) the theory of this effect. 

The interaction of a deuteron with a nucleus need not necessarily terminate 
in the capture of one of the nucleons; the interaction may simply lead to the 
splitting of a deuteron—the (d, np) reaction. This reaction can also take place 
outside the nucleus. The last reaction has been analyzed theoretically; how- 
ever, it remains almost entirely uninvestigated experimentally. 

At energies of several millions of electron volts the (d, p) and (d, n) reactions 
can be used to study characteristics of the levels of light nuclei. This possibility 
was first suggested by Butler* who also gave an analysis of the effect. Although 
the actual calculations carried out by Butler are unconvincing from a theoretical 
point of view and the theory of the effect has not been fully formulated, the 
qualitative results obtained by Butler are undoubtedly correct and are corrobo- 
rated by numerous experiments. It would seem that these results are highly in- 
sensitive to the actual assumptions used in the calculations. 

For example, consider the (d, p) reaction in a certain nucleus. From a 
knowledge of the energy of the proton after the reaction and the initial deuteron 
energy it is possible to determine at which level of the final nucleus the neutron is 
captured. If the nucleus is a light nucleus the spacing between levels is large 
and this identification is unique. Naturally this type of reaction can go via the 
intermediate nucleus; in this case, however, the protons are emitted uniformly 
in all directions. In the case in which only one neutron is “stuck,” the protons 
have a sharply defined angular distribution, the shape of which can be used to 
obtain information on the characteristics of the nuclear levels. The process in 
which the intermediate nucleus divides, which occurs in parallel with the above 
process, results in the appearance of an isotropic background which can be 
distinguished in principle. The angular distribution of the emitted protons is 
determined by the angular momentum carried into the nucleus by the neutron. 
If the neutron does not carry such momentum, the proton distribution has a 
narrow peak in the forward direction. This is the case which can most clearly 
be established in experiment. 

If the maximum in the angular distribution does not lie at zero angle, it means 


* S, T. Butler, Proc. Roy. Soc. (London) A 208, 559 (1951). 
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that the neutron has transferred to the nucleusa momentum Al> 0. The theory 
allows us, in a very rough way, to establish the position of the maximafor Al = 
1,2. Itis difficult to distinguish between higher values of Al. 

A large number of experiments have made it possible to obtain new (or to 
verify known) data on nuclear levels. There is considerable interest in the in- 
vestigation of successive (from nucleus to nucleus) nuclear levels in regions 
where shells have not been firmly established. It is clear that the same results 
can also be obtained in studies of (d, n) reactions in proton capture. 

The reaction which is the inverse of the one described above is also suitable 
for this purpose. Ifa proton enters a nucleus it can attract one of the neutrons 
in the nucleus, and a deuteron can be formed. As a result we have the (p, d) 
reaction. In the same way a neutron may eject a proton from the nucleus and 
cause the (n, d) reaction. These two reactions are called “pick up” reactions. 
These are the inverse reactions to the (d, p) and (d, n) reactions, and their proba- 
bilities can be expressed in terms of the probabilities for these latter reactions, 
using statistical relations (the principle of detailed balancing). Hence, all 
results which can be obtained from one reaction can also be obtained from the 
other. 

Similar reactions take place at high deuteron energies. In this region the 
(d, p) and (d, n) reactions are called stripping reactions. The nature of the 
effect is the same. However, the theoretical analysis is simpler. At small 
deuteron wavelengths the nuclear target can be considered a circular, opaque 
(or semi-opaque) object which absorbs all particles incident upon it—neutrons or 
protons. The stripping effect has been investigated in detail, both theoretically 
and experimentally. 

It is of interest to note certain features of (p, n) or (n, p) reactions in light nu- 
clei. These reactions are the inverse of each other (for a suitable choice of 
target), and hence their properties are similar. 

If, for example, a proton strikes a nucleus, then, as in the deuteron reaction, 
two reaction channels are possible. The reacton which proceeds via the inter- 
mediate nucleus will give neutrons which are distributed isotropically over 
angle; the other reaction is the direct ejection of a neutron from the nucleus. In 
this case conscrvation of energy requires (if the proton is not fast) that the in- 
coming proton stick in the nucleus. In this reaction the neutron angular dis- 
tribution depends on the proton state and the neutron state in the nucleus. 
The study of these reactions is of a special interest in nuclei in which these states 
are different since transitions between mirror nuclei lead to the trivial replace- 
ment of a proton by a neutron and vice versa—in this case the emitted neutron is 
essentially the “extension” of the proton, and the distribution exhibits a maxi- 
mum at zero degrees. 

Thus we see that the study of reactions which occur without the formation of 
an intermediate nucleus provides a unique method for studying the structure of 
the nucleus. 


LECTURE NINE 


w-Mesons 


Í n many cases the collision of a nucleon with a nucleus at high 
energies is characterized by the production of new particles—r-mesons. Proc- 
esses involving the production of r-mesons and the interaction of these particles 
with nucleons and nuclei have been intensely studied in recent years with the 
increasing availability of huge accelerators which can accelerate protons to 
energies greater than 10° ev. 

Before discussing the properties of r-mesons themselves, we shall establish 
certain general considerations which pertain to the production of particles and 
the transformation of particles from one type to another. 

Any process which takes place in a quantum system must conserve the total 
energy, the total angular momentum and the parity of the system. 

In computing the angular momentum of a system, in addition to the orbital 
momenta of the particles, we must take account of the inherent moments—the 
particle spins. In exactly the same way, in applying conservation of parity, we 
must consider not only the parity associated with the orbital momentum L, 
which is (—1)”, but also the (intrinsic) parity of the particles. Obviously, the 
intrinsic parity of particles can be considered only in cases in which the process 
does not result in the creation or annihilation of particles; in these cases the 
““orbital” parity is conserved separately. 

The intrinsic parity of particles is not uniquely defined. For example, it 1s 
possible, without violating any of the conservation laws, to change the parity for 
all particles in a system by multiplication by (—1)%, where Z is the charge 
expressed in units of electronic charge (by parity we are to understand the 
quantity +1 or —1, by which the wave function 1s multiplied in making a 
change from a right-handed coordinate system to a left-handed system). The 
total parity of all the particles before the reaction is multiplied by (~1)74, 
where the sum is taken over all particles. Similarly, the total parity of all the 
particles after the reaction is multiplied by (—1)*“". By virtue of charge con- 
servation EZ (before the reaction) is equal to ZZ’ (after the reaction); hence, 
this transformation does not violate the conservation law. 

There is, however, a case in which the parity of the particles can be given 
absolutely. We are referring to neutral particles, primarily photons. The 
parity of the m°-meson can also, in principle, be established from an investigation 
of the reaction in which it decays into two photons. 

In addition to the parity associated with coordinate reflection—space parity— 


87 


88 Lectures on Nuclear Theory 


we can consider the properties of a system under time reversal. In this case, in 
place of time reflection it is convenient to consider the simultaneous “reflection” 
of the space and time coordinates, rather than the time coordinates alone. 

We can classify particles in terms of behavior with respect to a change of sign 
for all coordinates (space and time). Exactly as in the case of space parity, we 
formulate a conservation law for the time parity of the system: the total time 
parity of the particles cannot be changed in any process. Here we are speaking 
of the total parity of the particles; in contrast with ordinary parity, there is no 
time parity associated with the orbital motion and we are discussing here only 
the conservation of the internal or intrinsic parity of the particles. 

The time parity of particles has not been studied to any great extent. 

We may say that the time parity of the photon is +1. This statement follows 
from the fact that the photon can appear without the annihilation or creation of 
any other particles in a system, 1. e., without a change in the time parity of these 
particles. The same considerations pertain to the m%-meson. 

It is of great interest to study the conservation law for time parity in order to 
understand the absence of a number of processes in nature which are formally 
possible. Such processes, for example, are those which lead to the annihilation 
of nucleons, say 

n+ p> rt +H 9” 


or 
p + e7 — 2y. 


These processes do not occur in nature; indeed, it is only because they do not 
occur that we have stable nuclei. 

It might be assumed that the hydrogen atom does not transform into two 
quanta (the last reaction) because the time parities of the electron and proton 
are different (similarly for the time parities of the proton and neutron). Then 
the hydrogen atom could not transform into an even system (2 quanta) since it, 
itself, would be an odd system. However, having prohibited this transformation 
on this basis we could not prohibit the decay of an even system—hydrogen 
molecule or deuterium atom—into a pair of photons. Thus, the simple ex- 
planation is unsatisfactory. 

Both parities which we have been discussing pertain to charged particles and 
neutral particles—nothing specifically related to the charge has been introduced 
in any of the above considerations. In neutral particles there is one other 
specific characteristic, related to the symmetry properties of the quantum- 
mechanical equations for neutral particles. 

Quantum-mechanical equations allow a unique transformation, called charge 
conjugation, in which the particle is replaced by the anti-particle. The anti- 
particle for the electron is the positron; the anti-particle for the r+-meson is the 
m~-meson. In October of 1955 the anti-proton was discovered by Wiegand, 
Ypsilantis, Segre and Chamberlain. The anti-neutron—a neutral particle with 
positive magnetic moment—has not yet been observed, but it is clear that it 
must exist. 


The property of the equations to which reference has been made is the fact 
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that the equation remains unchanged if we replace all particles of this system by 
their anti-particles and change the sign of the electromagnetic field. 

It is obvious when we speak of a particle—anti-particle pair, which of them we 
call the particle and which we call the anti-particle is really a question of 
notation and has no physical significance. 

We must be more specific in what is meant by a neutral particle. Strictly 
speaking, the neutron is not completely neutral—it has a magnetic moment. 
Hence, it is possible to have an anti-neutron—a particle with magnetic moment 
of the opposite sign. On the other hand, the photon is a neutral particle in the 
strict sense of the word. The 7°-meson is also a neutral particle. We may de- 
fine a neutral particle rigorously by saying that it is a particle which is identical 
with tis own anti-paritcle. 

According to this definition a particle may not be a true neutral particle even 
if it has no charge or magnetic moment. The simplest example of this type is 
the hydrogen atom in a state in which the total spin (electron—nucleus) is zero. 
Although the hydrogen atom in this state has neither total charge nor moment, 
it is still not identical with its own anti-particle—an anti-proton with a positron 
in the K-orbit. On the other hand, the photon has no anti-particle—it is 
identical with its own anti-particle and is a neutral particle in the true sense of 
the word. 

By this definition of a neutral particle we imply another property of neutral 
particles—charge parity. In making a transformation to the charge-conjugate 
system neutral particles are transformed into themselves. Hence (if we limit 
ourselves to particles with integral spin) the wave function for neutral particles 
may either change sign or not change sign. Thus two types of neutral particles 
arise—charge-odd and charge-even. 

The notion of charge parity has direct practical application. It is easily seen 
that the photon is a charge-odd particle. In making the transition to the 
charge-conjugate equation the field component changes sign. But the field 
components are the quantities which are associated with photons in quantum 
mechanics. Consequently, the wave function of the photon is characterized by 
sign as far as charge conjugation is concerned, and thus the photon is charge-odd. 

On the other hand, as will be seen below, the 7°-meson is charge-even. We 
may also note that a system of two photons or two mesons (or any two identical 
neutral particles) is always charge-even. 

Up to this point we have discussed only particles with integral spin. The case 
in which a particle has half-integral spin exhibits certain peculiarities. Particles 
of this type are described by special mathematical quantities called spinors. 
Whereas scalars, vectors and tensors come back to their original values when the 
coordinate system is rotated through 360°, a spinor reverses sign in such a 
transformation. This property of spinors does not lead to any inconsistency 
since the wave function itself is not measured in experiment; however, this 
property turns out to have an important effect on the reflection of wave func- 
tions. 

For example, consider reflection of the spatial coordinates. Suppose that 
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under a transformation of this kind the wave function is changed in accordance 
with some rule which we will write symbolically :* 


y => Iy. 


We repeat this transformation and return to the original coordinate system; 
the total operation can be written in the form: 


y > IUy) = Py. 


If we are dealing with a scalar, vector or tensor, we require that J? = 1, since 
the wave function must return to its original value. Whence it follows that two 
types of particles are possible, corresponding to the two signs for the root of 
unity. 

In the case of particles with half-integral spin the wave function, as has 
already been indicated, may either return to the initial value or may change sign. 
Correspondingly, there are two possible cases: either we can consider the 
double inversion as the identity operation, in which case J? = 1, or we can con- 
sider the case in which the wave function changes sign under double inversion, 
in which case J? = —1. Thus, in the case of spinors we have the possibility of 
four different inversion transformations and, correspondingly, the possibility of 
four “parities.” However, it is obvious that all four of these possibilities cannot 
be realized simultaneously. The square of the inversion operator must have a 
single definite value for all spinors; this, however, cannot be determined theo- 
retically. Each of the two possible values J? = +1 is consistent with the formal- 
ism of the theory. The determination of this fundamental property of the in- 
version operator can be made only on the basis of experiment. 

There is, however, an unexpected possibility for determining the value of J? 
from experiment. It turns out that if the square of the inversion operator J? = 1 
the wave functions of the particle and anti-particle transform differently under 
inversion; if, however, 1? = —1 the transformation laws are the same. This 
situation implies that a true neutral particle with half-integral spin can exist only 
if the second relation holds true. If the first law applies, obviously the particle 
and anti-particle cannot be identical since the transformation relations for their 
wave functions are different. 

Thus, to decide definitely which type of transformation to assign to the wave 
functions we must determine whether or not a true neutral particle exists in 
nature. A true neutral particle of this type could be the neutrino. The 
neutrino, by definition, has no charge; there is also a high probability that it 
has no magnetic moment. 

According to the usual 8-decay description, the neutrino (v) appears in the 
emission of a positron from the nucleus, while the anti-particle—the anti- 
neutrino ($) appears in the emission of an electron. If the neutrino is a true 
neutral particle the neutrino is identical with the anti-neutrino. (In this case, 


* In the case of a wave function which consists of one component (spin zero), J is 
simply a number. In the general case the different components of the wave function 
may be intermixed, and J, in general, is a matrix. 
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in particular, the magnetic moment of the neutrino must be identically zero.) 
There is a convincing method for deciding this question. This method is based 
on measurements of the lifetime for double B-decay. In principle, double £- 
decay takes place between two isobars which differ in atomic number by two 
units. The existence of a large number of such isobars in nature (with even 
numbers of protons and neutrons and spin zero) is due to the fact that either the 
intermediate nucleus (with odd numbers of protons and neutrons) has a large 
mass and ordinary ß-decay is impossible energetically, or it has high spin so that 
B-decay, while possible in a formal sense, is actually forbidden to such an extent 
that it becomes virtually unobservable. 
The usual mode of double 6-decay can be written schematically in the form 


2n—> 2p + 2e7 + v +20. 


If, however, the neutrino and anti-neutrino were identical, the decay of one 
neutron could occur with the emission of an anti-neutrino while the decay of the 
other neutron could occur with the emission of a neutrino or, what is essentially 
the same thing, with the absorpiton of the available anti-neutrino. Thus, the 
double @-decay process could take place without a neutrino at all: we would 
have the scheme 


2n —> 2p + 2e”. 


The theory indicates that the first process is approximately 10,000 times less 
probable than the second.* Since the second process (without neutrino) has a 
half-life greater than 101 years for real nuclei, the first mode of decay (anti- 
neutrino), in general, is unobservable with present-day experimental techniques. 
Hence, the observation of a double B-decay is possible only if it occurs according 
to the simple version in which the neutrino does not figure. 

Recently, McCarthy has suggested that he has observed double 6-decay 
(Ca® — Tis) with a half-life of approximately 10!” years. These results would 
mean that double ß-decay occurs according to the scheme in which the neutrino 
does not appear and that the neutrino is, correspondingly, a true neutral particle. 
Thus it would follow that the wave function of the neutrino transforms according 
to 1? = —1 under reflection of coordinates, i. e., its wave function changes sign 
upon double reflection. 

The considerations given above would have a wider meaning. It follows that 
all particles in nature would transform according to this scheme since we could 
not have particles which transform according to J? = 1 and particles which 
transform according to 7? = —1. The significant feature of these considerations 
is the fact that the transformation schemes for all spinors would be determined 
uniquely by the fact that there was even one true neutral particle. 

It is obvious, however, that if double B-decay could occur only with the emis- 
sion of the anti-neutrino it would still be impossible to draw any final conclu- 


* This is considered in greater detail in Zeldovich, Lukyanov and Smorodinsky, 
“Properties of the neutrino and double $-decay,’” Usp. Fiz. Nauk 54, 361 (1954). 
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sions, since this fact, in itself, would not tell whether or not there were true 
neutral particles in nature (even the neutrino). 

Having established the transformation scheme for the spinors (in particular 
I? = —1) we could now assert that all spinors can have two parities, correspond- 
ing to the two possible values of the root of —1. Whence it follows, for ex- 
ample, that a system consisting of two identical particles with half integral spin 
and without orbital moment is an odd system: the parity of such a system 1s 
determined by the products of the parities of the two particles; the reflection 
transformation will be the product of two identical transformations J, which as 
we now know is —1. If, however, the particles have different parity, the 
system is charge-even. In particular, a complex particle with integral spin 
which decays into two identical particles with half-integral spin and no orbital 
moment is necessarily an odd particle. 

Similar analyses can be carried out for time parity and charge parity (for 
neutral particles with spin !/2), but we will not stop here for this purpose. 

An interesting application of the theorem of charge parity occurs in the case of 
positronium—a bound system consisting of an electron and a positron, which is 
similar to the hydrogen atom. 

In accordance with our definition, positronium is a true neutral particle since 
under charge conjugation the positron becomes an electron while the electron 
becomes a positron; consequently, positronium is transformed into itself. The 
charge parity of positronium is determined by its orbital moment and spin. It 
can be shown that the charge parity of positronium is (— 1)” +5, where L is the 
orbital moment and S is the spin. Since the spin of positronium can be 0 
(singlet state) or 1 (triplet state), the parity of positronium is the same as the 
parity determined by Z in the singlet state and opposite to the parity determined 
by Z in the triplet state. 

We now consider the decay of positronium into photons. We have already 
indicated that the photon is charge-odd. This means that a system of an odd 
number of photons is charge-odd. Whence we have the following selection rule: 
positronium can decay into an even number of photons only if it is in the states 
15, 3P, 1D, etc.; it can decay into an odd number of photons in the states 4S, !P, 
3D, etc. 

The positronium states with the lowest energies are 1S and 38. These states 
are called para-positronium and ortho-positronium, respectively. We have seen 
that para-positronium (the 1S state) can decay only into an even number of 
photons, while ortho-positronium (the 35 state) can decay only into an odd 
number. Since decay to one quantum is impossible because of conservation of 
energy and momentum, we can conclude that the most probable decay schemes 
for positronium are 


150 e 27, 3S) =% 37. 


It follows that the decay of ortho-positronium will be less probable than the 
decay of para-positronium. The experimental work which has been carried out 
gives the half-lives for these decays as 7:1078 and 2-10" sec, respectively (in 
agreement with theory). 
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It is interesting to note the fact that the decay of ortho-positronium into two 
photons is forbidden as a consequence of only one theorem (namely, the 
theorem that a system of two photons cannot have a total moment equal to 
unity). 

All the general considerations which we have derived refer to any particle; 
in particular they refer to all mesons which have been recently discovered. 

In these lectures we are considering in detail only the properties of one kind of 
meson—the r-meson, since its properties have been investigated most com- 
pletely. This situation is a result of the fact that r-mesons are the only mesons 
of which intense beams have been obtained in the laboratory with accelerators.* 

T-mesons were originally discovered by Powell in cosmic rays and it was only 
later that they were produced artificially in large accelerators. Three types of 
T-mesons are produced in collisions of high-energy nucleons with protons or 
complex nuclei: charged (r* and m”) mesons and neutral (m?) mesons. The 
properties of all r-mesons are very similar. 

The similarity between r*- and r”-mesons is obvious from the production 
processes 


bpn + rt, nop+ rm. 


It is clear that transformations of this kind could not occur in a free nucleon, 
since they would violate conservation of energy and momentum. Processes of 
this type can occur only in the collision of a nucleon with another nucleon or 
nucleus. 

The 7°-meson is produced in a similar way: 


b>p+ T, n—on+ w. 


All four processes, as we shall see later, are connected by certain relations 
which appear as a consequence of the isotopic invariance of r-mesons; these 
relations allow us to consider all three r-mesons in terms of a single particle in 
three charge states (just as the proton and neutron are considered in terms of one 
particle). 

m-mesons are unstable particles. They decay in the free state. The decay 
schemes and lifetimes are different for charged r-mesons and neutral r-mesons. 
Charged r-mesons decay into a y-meson and a neutrino 


+ 
m ou + p. 


Since the mass of the u-meson is 206 me, the decay of the r-meson is characterized 
by the appearance of an energy of 35 Mev. Because of conservation of momen- 
tum the u-meson acquires 5 Mev while the remaining 30 Mev is carried away by 


* In recent years u-mesons have become available in the laboratory as products of 
the decay of r-mesons. However, the investigation of the properties of these parti- 
cles is just starting. In 1955, the Bevatron at Berkeley, which yields protons with 
energies of more than 6 Bev, was able to furnish a beam of K-mesons—particles with 
a mass equal to 970 electron masses. Only preliminary reports of these experiments 
have appeared so far. 
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the neutrino. The half-period of r -mesons is 
Try 7°) = (3.6 + 0.3)10°% sec. 
The neutral 7°-meson decays into two photons 
ey yy; 
the half-period for this process 1s 
T1/,( 7°) = (5 + 3)10~% sec. 


It should not be thought that the difference in decay modes and lifetimes 
violates the charge symmetry between m -mesons and 7°-mesons. The decay 
into two photons has no analog in the charged particle case. However, the 
process 7° — u? + v, if it exists (up to the present time there has been no indica- 
tion of the existence of the u°-meson), should have a probability approximately 10% 
times smaller than the probability of decay into two photons. 

The lifetime of the 7°-meson is itself short compared with the measured life- 
times of elementary particles. The path traversed by the 7°-meson in its life- 
time is of the order of 10” cm. This time, however, is large enough so that 
we can still speak of the -meson as an independent particle. The character- 
istic time of a particle (in order of magnitude) is h/uc?, since this is the single 
combination having the dimensions of time which can be formed from the atomic 
constants. In the case of the 7°-meson this time is 5:107? sec; in comparison 
with this quantity, the lifetime of the m°-meson is very large. Hence, the reality 
of the 7°-meson is as well founded as that of any other particle. 

With the exception of decay processes the similarities in the properties of the 
three mesons are very marked. Leaving out electric effects, which are important 
only in the study of slow mesons, the production of different r-mesons and the 
scattering of these particles on nucleons and nuclei are very similar. 

This similarity can be noticed immediately from the masses of r-mesons 


mn) = (272 + 0.3) electron masses, 
m(m) = (263.9 + 0.9) electron masses. 


Although the relative difference in mass is larger than in the case of the neutron 
and proton (approximately 0.15%), it is small compared with the mass of the 
meson itself. The similarity in the behavior of r-mesons at the present time 
has been taken to mean that we can apply the idea of charge invariance to r- 
mesons and consider them as parts of a triad. 

In terms of isotopic space the existence of three states means that the isotopic 
spin of the r-meson is 1 and that it can have three different projections on the 
g-axis: +1,0and —1. The three projections of the isotopic spin on the ¿-axis 
correspond to the three different r-mesons: 


tt = +1 (rt), re =0 (x), Tt = -—1 (rD), 


These conditions correspond to those used earlier in choosing the isotopic spin to 
be associated with the proton as +. Naturally, the isotopic invariance 
hypothesis for r-mesons requires experimental verification. 
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It is immediately obvious that considerations of charge symmetry can yield 
certain relations between the cross sections; for example, it may be asserted 
that in the collision of a proton and a neutron the number of charged mesons 
(of the same energy and in the same solid angle) of both signs must be equal 


olp + nn + n + rt) = p+ n> p+ p+r). 


Relations of this kind are essentially trivial. However, isotopic invariance leads 
to a more interesting relation. 

At the present time only one relation of this kind has been verified experi- 
mentally. There is a well-known reaction in which a deuteron and 7+-meson are 
produced in the collision of two protons 


p+pomarttd. 


We will encounter this reaction again later on. An analogous reaction is 
observed in the collision of a neutron with a proton: 


n+p> m +d. 


It can be shown that charge invariance implies that the cross section for the first 
of these reactions must be twice as large as the cross section for the second: 


olp + p—> rt +d) = 20(n + p > m + d). 


This relation is easily established if one notes that the final products of the 
reaction have a total isotopic spin equal to 1 (Tr = 1, Tg = 0). The isotopic 
spin of the two protons is 1, the isotopic spin of the n-p system is either 0 or 1 
(the projection Ty = 0). Both of these values are equally probable. Hence, 
the reaction in question is possible in only half the states of the n-p system. 
This relation has been verified experimentally. 

Interesting relations arise in the production of mesons in the collision of 
nucleons with nuclei whose isotopic spin is 0. Such nuclei are H?, He4, C?!?, 
and all nuclei with an even number of protons and neutrons (in the ground 
state). All three types of mesons can be produced in these collisions. The 
following relation obtains between the meson-production cross sections: 


at + o7 = 20%; 


this relation applies for a given meson energy and angle. The relation has not 
as yet been verified experimentally. 

A further problem is that of the spin of r-mesons. Direct measurements of 
the spin of these particles are extremely difficult. There is, however, a simple 
and extremely convincing method of indirect measurement. This method is 
based on an investigation of the relations between the reaction p + p => trt +d 
and its inverse—the capture of a positive meson by a deuteron. 

Using the principle of detailed balance we can establish a relation between 
these reactions which is independent of any assumption as to the mechanism 
involved. If we use the symbols p, and p, to denote the momentum of the 
proton and r-meson in the center-of-mass system, and the symbol s to denote the 
spin of the 7-meson 
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op + p — rt +d) = ?/{2s + 1)(pp/pr)?o(d + r* > p +p). 


Because of the factor (2s + 1), the ratio between these reactions will change 
markedly depending on whether the spin of the meson is zero (2s + 1 = 1) or 
unity (2s + 1 = 3). The experimental data definitely indicate that the spin 
is Zero. 

It can be independently established that, in any case, the spin of the neutral 
meson is not unity. This statement follows from the fact, mentioned above, 
that the 7°-meson decays into two photons; as we have already mentioned, a 
particle with spin 1 (for example, ortho-positronium) cannot decay into two 
photons. 

The next characteristic of the r-meson which must be established is its parity 
under coordinate reflection. The parity can be established from the reaction in 
which a negative meson is captured by a deuteron 


a .d>n>+m0m 


This reaction has a high probability and is easily observed. When a r”-meson 
enters a medium in which there are deuterium nuclei it is slowed down and 
captured by the deuteron in the K-orbit, forming a neutral system similar to 
hydrogen. Since capture of the r-meson takes place in the K-orbit, our analysis 
of the problem is simplified since it is known beforehand that the momentum of 
this system in the initial state is 1 (deuteron spin). Since the intrinsic parity of 
the deuteron is the same as the intrinsic parity of a system composed of two 
neutrons,* the parity of the system is equal to the intrinsic parity of the r-meson. 

We consider the final state of the system. A system of two nucleons can have 
a total moment equal to 1 in the four following states: 7S), ?P,, Pi, ?Dı. Of 
these states, only one is available to two identical nucleons, i. e., the $P, state. 
Thus, the parity of the system is determined uniquely. The state is odd; 
consequently, the r”-meson is an odd particle. 

It is reasonable to assume that the m*-meson is also an odd particle. In 
principle, there is a possibility for a direct measurement of the absolute parity 
of the 7°-meson, in contrast with the parity determination of the r=-mesons. 
The theory indicates that if the m°-meson is even it will decay into two photons, 
with parallel polarization (in the system in which the 7°-meson is at rest); 
if the m-meson is odd, the polarization planes for the two photons are 
mutually perpendicular. This experiment has not been carried out because 
of the serious experimental difficulties involved. 

Usually, in speaking of the parity of mesons with respect to coordinate reflec- 
tion, instead of the designation “odd r-meson” it is customary to say that the 
meson is a pseudoscalar particle, since an odd particle with spin zero is described 
by a one-component wave function which changes sign under coordinate reflec- 
tion and is therefore a pseudoscalar. 

It is also easy to establish the fact that r-mesons are even particles under time 





* We consider the intrinsic parity of the proton and neutron to be the same, 
thus determining the parity of the r-meson with respect to the p-n system. 
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reflection. This result is obvious directly from an examination of the processes 
in which these particles are formed in collisions between two nucleons. In such 
processes the number of nucleons remains unchanged; consequently, their time 
parity remains unchanged. Hence, the r-meson is even with respect to time 
reflection. 

i Finally, from the fact that the r’-meson decays into two photons we can 
conclude that the 7°-meson is an even particle under charge conjugation; that 
is to say, its wave function is not changed under charge conjugation. 


LECTURE TEN 


Interaction of w-Mesons with Nucleons 


T; problem which arises, once the basic properties of r- 
mesons have been established, is the examination of the laws which govern the 
interaction of these particles with nucleons. It turns out that the situation here 
is fundamentally different from that which obtains in electrodynamics. In 
particular, in contrast with the interaction of an electron with the electro- 
magnetic field, the interaction of mesons with nucleons is “strong” rather than 
“weak.” 

The strength of the interaction between the electron and the electromagnetic 
field is determined by the magnitude of the charge. We must not forget in 
discussing the numerical value of any physical quantity in theoretical physics 
that the value depends on which particular set of units is used to express the 
measurement. Hence, the notion of a “large quantity” or “small quantity” 
can be formulated only when we indicate with what the quantity in question is 
being compared; when this concept is applied to dimensionless quantities 
which are independent of a choice of units, a comparison is made with unity. 

Electric charge is a quantity with dimensions. However, we can form a 
dimensionless constant from the electric charge aud other constants which 
appear in the quantum mechanical equations: 

e? 
T he 
No other dimensionless constants can be formed from the parameters of a 
charged particle and the electromagnetic field; the constant a then can be taken 
as a measure of the strength of the interaction of the charged particle and the 
electromagnetic field. As is well known, its magnitude is small and is numeri- 
cally approximately */137. In many cases a still smaller quantity is important: 
e/2rhc ~ Viooo. Hence electrodynamics is an example of a weak-coupling 
theory. 

Because the quantity e?/Ac is so small, it is possible, using present-day per- 
turbation-theory techniques, to calculate any electrodynamic effect to any 
approximation. In recent years there have been brilliant demonstrations of the 
agreement between experiment and the theoretical calculations of the corrections 
to the energy levels in the hydrogen atom (the so-called Lamb shift) and the 
calculations of the intrinsic magnetic moment of the electron, which is eh/2mc 
in the first approximation; in the second approximation it 1s 
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eh (: m z) 
2mc Qn) 


In contrast with the weak coupling which is characteristic of electrodynamics, 
it turns out that meson interactions are strong. 

The clearest indication of the strong interaction appears in the production of 
m-mmesons by photons in nuclei (meson photoeffect). If heavy nuclei are bom- 
barded by photons with sufficiently high energy, the cross section for the produc- 
tion of r-mesons does not increase in proportion to the number of particles in 
the nucleus, i. e., in proportion to the atomic weight A; it goes approximately 
as A”. 

It is apparent that the probability for the production of mesons should be 
proportional to the number of particles in the nucleus—the interaction between 
the photon and nucleus is small and the probability that a photon is absorbed 
in the nucleus is small; hence, the photon intensity is uniform over the nucleus 
and the nucleons absorb independently. However, not all the mesons which are 
produced in the nucleus can escape. Because of the strong meson interaction 
there is a high probability that mesons are absorbed by one of the nucleons in 
the nucleus; only a fraction of the mesons escape: namely, those formed near 
the surface of the nucleus. It is obvious that in this case the number of mesons 
should be proportional to the nuclear surface area, i. e., proportional to 4”. 
Thus, the observed dependence of the meson photoeffect on atomic weight is an 
indication of the small meson range in nuclear matter; in turn, this fact is 
evidence of the existence of a strong interaction. 

Perturbation theory cannot be applied in the case of a strong interaction. On 
the other hand, the practical approach to calculating these effects is based on 
perturbation theory. In addition, Pomeranchuk has shown that all forms of the 
interaction which have been considered up to now are useless when examined in 
detail. 

Thus, at the present time we have no theory for meson interactions; moreover 
it is clear that the formulation of such a theory will be impossible without radical 
modification of some of the fundamental concepts. 

Notwithstanding the absence of any quantitative meson theory there is no 
doubt that r-mesons play an important role in nuclear interactions and are 
responsible for many of the properties of nucleons. 

In our first lecture we have noted that the magnetic moment of the proton is 
2.79 (in nuclear magnetons), while the magnetic moment of the neutron is —1.91. 
According to the Dirac equation, however, the magnetic moments of these 
nucleons should be 1 and 0, respectively. 

Hence, the proton has an excess positive magnetic moment of 1.79, while the 
neutron has an excess negative moment of 1.91. The fact that these two num- 
bers are so close in value is far from accidental. There is little doubt that this 
additional magnetic moment results from the interaction of the nucleon with 
the meson field. Charge invariance considerations allow us to state that the 
magnetic moments associated with the mesons should be the same in magnitude 
and opposite in sign in these nucleons since the proton and r”-meson form a 
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pair similar to the neutron and r*-meson. However, the lack of any theory 
makes it impossible to give a quantitative analysis of this effect. 

It is also undoubtedly true that r-mesons play an important role in nuclear 
forces. From the point of view of contemporary theoretical physics the inter- 
action of protons and neutrons can be described in terms of meson exchange. If, 
for example, we have a proton and neutron, the proton can emit a r*-meson, 
thereby becoming a neutron, while the emitted meson can be absorbed by the 
neutron, which is then transformed into a proton. This meson exchange auto- 
matically leads to nucleon interaction, since radiation and absorption (even 
virtual radiation and absorption) are accompanied by a change in the nucleon 
momentum. In the exchange of charged mesons the charge of the particle is 
changed; as a result the proton becomes a neutron and the neutron becomes a 
proton, and the interaction is of exchange character. On the other hand, the 
exchange of neutral mesons leads to the usual (nonexchange) forces; thus the 
meson mechanism provides for both types of forces. 

To estimate the range of the forces which are to be associated with meson 
exchange we must return to dimensionality considerations. The range can 
depend only on the meson mass. The mass of the nucleon does not appear in 
an important way in the range expression since this range must be finite for a 
nucleon at rest (infinitely heavy nucleon). Using the mass of the r-meson (pu) 
and the constants % and ¢ we can form only one constant which has the dimen- 
sions of length—the Compton wavelength of the meson* 


ro m a = 1.4. 10723 cm. 
uc 
This quantity then determines the range of the interaction. The dimensional 
estimate, as we have seen, agrees with the range of nuclear forces (~107! cm), 
thereby tending to support the mesonic origin of these forces. 

It is possible that the heavier mesons which have been discovered in recent 
years may also make a contribution to nuclear forces. However, the experi- 
mental data on these particles are as yet very incomplete, so that it is impossible 
to say anything more definite about them. 

We now consider the existing experimental data on the interaction of r- 
mesons with nucleons. Of the enormous amount of experimental data which is 
available we shall be concerned only with those effects which allow us to draw 
direct conclusions as to properties of r-mesons. 

We start with the simplest process: the scattering of r-mesons by nucleons. 
This process has been studied over a wide meson-energy range: from 20 to 1500 
Mev; research has been done on the scattering of both +* and +~-mesons by 
protons. Meson charge exchange (m7 + p— =’ + n) has also been studied up 
to several hundreds of millions of electron volts. 

The energy dependence of the cross section for these reactions is shown in 
Fig. 9. The cross section for the scattering of r+ on protons (equal to the cross 


* If the mass of the meson were zero (as is that of the photon) it would be im- 
possible to form a constant with the dimensions of length. In this case the forces 
would have an infinite range (like Coulomb forces). 
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Fig. 9. Total cross sections for the interaction of positive and negative r- 
mesons with protons (the figure is from an unpublished survey by Barkov 
and Nikolsky). The interaction curves in the states with isotopic spin 
T = 1/,are obtained from the experimental data using the formula 20,(T = 
1/2) = 30 r7) — or); Ez is the energy of the z-mesons in the laboratory 
system. 


section for r” on neutrons) reaches a maximum of approximately 200 mbarns at 
an energy of 200 Mev. The total cross section for m” scattering (scattering and 
charge exchange) on protons also reaches a maximum in this region, but its 
magnitude is smaller, approximately 65 mbarns; approximately 1/3 of this cross 
section is due to elastic scattering. As the energy is increased the cross section 
falls off rapidly (at an energy of 360 Mev, o:(4*, p) = 50 mbarns) while at the 
highest energies which have been achieved, 1500 Mev, the cross sections are 
approximately 35 mbarns for o,(#”, p) and 30 mbarns for omt, p). 

It is interesting to note that if only the state with isotopic spin 3/2 participates 
in scattering* the scattering cross sections should obey the following relation: 


olp, rt — p, wt): o(p, t7 — p, w):6(p, 77 > p, 77) = 9:2:1. 


We see that this relation is approximately satisfied in the region of the maximum 
(at higher energies the cross sections are almost the same). Thus, in this region 
the interaction 1s a strong function of the isotopic spin—a property which we 
have already noted in neutron interactions. 

The pattern seems to be associated with scattering in the “level” with T = 
3/2, and the phase analysis indicates that the corresponding state is ps,,. 

It 1s extremely interesting to examine the cross section for the meson photo- 
effect in the proton. This effect is also characterized by a maximum cross sec- 
tion at a photon energy of approximately 200 Mev, indicating a strong inter- 


* Obviously, the r-meson and the nucleon can have a total isotopic spin of either 
1/2 or 3/2, 
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action in the T = ®/, state; unfortunately the absence of any theory prevents 
our carrying out a more thorough analysis of the experimental data. 

Interesting effects are observed when mesons are produced as a result of 
nucleon collisions. In proton collisions, both r+- and m%-mesons can be produced 
in principle: 


pPtrpmptntaet, 
PPPS Pea m. 


Actually, however, the cross section for the production of r'-mesons is very 
small and there are virtually no 7°-mesons near threshold. Thus, at an energy 
of 460 Mev the cross section is divided between three possible reactions, as 
follows: 


elastic (p, $) 20 mbarns 
rr * production 4.5 mbarns 
m’ production 0.4 mbarns 


Total cross section ~ 25 mbarns 


Even at 660 Mev the fraction of +°-mesons produced in p-p collisions amounts 
to only 10% of the total cross section. 

The small probability for the production of m’-mesons can be explained if we 
assume that mesons are produced in singlet states of the p-p system. Since 
there is little excess energy close to threshold, the final state will most probably 
have L = 0 (35, or 150); the two protons can occupy only the 4S) state, but both 
states are available to a proton and neutron. It is easy to see that the rcaction 
is forbidden for two protons in a final state IS, if the initial states are singlet 
states. The singlet states of the p-p system (150, D2) have even momenta and 
are even. Since the final state of the p-p system has zero moment and is even, 
all the initial momentum is transferred to the meson. But the meson is odd; 
hence, having even moment it must be in an odd state. Consequently, the 
process is forbidden by conservation laws. 

There is one further restriction on the production of mesons in neutron-proton 
collisions. The following reactions are possible: 


»t pr ey a ae 
n+p>n+p+ ow. 


The first reaction has a rather small probability. This effect can be explained 
in the same way if it is assumed that mesons are produced only in singlet states 
of the colliding nucleons. 

As the energy increases the meson-production cross section also increases. 
At 850 Mev the r-meson production cross section in proton-proton interactions 
reaches a value of approximately 25 mbarns and remains more or less constant 
up to 1275 Mev. As we have already indicated, this cross section coincides with 
the elastic cross section and corresponds to the interaction of two absolutely 
absorbing spheres—protons with a radius of 4.5:1071 cm. 

The extent of the energy region for which this simple model applies is as yet 
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unknown. Itis interesting that r-meson scattering on protons at 1400 Mev can 
also be described by a model in which the proton is considered a “black” sphere 
with the same radius.* 

There is one more effect, first observed by Migdal,f which is observed in meson 
production in nucleon-nucleon collisions near threshold. The meson spectrum 
indicates that the mesons carry away a large part of the energy so that the 
fraction left to the nucleons is relatively small. This energy distribution is 
related to the nature of the nucleon interaction. We have seen that this inter- 
action exhibits a resonance effect at low energies; it follows, as can be demon- 
strated rigorously, that the nucleons have high probability of having small 
relative velocities (strong interaction). 

This effect is so strong that deuteron formation is very probable in the final 
state; at these high energies this effect is very strange. At an energy of 460 Mev 
a deuteron is formed in approximately !/; of the cases in which r*-mesons are 
produced (in proton-proton collisions); at 660 Mev the number of deuterons is 
still sizable. Deuteron formation due to resonance effects is observed not only 
in meson production but also in other nuclear reactions. 

As the nucleon energy available to experimenters has increased, new effects 
have been observed. At the present time experiments are reported in which 
studies have been made of events in which several r-mesons are produced 
in meson-meson interactions. However, as yet these are very rare and an 
analysis would not be fruitful. 

New effects are observed at still higher energies; as the energy is increased the 
probability of multiple-meson production becomes greater than the probability 
of single-meson production. 

When nucleons with energies greater than 10,000 Mev collide with nuclei, 
nuclear disintegration events (‘‘stars’’) occur in which the majority of particles 
are mesons. The cross section for star production is comparable with the 
geometric dimension of the nucleus— TR? (R is the nuclear radius). 

This property of the interaction distinguishes the “strong” nuclear interaction 
from the “weak” electromagnetic interaction. The formation of two photons of 
comparable energy is always less probable (by a factor of a) than the formation 
of a single photon;{ the cross section for the formation of three photons is still 
less probable. No matter how high the energy, there is never multiple produc- 
tion of high-energy photons solely as a result of an electromagnetic interaction.$ 
On the other hand, the intensity of the nuclear interaction increases with energy, 
as does the number of r-mesons produced in collisions. 


* W. B. Fowler, R. M. Lee, W. D. Shephard, R. P. Shutt, A. M. Thorndike and 
A. Winston, Phys. Rev. 97, 809 (1955). 

TA. B. Migdal, J. Exptl.-Theoret. Phys. 28, 3 (1955). 

I It should be emphasized, however, that this statement does not apply to low- 
energy photons. As is indicated by the theory, any process is accompanied by the 
radiation of photons with very low energy; these do not play any significant role in 
the majority of cases. This effect is called the infrared catastrophe. 

$ In stars many photons arise as a result of m-meson decay; the production of 
these mesons is of multiple character. 
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In the absence of any meson theory we cannot offer any analysis of multiple- 
meson production; the extent to which the single-meson production mechanism 
applies in this case is also not clear. 

The fact that the interaction increases with energy, however, does allow us to 
formulate a quantitative theory for the effects in the high energy region, in 
Which the number of mesons becomes so large that thermodynamic consider- 
ations can be applied. The idea of applying thermodynamics in multiple pro- 
duction of mesons was first suggested by Fermi.* The method proposed by 
Fermi is based on the assumption that in the collision of high-energy nucleons all 
the energy appears instantaneously in a small volume and that both nucleons 
come to rest in the center-of-mass system. 

Because of the Lorentz contraction this volume is not spherical but is com- 
pressed in the direction of motion of the nucleons prior to the collision. The 
amount of compression is given by the well-known factor 


1 v Me? 
È E” 


where M is the mass of the nucleon and E” its energy in the center-of-mass sys- 
tem. The latter is related to the energy of the incoming nucleon E by the 
relation 








2E” 
— Me?" 


The volume of the region in which energy appears is (order of magnitude) 


Mc? ( 3 
== bs 
E” Nuc 


since h/uc (u is the meson mass) is of the order of the interaction range. We 
assume that the interaction is so strong that there is time for statistical equi- 
librium to be established during the collision; thus, the entire system can be 
characterized by a temperature To (which, of course, is very high). 

If To is so large that 


II 


To > pe?, 


many mesons can be produced in the system. However, the conditions under 
which we may speak of a definite number of mesons in a system are not clear; 
at high densities, when the dimensions of the system and the interaction range 
are of the same order of magnitude the notion of what we mean by individual 
particles, as far as the principles of quantum mechanics are concerned, is not 
obvious. Hence, one does not expect the number of mesons observed in a, "meson 
star to be determined by the equilibrium conditions in the initial stages ‘of the 
collision. 

* A translation of the Fermi paper is given in Usp. Fiz. Nauk 46, 71 (1952). A 


detailed presentation of these problems is given in a paper by S. Z. Belenky and L. D. 
Landau, Usp. Fiz. Nauk 56, 309 (1955). 
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During the entire time in which the system is in a highly compressed state, 
new particles are being produced and old ones are being absorbed. These 
processes come to a halt only when the system has expanded to such an extent 
that the interactions become small. The number of particles which remain in 
the system at this time will be different from the number of particles at the initial 
stage of the collision and will be more or less determined by the expansion process. 
Hence, it is obvious that thermodynamic considerations alone are not sufficient 
for a description of the system; we must also be concerned with the hydro- 
dynamics involved in the expansion of a system in thermal equilibrium. This 
problem, however, is not analogous to the classical hydrodynamic problem since 
the velocity of the particles is very high and relativistic equations must be used. 

There are certain general considerations, derived from relativity theory, which 
allow us to draw conclusions as to the dependence of the number of particles on 
the energy of the initial incoming particle. In this problem an important part is 
played by the equation of state of the system—the relation between the energy 
density and the pressure. This equation cannot be introduced in general form— 
it is determined by the actual properties of the system and is different for differ- 
ent systems. However, using relativistic mechanics, we can indicate a certain 
limiting form for the equation of state. Specifically, it can be shown that the 
energy density 6 and the pressure p must satisfy the inequality 


p S 8/3. 


Asis known, in the case of light the relation p = 8/3 applies; the same relation 
is valid for an ideal relativistic gas (i. e., a gas of noninteracting particles moving 
with high velocity). It may be assumed that in the case of a highly interacting 
gas the pressure is close to the limiting value; hence, we shall postulate that the 
equation of state for the system being considered is of the form 


p = 8/3. 


Now, in considering the rapid dispersion of the particles we must take account 
of the fact that the temperature of the system is reduced as the particles fly 
apart. As the temperature decreases the density of the particles decreases, and 
there is a resultant reduction in interaction, 1. e., the probability of creation and 
annihilation of particles is also decreased. It is clear that the final number of 
outgoing particles is determined when the temperature of the system falls to a 
value such that a further change in the number of particles becomes very im- 
probable. In order-of-magnitude terms this temperature Ty is determined by 
the relation 


Ty ~N pc, 


The temperature 7', should be essentially independent of the properties of the 
system because the equilibrium density of the particles is a rapidly varying 
function of temperature. During the period of hydrodynamic expansion, 
because of the strong interaction between particles the path lengths remain small 
so that we can neglect viscosity and thermal-conduction effects and consider the 
hydrodynamic stage of the expansion as an adiabatic process. At the same 
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time, the ratio of the density of particles to the entropy density of the system is a 
very slowly varying function of temperature. Since the temperature is approxi- 
mately constant during expansion we may consider this ratio as constant and 
independent of the properties of the system. Consequently, the entropy 
density and the density of the number of particles are proportional: 


n = const s, 


where s is in dimensionless units; the constant, from dimensionality consider- 
ations, must be of the order of (uc/h)*. These relations also allow us to deter- 
mine the approximate dependence of the number of particles on the energy of the 
incoming nucleon. 

Since we have assumed that the pressure is 1/3 of the energy density, it follows 
from well-known thermodynamic relations that 


smg”, 
This can be shown most simply by an analogy with black-body radiation, in 


which the same law applies (6 = 3p). In the case of a blackbody € ~ T‘, and 
consequently 


Tre’; 
on the other hand, for a unit volume 
ds = Tds, 


whence it follows that 8“=s. 

As we have indicated, the volume of the system 2 is equal to the volume of a 
sphere of radius a ~ (h/uc), reduced by the factor Mc?*/E”. All the energy 
which appears in the collision, i. e., 2E’, appears in the volume of this ellipsoid. 
Consequently, the energy density is 


2E” 2E' 
— SS SS a EN wh" 
Q 4/zra?( Mc?/E’) 
Therefore, the entropy density is 
sE mE"? 


(E' is the energy in the center-of-mass system). Since the density of particles is 
proportional to the entropy density, 


an E"?, 


and the total number of particles is 
1 
N = no~ E", = E"? 


since 2 is inversely proportional to the energy. 


In order to transform to energy in the laboratory system, E, we must sub- 
stitute (E^)? = 1/¿Mc?E. Finally we obtain 
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N~E™, 


This expression coincides with the relation obtained by Fermi (although he 
started from the incorrect assumption that N is the number of particles created 
at the instant of collision). 

A noteworthy feature of this formula is the fact that it predicts a slow growth 
in creation multiplicity with energy of the incoming particle. This result ex- 
plains the well-known experimental fact that stars with more than a hundred 
particles are observed very rarely and that the majority of observed stars have a 
comparatively small number of high-energy particles (the mean particle energy 
Emean™E/N~E™"), 

Strictly speaking, collisions which are observed in cosmic rays are not colli- 
sions between two nucleons, but collisions between nucleons and nuclei or be- 
tween nuclei and nuclei. The picture for a collision between a nucleon and a 
nucleus is virtually the same as that for a collision between two nucleons, except 
that the nucleon has a radius much smaller than the nucleus, so that it strikes 
only those nucleons in the nucleus which lie in its path. Hence, the number of 
particles created in such collision will be essentially the same as the number 
created in collisions between two nucleons, and will be relatively independent of 
the atomic weight of the nucleus. 

If, however, both nuclei are heavy, the number of created particles will be a 
strong function of the dimensions of the nucleus. In this case the number of 
created particles is proportional to the atomic weight raised to the three-fourths 
power: 


NWEYıa 


Thus, a nucleus turns out to be much more effective, in terms of creation of 
new particles, than a proton with the same energy. 

We may also consider the angular distribution of the created particles. Fermi 
assumed that the particles are emitted isotropically in the center-of-mass system. 
However, there is no basis for this assumption. A careful analysis of the angular 
distribution requires a complete solution of the rather complicated hydro- 
dynamic problem; qualitative considerations alone are not very useful. 

This solution can, in fact, be obtained. It leads to the following picture of 
the angular distribution. In the center-of-mass system the majority of particles 
are emitted at large angles with respect to the direction of the incoming particle. 
However, the density of particles at small angles is larger, so that the angular 
distribution curve is highly anisotropic. In the laboratory system the distribu- 
tion also has a peak at angles close to 0°. 

This type of angular distribution leads to a peculiar energy distribution. In 
this distribution a comparatively large number of particles have energies higher 
than the mean energy and a large number of particles have energies considerably 
lower than the mean value. 

A comprehensive comparison of the statistical theory of multiple production 
of r-mesons with experiments cannot be carried out since only a few stars have 
been recorded, and these are not suitable for analysis. 


